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Executive  Summary  of  "Vibration  Dynamics  and  Control  of  Bladed  Disk  Assemblies" 

This  final  report  documents  the  work  performed  at  Purdue  University  during  the  period  of 
November  1988  to  December  1990.  The  original  APOSR  Contracts  (#AFOSR- 89-0002, 
AFOSR-89-0014)  were  written  for  two  years.  Consequently  this  research  v. as  partially  funded 
from  Professor  Nwokahs’  PRF  grant  #670-1667.  The  objective  of  the  proposed  research  was  to 
gain  a  fundamental  understanding  of  how  and  why  periodically  configured  mechanical  and 
structural  systems,  (in  particular  bladed-disk  assemblies)  with  cyclic  symmetry  and  nominally 
identical  sub-structures  can  display  non-uniform  amplitudes  of  vibration  when  subjected  to  small 
but  random  parameter  perturbations  that  are  often  within  the  component  manufacturing 
tolerances.  A  secondary  aim  of  the  proposal  was  to  determine  ways  of  passively/actively  (if 
possible)  controlling  these  uneven  vibration  amplitudes.  This  work  specifically  dealt  with  the 
influence  of  the  double  (degenerate)  eigenvalues  present  in  every  cyclic  mechanical  system  and 
their  subsequent  splitting  under  small  perturbations,  on  the  uneven  vibration  amplitudes  of  the 
components. 

Status: 

The  work  associated  with  the  principal  objectives  of  the  project  is  almost  completed  and  is 
included  in  this  final  report.  The  procedure  for  detecting  a  priori  which  degenerate  eigenvalue 
pairs  will  split  under  given  parameter  perturbations  has  been  formalized  by  use  of  finite  group 
representation  theory  and  is  presented  in  Appendix  1.  The  procedure  for  accurately  unfo’ding 
the  singularities  induced  by  the  splitting  of  the  double  modes  of  cyclic  systems  has  been 
formalized  by  the  use  of  a  singular  perturbation  analysis  technique  which  is  valid  for  any  finite 
order  cyclic  system  and  is  included  in  Appendix  2. 
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The  topological  basis  for  the  singularities  induced  by  the  double  modes  and  the 
consequences  there  of  are  carefully  examined  and  detailed  in  Appendix  3. 

In  contra-distinction  from  recent  work  in  the  bladed-disk  research  literature,  numerical 
studies  which  show  that  uneven  amplitudes  of  vibration  in  perturbed  cyclic  systems  can  arise 
both  under  strong  coupling  as  well  as  the  weak  coupling  conditions  is  included  in  Appendix  4. 

A  systematic  framework  has  now  been  established  for  a  detailed  study  of  perturbed  cyclic 
systems.  Future  efforts  will  be  aimed  at  completing  any  remaining  theoretical  analysis, 
development  of  computational  algorithms  for  such  analysis,  passive  structural  redesign  to  avoid 
localized  high  vibration  amplitudes  and  experimental  validation  of  the  analytical  results. 

Publications 

Six  papers  have  been  developed  from  this  work.  The  first  is  essentially  Appendix  1.  The 
second  and  third  are  included  in  Appendix  2.  The  fourth  is  given  in  Appendix  3,  while  the  rest 
are  in  conference  proceedings  as  given  below. 

1.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  On  the  dynamics  of  perturbed  symmetric 
systems.  Accepted  for  presentation  and  publication  in  conference  proceedings  for  13th 
Biennial  ASME  Conference  on  Mechanical  Vibration  and  Noise,  September  22-25,  1991, 
Miami,  Florida. 

2.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  A  singular  perturbation  perspective  on 
mode  localization.  J.  Sound  and  Vibration,  (To  appear). 

3.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  A  singular  perturbation  analysis  of 
eigenvalue  veering  and  mode  localization  in  linear  chain  and  cyclic  systems. 
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J.  Sound  and  Vibration,  (Submitted). 

4.  Nwokah,  O.D.I.,  Afolabi,  D.,  Damra,  F.M.,  On  the  modal  stability  of  imperfect  cyclic 
systems.  Control  and  Dynamic  Systems,  35,  137-164,  1990. 

5.  Afolabi,  D.,  Nwokah,  O.D.I.,  The  frequency  response  of  mistuned  periodic  systems.  Proc. 
12th  ASME  Biennial  Vibrations  Conference,  Montreal,  Canada,  September  i9*T 

6.  Afolabi,  D.,  Nwokah,  O.D.I.,  Effects  of  mild  perturbations  on  the  dynamics  of  structures 
with  circulant  matrices.  Proc.  ALAA/ASME/ASCE/AHS  Structures,  Structural  Dynamics 
and  Materials  Conference,  Long  Beach,  CA,  1990. 

Personnel 

Three  faculty  members  and  one  graduate  student  were  funded  by  this  contract. 

Professor  Osita  D.I.  Nwokah  (Resume  at  the  back). 

Professor  Anil  K.  Bajaj  (Resume  at  the  back). 

Professor  Dare  Afolabi  (Resume  at  the  back). 

Gemunu  S.  Happawana  (Doctoral  Student). 

A  doctorate  degree  is  expected  to  be  awarded  for  this  work  in  the  next  24  months. 

Presentations 

Several  seminar  presentations  resulted  from  this  work.  Details  are  contained  in  the 


individual  Professors’  resumes. 
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1.  INTRODUCTION 

1.1  Problem  Statement 

The  central  aim  of  the  present  project  has  been  to: 

(i)  Gain  a  fundamental  understanding  of  how  and  why  periodically  configured  mechanical 
and  structural  systems  with  cyclic  symmetry  and  nominally  identical  sub-structures  can 
display  non-uniform  amplitudes  of  vibration  under  differential  (i.e.,  small)  parameter 
perturbations  that  are  often  within  the  component  manufacturing  tolerances. 

(ii)  Design  passive  and/or  active  control  mechanisms  to  overcome  such  possible  uneven 
amplitudes  of  vibrations. 

1.2  Background  and  Overview 

The  study  of  cyclically  configured  dynamical  systems,  otherwise  known  as  bladed-disk 
assemblies,  has  been  a  very  active  area  of  research  in  structural  dynamics  over  the  last  25  years. 
It  is  a  measure  of  the  theoretical  difficulties  involved  in  an  accurate  analysis  that  at  the  present 
time  there  is  no  general  agreement  in  the  literature  as  to  either  the  causes  of  the  uneven 
component  vibration  amplitudes  or  as  to  which  component  will  vibrate  with  the  highest 
amplitude  under  parameter  variations.  This  is  a  sine'  qua-non  to  establishing  benchmark 
specifications  for  component  vibration  control.  The  work  performed  under  this  grant  in  the  last 
two  years  clearly  indicates  that: 

(i)  Uneven  amplitudes  of  vibration  are  caused  by  the  modal  bifurcation  phenomenon  or  the 
sensitive  dependence  of  eigenvectors  on  small  parameter  variations  under  some  clearly 


defined  conditions.  [1,2,31 
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(ii)  Modal  bifurcations  in  turn  are  caused  by  mistuning  or  small  parameter  variations  from 
nominal  design  values,  which  are  often  within  the  component  manufacturing  tolerances. 
[4] 

(iii)  Extreme  cases  of  the  uneven  amplitudes  of  component  vibrations  produce  the  mode 
localization  phenomenon. 

We  had  shown  in  the  paper  in  Appendix  3,  that  very  useful  qualitative  information  on 
the  blade  mistuning  could  be  obtained  by  application  of  the  methodology  of  singularity 
theory  to  this  problem. 

To  understand  mode  localization,  one  must  first  study  modal  bifurcations. 

Let  )i  e  T  c  Rr,  where  Rr  is  an  r-dimensional  parameter  space.  If  a  given  structural  system 
has  n  degrees  of  freedom,  then  the  characteristic  equation  for  natural  frequencies,  to2  =  X,  can  be 
written  as  the  n-th  order  polynomial  equation: 

G(X  ,  p.)  =  0 ,  for  some  p  e  T  . 

It  turns  out  that  the  characteristic  polynomial  under  appropriate  modifications  behaves  like  the 
potential  function  in  singularity  (or  catastrophe)  theory.  [5]  Hence  the  degenerate  critical  points 

of  G(\  ,  p)  =  0  correspond  to  the  repeated  roots  (i.e.,  the  repeated  eigenvalues)  of  G(X  ,  p).  By 

5G 

studying  -  together  with  appropriate  higher  order  differentials,  and  G(^. ,  p)  =  0,  we  can 
oh 

determine  the  set  of  all  pi  e  T  at  which  G(X  ,  p)  has  degenerate  eigenvalues.  This  set,  which  is 
called  the  bifurcation  (or  catastrophe)  set,  partitions  the  parameter  space  T  into  distinct 
submanifolds  whose  boundary  is  the  bifurcation  set. 

We  can  conclude  from  the  basic  theorems  and  results  from  singularity  theory  and  catastrophe 
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theory  [6,7,8]  that  the  modal  behavior  of  our  structural  system  displays  sensitive  dependence  on 
parameters  only  in  the  neighborhood  of  the  bifurcation  set.  We  have  identified  two  distinct 
degenerate  behavior  patterns  in  structural  systems,  namely: 

a)  Coupling  induced  degeneracy, 

b)  Geometry  or  symmetry  induced  degeneracy. 

Furthermore,  we  have  noted  that  one  dimensional  lattice  type  periodic  structures  need  to  be 
divided  into  two  main  classes: 

(i)  The  Linear  Chain, 

(ii)  The  Cyclic  Chain. 

Each  of  these  classes  has  its  own  peculiar  characteristics  which  are  dictated  both  by  the 
geometry  (boundary  conditions)  and  the  physics  of  the  system.  For  example,  in  the  linear  chain, 
degenerate  and  therefore  ‘seemingly’  unpredictable  behavior  under  perturbations  appears  to 
occur  only  under  very  weak  coupling  conditions.  Topologically,  this  behavior  is  equivalent  to 
an  unfolding  of  the  m-fold  (here  m  is  the  number  of  nominally  identical  subsystems  which  arc 
weakly  coupled)  degeneracy:  (X  -  co2)171  =  0.  This  corresponds  to  what  Pierre  has,  in  a  series  of 
papers,  consistently  referred  to  as  a  perturbation  of  the  uncoupled  system  behavior.  [9,10,11] 
This  behavior  does  not  exist  under  strong  coupling  conditions. 

On  the  other  hand  for  cyclic  systems,  even  under  very  strong  coupling  conditions,  extra 
degeneracy  is  induced  by  the  cyclic  symmetric  nature  of  the  system  matrices.  It  is  then  well 
known  that  cyclic  systems  have  several  pairs  of  degenerate  (coincident)  eigenvalues  which  is 
distinct  from  the  case  of  linear  chains  where  no  degeneracy  or  multiplicity  of  eigenvalues  arises 
under  strong  coupling  conditions.  The  crucial  observation  is  that  because  of  the  coincidence  of 
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eigenvalues,  and  the  continuity  of  eigenvalues  with  respect  to  parameters,  cyclic  systems  (to 
which  bladed-disk  assemblies  belong)  always  operate  in  the  neighborhood  of  the  bifurcation  set. 
For  cyclic  systems  it  is  therefore  of  great  interest  to  determine  the  relative  influence  of  coupling 
and  geometry  in  the  subsequent  degenerate  system  behavior.  Since  a  tuned  bladed-disk  assembly 
has  pairs  of  degenerate  eigenvalues,  the  parameters  corresponding  to  the  tuned  state  are  clearly  a 
subset  of  the  bifurcation  set. 

The  number  of  degenerate  pairs  of  eigenvalues  as  well  as  the  effect  of  different  types  of 
perturbations  depends  on  the  nature  of  the  symmetry.  Some  of  the  qualitative  ramifications  of 
the  geometric  symmetry  can  be  studied  using  the  theory  of  groups.  The  effects  of  symmetry 
preserving  and  symmetry  breaking  perturbations  can  be  qualitatively  studied  using  the  ideas 
from  perturbation  of  group  action  as  well  as  the  singularity  theory  for  symmetric  systems.  While 
the  results  for  universal  unfolding  of  positive  definite  matrices  and  the  behavior  of  eigenvalues 
for  symmetry  preserving  perturbations  are  available  [15],  those  for  perturbations  that  destroy 
symmetry  are  not,  and  we  will  later  present  some  examples  displaying  the  interesting 
consequences  of  various  types  of  perturbations.  Finally,  neither  the  giuup  theory,  nor  the 
singularity  theory,  provide  quantitative  results  such  as  formulas  for  the  computation  of  the 
perturbed  eigenvalues  and  eigenvectors  as  a  function  of  the  perturbation  parameters.  Only  such 
information  can  provide  the  measures  for  eigenvalue  loci  veering  and  mode  localization,  and  one 
possible  tool  for  developing  these  expressions/results  is  the  singular  perturbation  theory.  Thus, 
tools  or  ideas  from  the  disciplines  of  group  theory,  singularity  theory,  and  singular  perturbation 
theory,  are  all  needed  to  make  a  strong  headway  in  understanding  the  phenomenon  of  mode 


localization. 
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2.  GROUP  THEORY  AND  CYCLIC  SYMMETRY 

Although  it  had  been  observed  that  turned  bladed-disk  assemblies  always  have  many  pairs  of 
degenerate  eigenvalues,  no  theoretical  justification  for  this  phenomenon  was  available  in  the 
bladed-disk  literature.  Our  first  order  of  business  in  the  investigation  was  therefore  to  obtain  a 
formal  explanation  for  this  phenomenon.  The  coefficient  matrices  in  the  equations  of  motion  of 
forced  bladed-disk  assemblies  as  well  as  the  dynamic  stiffness  matrices  are  always  banded 
circulant  matrices  [12].  These  matrices  have  unique  symmetry  properties  [13]  which 
immediately  indicate  that  group  theory  would  be  applicable.  It  turns  out  that  the  set  of  allowed 
symmetry  operations  in  a  bladed-disk  assembly  namely:  rotations  about  a  fixed  axis,  reflections 
about  a  fixed  axis  and  vibrations  about  a  reference  point,  can  be  captured  by  the  operations  of  the 
Dihedral  group  Dn  [14].  By  purely  formal  arguments  from  group  theory  and  standard  results  for 
the  Dihedral  group,  we  are  able  to  show  the  number  and  order  of  degenerate  eigenvalues  which 
any  finite  order  bladed-disk  assembly  can  have.  Furthermore  by  considering  the  irreducible 
representations  to  which  the  translational,  rotational,  and  vibrational  modes  belong,  along  with 
the  corresponding  Hamiltonians,  we  can  sufficiently  study  the  effects  of  mild  perturbations  on 
these  degenerate  doublets.  For  example  under  a  given  parameter  variation,  the  symmetry 
operations  generate  a  new  group  which  is  necessarily  a  subgroup  of  the  original  group  Dn.  By 
comparing  the  properties  of  this  new  sub-group  witn  those  of  the  original  group,  we  are  able  to 
determine  if  such  a  perturbation  would  lead  to  a  splitting  of  any  of  the  degenerate  pairs  of 
eigenvalues.  It  may  then  be  possible  to  determine  the  minimum  number  of  parameters  which 
must  be  varied  simultaneously  in  order  for  a  certain  number  of  degenerate  eigenvalue  pairs  to  be 
split  at  the  same  time.  By  now  concentrating  on  those  perturbations  that  lead  to  splitting  of 
degenerate  pairs  we  can  more  fully  study  the  effects  of  these  perturbations  on  the  forced 
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amplitude  response  of  the  assembly.  These  results  are  summarized  in  the  paper  in  Appendix  1. 

3.  SINGULARITY  THEORY  AND  CYCLIC  SYMMETRY 

It  is  known  from  singularity  theory  that  the  splitting  of  eigenvalues  of  a  matrix  can  lead  to 
rapid  changes  in  the  eigenvectors,  which  in  turn  can  result  in  significant  changes  in  the  forced 
amplitude  of  response  of  the  assembly  to  external  aerodynamic  loading.  Since  the  group  theory 
results  indicate  that  bladed-disk  assemblies  could  only  have  degenerate  pairs  of  eigenvalues,  the 
simplest  essential  properties  of  any  finite  order  bladed-disk  assembly  are  capture  i  '  v  the 
properties  of  an  assembly  of  order  3.  Note  that  we  need  a  3rd  order  assembly  in  order  to  inscribe 
a  circle  and  hence  obtain  a  cyclic  system.  A  third  order  tuned  assembly  would  thus  have  a 
degenerate  pair  of  eigenvalues  and  an  isolated  eigenvalue.  We  may  therefore  study  the  influence 
of  perturbations  in  the  masses,  ground  springs  and  coupling  springs  on  the  dynamics  of  this 
system.  From  Arnold’s  results  in  singularity  theory  [16],  it  is  self  evident  that  under  mild 
parameter  variations  interest  should  be  concentrated  not  on  the  isolated  eigenvalue  but  only  on 
the  subsequent  behavior  of  the  degenerate  doublet.  To  understand  its  behavior  it  is  necessary  to 
study  the  behavior  of  any  arbitrary  doublet  and  the  subsequent  eigenloci  as  a  function  of 
parameters  in  a  manner  reminiscent  of  root  loci  behavior  in  classical  control  theory’.  The 
simplest  doublet  which  contains  the  essential  ingredients  of  the  problem  turns  out  to  be  the 
symmetric,  coupled  double  pendulum  shown  in  figure  1  in  Appendix  2.  The  essense  of  this 
study  was  to  discover  the  relation  of  the  eigenloci  to  parameter  variations  and  the  corresponding 
eigenvectors.  We  had  conjectured  that: 

(i)  Uneven  amplitudes  of  vibration  in  symmetric  structural  systems  are  caused  by  the 
sensitive  dependence  of  system  eigenvectors  or  parameters. 
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(ii)  Sensitive  dependence  of  eigenvectors  on  parameters  gives  rise  under  appropriate 
conditions  to  the  mode  localization  phenomenon. 

(iii)  Rapid  convergence-divergence  (veering)  of  eigenvalues  is  a  signature  for  the  sensitive 
dependence  of  eigenvectors  on  parameters  and  hence  of  the  possible  existence  of  mode 
localization  under  appropriate  conditions. 

If  the  conjecture  were  to  be  true,  we  hoped  to  be  able  to  obtain  an  estimate  for  the 
eigenvector  sensitivity  measure  in  an  appropriate  manner  as  well  as  an  estimate  for  the 
eigenvector  rotations  resulting  from  any  mild  perturbation.  If  extreme  imperfection  sensitivity 
were  present,  it  was  expected  that  both  measures  would  show  a  singularity  which  is  an  indication 
of  imperfection  sensitivity.  Furthermore  these  were  expected  to  occur  at  the  parameter  values  of 
maximum  curvature  of  the  eigenloci.  If  the  double  pendulum  were  decoupled  (no  coupling 
spring)  there  then  would  exist  two  independent  but  equal  vibration  frequencies.  By  including 
very  weak  coupling  between  the  masses  we  could  study  system  behavior  in  the  neighborhood  of 
the  erstwhile  equal  eigenvalues.  The  study  of  weakly  coupled  systems  is  very  important  since  in 
practice  the  aim  has  always  been  towards  use  of  rigid  disks,  in  effect  making  the  inter-blade 
coupling  very  weak  indeed.  By  assuming  that  the  imperfection  parameter  is  a  slight  difference 
in  the  length  of  the  two  pendula,  we  could  then  study  the  behavior  of  the  eigenvalues  and 
eigenvectors  of  this  simple  symmetric  system  under  slight  changes  in  coupling  and  disorder. 

4.  QUANTITATIVE  UNFOLDING  OF  THE  MODAL  SINGULARITIES  BY  SINGULAR 
PERTURBATION  ANALYSIS 

We  may  therefore  write  down  the  characteristic  polynomial  as  a  function  of  both  the 
coupling  parameter  and  the  imperfection  parameter.  The  characteristic  polynomial  in  turn 


-  12- 


behaves  identically  to  a  potential  function  in  singularity  theory  [17].  Thus  the  degenerate 
critical  points  of  this  function  correspond  to  the  repeated  eigenvalues  if  any.  By  writing  down 
the  expressions  for  the  eigenvalues  as  functions  of  the  two  parameters,  the  detailed  behavior  of 
the  eigenloci  in  any  neighborhood  can  be  obtained.  A  regular  expansion  of  these  eigenvalues  as 
a  function  of  the  two  parameters  breaks  down  (loses  uniformity)  in  the  neighborhood  of  the 
critical  point  (where  the  eigenvalues  are  coincident).  By  applying  the  techniques  of  singular 
perturbation  analysis  and  appropriate  stretching  transformations  it  became  possible  to  obtain  the 
eigenloci  expressions  which  were  uniformly  valid  over  the  domain  of  definition  of  the  small 
parameters  and  whose  loci  clearly  indicated  the  veering  phenomenon.  The  same  technique  was 
also  applied  to  the  eigenvector  expressions  as  functions  of  the  parameters.  From  these 
expansions,  expressions  were  obtained  both  for  the  modal  sensitivity  measure  and  the 
eigenvector  rotation  measure  under  slight  parameter  variations.  All  the  results  obtained, 
confirmed  the  conjecture.  The  first  part  of  these  results  are  to  appear  in  the  Journal  of  Sound  and 
Vibration  while  the  second  part  involving  the  full  eigenvector  work  has  been  submitted  to  the 
Journal  of  Sound  and  Vibration.  These  manuscript  preprints  are  enclosed  in  Appendix  2.  We 
can  now  claim  that  we  understand  fairly  well  the  causes  of  localization  phenomenon  for  simple 
doublets.  We  note  however  that  the  double  pendulum  analysis  displayed  the  noted  strange 
behavior  only  under  very  weak  coupling  conditions.  Our  work  had  shown  earlier  (see  Appendix 
4)  that  for  bladed-disk  assemblies  uneven  amplitudes  of  vibration  and  hence  mode  localization 
could  occur  even  under  very  strong  coupling  conditions.  We  therefore  had  to  discover  under 
what  conditions  the  double  pendulum  analysis  remained  valid  also  for  the  bladed-disk  assembly. 
It  turns  out  that  the  pendulum  analysis  remains  valid  for  the  bladed-disk,  irrespective  of  the 
coupling  strength.  However  by  systematically  reducing  the  coupling,  more  complicated 
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singularities  appear.  This  is  because  at  very  low  coupling  the  bladed-disk  behaves  like  a 
perturbation  of  a  triple  degeneracy  which  in  general  requires  at  least  three  parameters  to  unfold 
(completely  analyse).  On  the  other  hand  for  the  double  pendulum,  two  parameters  were  enough 
to  unfold  the  doublet  degeneracy. 

The  lessons  learned  so  far  are  thus  that  for  the  bladed-disks  even  under  strong  coupling 
conditions  mode  localization  can  occur.  On  the  other  hand  for  the  double  pendulum  or  linear 
chains  in  general,  mode  localization  occurs  only  under  weak  coupling.  This  finding  contradicts 
the  current  view  in  bladed-disk  research  [11],  which  holds  that  in  both  linear  chains  (coupled 
pendula)  and  cyclic  chains  (bladed-disk  assemblies)  mode  localization  only  occurs  under  very 
weak  coupling.  What  is  however  true  for  bladed-disk  assemblies,  is  that  under  very  weak 
coupling  new  singularity  types  (which  do  not  exist  under  strong  coupling)  appear.  We  do  not  yet 
understand  the  full  effects  of  these  new  singularity  types.  We  however  conjecture  that  they  will 
further  complicate  the  modal  behavior  of  the  assembly  under  aerodynamic  loading.  The  key 
question  we  seek  to  answer  presently  is  which  of  either  symmetry  breaking  perturbations  or 
coupling  induced  perturbations  have  more  influence  on  the  modal  behavior  of  a  bladed-disk 
assembly.  Are  there  regimes  where  each  has  more  influence  than  the  other  and  if  so,  what  is  the 
transition  region?  If  we  could  answer  these  questions  then  we  could  specify  apriori  the 
acceptable  range  of  coupling  so  that  design  effort  could  be  concentrated  on  symmetry  breaking 
bifurcations  and  how  to  prevent  their  effects  from  being  felt  at  the  blade  amplitudes.  The 
singular  perturbation  analysis  acts  as  an  unfolding  of  the  singularities  involved,  since  by  this 
methodology  we  are  able  to  obtain  detailed  information  on  the  modal  behavior  of  the  structure  in 
the  neighborhood  of  the  singularities. 
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5.  THE  CONTROL  PROBLEM 

We  have  not  carried  out  the  control  design  component  of  the  project  as  stated  in  the 
statement  of  work  because  it  has  only  been  in  the  last  six  to  eight  months  that  a  thorough 
understanding  of  the  structural  dynamics  has  emerged.  We  however  are  clear  on  the  work  that 
needs  to  be  done.  The  next  stage  of  our  work  will  involve  classification  of  different 
perturbations  with  the  corresponding  amplitudes  of  vibration.  The  control  problem  in  one 
possible  approach  is  a  structural  redesign  that  deliberately  breaks  the  symmetry  by  splitting  the 
degenerate  eigenvalues  with  only  those  perturbations  that  do  not  lead  to  amplification  of 
vibration  amplitudes.  Provided  the  split  eigenvalues  are  not  in  the  neighborhood  of  the 
bifurcation  set,  all  further  slight  perturbations  would  not  be  expected  to  display  extreme 
imperfection  sensitivity.  Another  alternative  control  methodology  which  we  are  presently 
considering  is  a  regular  adapative  control  scheme  that  seeks  by  means  of  active  addition  or 
subtraction  of  control  masses  and  springs  to  restore  symmetry  whenever  the  symmetry  breaking 
signature  is  observed.  Under  this  scheme  the  degree  of  sensitivity  and  eigenvector  rotation  will 
determine  the  amount  of  modification  called  for  and  the  location  where  to  apply  it.  However 
this  kind  of  scheme  seems  to  us  to  be  more  appropriate  for  aerospace  structural  systems  than  to 
turbine  rotor  disks. 

6,  CONCLUSIONS 

The  results  obtained  from  the  study  in  the  last  two  years  have  helped  to  clarify  and  unify 
several  conflicting  viewpoints  within  the  bladed-disk  research  community.  What  is  more 
significant  is  that  it  has  led  to  a  better  understanding  of  the  potentially  very  complicated 
dynamical  structure  which  ensues  when  geometric  (spatial)  symmetry  interacts  with  weak 
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coupling  in  periodic  structures.  Without  this  understanding  any  attempts  at  either  structural 
redesign  or  structural  control  of  such  systems  would  inevitably  be  frought  with  danger.  We  are 
currently  continuing  work  on  the  forced  response  of  cyclic  systems  with  a  view  to  a  more 
complete  mathematical  characterization  of  the  relationship  between  amplitudes  of  vibration, 
mode  localization,  and  perturbation  type.  We  are  also  generalizing  the  singulai  perturbation 
approach  to  linear  chains  and  cyclic  symmetric  systems  of  any  finite  order.  We  believe  that  the 
development  of  the  structural  control  schemes  would  be  worthless  without  this  full 
understanding. 
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STATEMENT  OF  WORK 

The  principal  aim  of  the  proposed  research  is  to  carry  out  an  in-depth  mathematical  and 
numerical  investigation  of  the  dynamics  of  mistuned  cyclic  systems,  by  use  of  some  new  and 
extremely  powerful  topological  theory  of  dynamics,  and  to  develop  simple  control  schemes  for 
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preventing  unacceptable  vibration  characteristics  in  such  systems.  To  accomplish  this  task,  we 
will: 

(i)  Identify  the  topological  structure  of  nominally  tuned  bladed-disk  assemblies,  the  order  of 
the  degeneracy  in  the  natural  frequencies,  the  minimum  number  of  canonical  parameters 
needed  to  unfold  the  degeneracy,  and  the  classification  of  the  bifurcation  set  in  the 
parameter  space. 

(ii)  Use  the  Jordan-Amold  canonical  structure  theory  to  completely  characterize  all  the  blade 
motion  forms  expected  when  a  given  nominally  tuned  system  is  generically  mistuned. 

(iii)  Relate  the  canonical  unfolding  parameters  to  the  disk  assembly  elements  of  mass, 
generalized  damping  and  generalized  stiffness;  and  hence  determine  which  mistuning 
parameters  or  combinations  thereof,  govern  the  escalation  of  forced  response  amplitudes 
and/or  unacceptable  blade  motions. 

(iv)  Employ  the  control  methodologies  of  either  entire  eigenstructure  assignment  or  quadratic 
optimization  to  deliberately  mistune  the  blade  assembly  passively  so  that  eigenvalue 
degeneracy  under  slight  parameter  variations  are  avoided  and  at  the  same  time  the 
parameter  combinations  which  lead  the  assembly  to  unacceptable  blade  motions  (the 
bifurcation  set)  are  never  allowed  to  occur. 

(v)  Carry  out  a  thorough  numerical  simulation  on  typical  nominal  and  perturbed  bladed-disk 
assemblies  to  verify  and  validate  the  predictions  of  the  new  topological  theory. 

The  above  will  set  the  stage  for  a  controlled  laboratory  hardware  experimental  verification, 
which  we  hope  to  undertake  in  a  follow-up  project. 
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ABSTRACT 

In  this  work,  we  consider  the  dynamics  of  linear  mechanical  systems  possessing  geometri¬ 
cal  symmetry  subject  to  differential  or  small  parameter  variations.  The  machinery  of  group 
theory  including  the  irreducible  group  representations,  and  the  consideration  of  representations 
to  which  the  translational,  rotational  and  vibrational  modes  belong,  allow  us  to  predict  apriori, 
the  number  and  the  order  of  degenerate  eigenvalues  in  the  symmetric  system.  By  considering 
the  resultant  Hamiltonians  of  the  perturbed  symmetric  system,  we  show  further  the  effects  of  the 
perturbations  on  the  eigenvalues  and  their  degeneracies.  Since  the  vibration  modes  of  systems 
with  degenerate  eigenvalues  are  known  to  display  sensitive  dependence  on  parameters,  we  may 
use  these  techniques  to  identify  in  principle  the  possibility  of  maximum  vibration  amplitudes  and 
where  they  are  likely  to  occur.  Applications  of  these  ideas  include  the  mistuned  turbine  rotor 


bladed  disk  assemblies. 
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LIST  OF  SPECIAL  SYMBOLS 

T  ^  (R)  Reduced  cartesian  representation  of  a  group  element  R. 

T  t”d  (R)  Reduced  translational  representation  of  a  group  element  R. 

fr”d  Reduced  Rotational  representation. 

T  Vrad  Reduced  Vibrational  representation. 

HR),  r‘(R)  Matrix  representations. 
r‘(R)  ^representation. 

r^(R)  Matrix  element  of  the  IIth  row  and  the  Vth  column  of  the  matrix  representing  the 

group  element  R  in  the  i*  representation. 

r£v  (R)  Complex  conjugate  of  T^R). 

^*(R)  Character  of  a  group  element  R  in  the  i*  matrix  representation. 

a{  Number  of  times  T  '(R)  appears  in  the  reducible  representation, 

e  Row  vector. 

M  Mass  matrix. 

K  Stiffness  matrix. 
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1.  INTRODUCTION 

Eigenvalues  and  eigenvectors  of  a  vibrating  system  are  important  for  characterizing  its 
dynamical  response.  The  eigenvalues  are  related  to  natural  frequencies  whereas  the  eigenvectors 
correspond  to  special  forms  of  displacements  when  vibrating  at  a  natural  frequency.  Exact 
evaluation  of  the  eigenvalues  of  higher  order  vibrating  systems  in  general  involves  considerable 
effort  and  is  time  consuming.  Most  cyclic  symmetric  systems  possess  degenerate  eigenvalues 
[1,2].  Systems  with  degenerate  eigenvalues  are  expected  to  display  severe  sensitive  dependence 
on  parameters  [3]  that  destroy  the  symmetry  or  degeneracy. 

In  the  eigenvalue  problem  if  there  is  any  symmetry  of  the  system,  the  application  of  group 
theory  enables  us  to  decide,  at  the  outset,  exactly  the  number  of  distinct  eigenvalues  together 
with  their  respective  degrees  and  degeneracies. 

By  considering  the  symmetry  operations  of  the  physical  system  at  the  equilibrium  points, 
the  representing  group  can  be  formulated.  Using  group  theoretical  ideas,  we  can  predict  apriori 
the  degeneracy  of  the  eigenvalues.  This  is  accomplished  by  the  use  of  the  irreducible  represen¬ 
tations  of  this  group  which  is  obtained  by  using  the  orthogonality  theorem  and  the  reduction  for¬ 
mula  [3].  Once  the  irreducible  representations  are  known,  we  can  find  the  translational,  rota¬ 
tional  and  vibrational  modes  of  the  system.  These  results  are  well  known  in  the  physics  litera¬ 
ture  on  group  theory  but  have  not  been  used  sufficiently  effectively  in  the  vibration  community. 
The  essential  purpose  of  this  work  is  to  summarize  some  of  these  results  and  show  some  applica¬ 
tions  as  they  relate  to  the  symmetric  bladed  disk  assemblies. 

In  general,  a  reduced  cartesian  representation  of  a  group  element  R,  rc,-(R),  can  be  writ¬ 
ten  as 
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(R)  =  rt~*  (R)©rr-  (R)©rVi-  (R), 

where  rt~i(R),  F^fR),  and  rv"d(R)  arc  translational,  rotational  and  vibrational  reduced  carte¬ 
sian  representations  of  a  group  element  R.  In  three  dimensional  vibrating  systems,  six  of  the 
normal  modes  belong  to  the  zero  frequency  modes  and  correspond  to  pure  translations  and  pure 
rotations.  Since  we  are  primarily  interested  in  vibrational  modes  (non  zero  frequency)  the  zero 
frequency  modes  are  not  discussed  further  in  this  work. 

Some  of  the  symmetry  of  the  physical  system  may  be  lost  once  the  system  is  subjected  to 
parameter  perturbations.  The  perturbed  system  may,  however,  still  possess  some  symmetry 
which  may  be  considered  to  be  a  subgroup  of  the  group  characterizing  the  unperturbed  or  origi¬ 
nal  system.  Applying  group  theoretic  ideas  now  to  this  subgroup  we  can  predict  the  splitting  of 
the  degeneracy  of  eigenvalues.  As  a  result  we  can  see  whether  the  degeneracy  of  some  of  the 
eigenvalues  has  or  has  not  been  removed. 

Ideas  similar  to  the  ones  proposed  in  this  work  were  used  by  R.  Perrin  [4]  in  1971  for  a 
thin  circular  ring.  In  his  paper  group  theoretical  arguments  were  applied  to  a  ring  where  pertur¬ 
bation  was  applied  in  the  form  of  equal  masses  attached  to  the  ring  at  the  vertices  of  an  inscribed 
regular  n^  order  polygon.  Further,  eigenfrequencies  and  eigenfunctions  for  the  unperturbed  ring 
were  assumed  to  be  known  apriori.  Knowledge  of  these  degenerate  pairs  of  eigenfunctions  was 
used  to  find  the  characters  of  each  irreducible  representation  of  the  corresponding  Dn  group  for 
the  perturbed  system.  In  the  present  work,  group  theoretical  techniques  are  developed  without 
apriori  assuming  any  knowledge  of  either  the  eigenfrequencies  or  the  eigenfunctions  for  the  cir- 
culant  symmetric  system.  Also,  parametrically  perturbed  cases  were  not  discussed  by  Perrin. 
Parametric  perturbations  are  important  in  turbine  blade  vibration  problems  where  a  slight  pertur¬ 
bation  can  lead  to  loss  of  cyclic  symmetry,  which  in  turn  can  induce  rogue  blade  failure  under 
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certain  circumstances  [5,6]. 

2.  THE  GROUP  THEORETICAL  CONCEPTS 

We  first  define  some  standard  terms  from  the  literature  on  group  theory.  One  of  the  stan¬ 
dard  references  is  the  text  by  Hamermesh  [1].  Following  definitions  and  theorems  are  obtained 
from  [2]. 

2.1  Definition  1.  Symmetry  operations:  All  the  operations  which  leave  a  system 
configuration  unchanged  are  called  symmetry  operations. 

In  physical  terms,  this  refers  to  the  movement  of  a  system  in  such  a 
way  that  it  interchanges  the  positions  of  various  particles  of  the  system  but 
results  in  the  system  looking  exactly  the  same  as  before  the  symmetry 
operation.  For  instance  some  of  the  symmetry  operations  are  defined  as 
follows: 

E:  Identity.  The  system  is  not  rotated  at  all  or  rotated  by  2k  about  any 

axis. 

Cn  Rotation:  This  is  an  operation  which  effects  rotation  through  an 
angle  2rt/n  about  an  axis,  fixed  in  space,  where  n  is  an  integer.  In 
addition  we  can  have  C„,  which  is  Q  raised  to  the  power  k,  that  is, 
a  rotation  through  an  angle  2nkjn  about  the  same  axis.  C"  is  a  rota¬ 
tion  through  an  angle  2k  and  is  the  identity  operation,  since  a  rota¬ 
tion  through  2n  leaves  the  object  unchanged,  n  is  known  as  the 
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Definition  2. 


Definition  3. 


multiplicity  of  the  axis,  and  the  latter  is  called  on  n  fold  axis.  If 
n=2,3,...  then,  respectively,  we  get  2-fold,  3-fold...  axes.  If  a  sys¬ 
tem  has  more  than  one  axis  of  symmetry  then  the  axis  with  the 
highest  value  of  n  is  called  the  principal  axis. 

Group:  A  set  of  elements  {a,b,c..}  is  called  a  group  G,  if  a  multiplication 
rule  is  defined  for  any  two  elements  so  that  the  product  ab  has  a  definite 
meaning  and  the  following  four  postulates  are  satisfied: 

1.  Closure:  If  a  and  b  belong  to  the  set,  then  ab  also  belongs  to  the  set. 

2.  Associativity:  a 'be)  =  (ab)  c. 

3.  There  exists  the  identity  element  e  such  that  ae  =  ea  =  a  for  any  a 
belonging  to  G. 

4.  There  exists  the  inverse  element,  i.e.,  for  each  element  a,  there  is  a 
corresponding  element  b  such  that  ab  =  ba  =  e.  b  is  called  the 
inverse  element  of  a  and  is  denoted  by  b  =  a-1 . 

D„  group:  This  group  concerns  a  system  possessing  one  n-fold  axis  called 
the  principal  axis  and  n  2-fold  axes  symmetrically  placed  in  a  plane  per¬ 
pendicular  to  the  principal  axis.  The  n-fold  axis  provides  the  n  elements 
of  the  cyclic  group  C„.  The  group  also  contains  one  C2  element  provided 
by  every  perpendicular  2-fold  axis  where  we  do  not  count  C2  =  E  because 
E  only  occurs  once  in  a  set  of  group  elements  and  it  has  already  appeared 
in  the  n  elements  of  Q,.  The  group  Dn  therefore  contains  a  total  of  2n  ele¬ 


ments. 
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Definition  4.  Equivalent  and  reducible  representations:  Two  representations  are  said 

to  be  equivalent  if  the  two  matrices  representing  any  element  R  of  the 
group  are  related  by  the  equation 

r'  (R)  =  t-1  r  (R)  t,  (i) 


where  T  is  any  nonsingular  square  matrix  (operator).  However,  if  there 
does  not  exist  any  matrix  T  which  transforms  T '  (R)  into  T  (R),  then 
f '  (R)  and  T  (R)  are  said  to  be  inequivalent. 


A  reduced  representation  of  a  group  element  R,  Tred  (R),  is  com¬ 
posed  of  two  or  more  irreducible  representations: 


■red 


(R)  = 


rl(R) 

0 


0 

P  (R), 


We  write  this  by  using  the  symbol  ©  and 
rred  (R)  =  ri  (R)@rJ  (R). 


(2) 


(3) 


Definition  5.  X*(r)>  Character  of  a  group:  The  character  of  a  group  element  R  in  the  i* 

matrix  representation  of  the  group  element  is  the  trace  (sum  of  diagonal 
terms)  of  the  matrix. 


2.2.1  Orthogonality  Theorem:  All  the  vectors  formed  by  the  inequivalent  irreducible  unitary 
representations  are  orthogonal  to  each  other,  or: 

s  rj;v  (R)  riv  (R)  =  f 

R  *> 


where  i  and  j  denote  the  representation,  p.  and  p'  denote  rows  of  the  matrix  elements  and  v  and 
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v '  denote  the  columns  of  the  matrix  elements,  g  is  the  order  of  the  group  and  lj  is  the  dimen¬ 
sionality  of  the  i*  representation.  r^v  (R)  is  the  matrix  element  of  the  (Xth  row  and  Vth  column 
of  the  matrix  representing  the  group  element  R  in  the  i*  representation,  and  *  denote  its  complex 
conjugate  equivalent  to  T1  (R) 

2.2.2  Character  table:  The  character  table  is  formed  by  considering  the  characters  of  group 
elements.  The  character  of  a  group  element  is  important  because  the  character  is  unaltered  by  a 
similarity  transformation.  On  the  other  hand  since  the  character  of  equivalent  irreducible 
representations  are  identical  a  table  of  characters  is  a  unique  way  to  characterize  a  group.  The 
general  form  of  a  character  table  is: 


where 

c(r):  the  nature  of  elements  in  the  class, 

Nj :  number  of  elements  in  the  class,  and 
T 1 :  i111  irreducible  representation. 

The  character  tables  of  the  groups  can  be  obtained  by  the  application  of  the  following  rules: 
(1)  Number  of  inequivalent  representations  is  equal  to  the  number  of  classes. 
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7 

(2)  X  m  =  g>  n-  number  of  irreducible  representations. 

«=i 

(3)  X  X1*  OR)  Xj  (R)  =  g  8ij,  where  %‘(R)  denote  the  character  of  a  group  element  R  in  the  i* 

R 

matrix  representations. 

(4)  X  X**  (Ck)  X*  (Ck)  =  -rj~  5^. 

i=l  N* 


In  labeling  the  rows  of  the  character  table,  the  foUowing  standard  notation  is  used. 

(1)  One-dimensional  representations  are  labeled  as  A  if  the  character  of  the  elements  C„  about 
the  principal  rotation  axis  are  +1  for  all  k,  and  as  B  if  the  characters  C„  are  (-l)k  for  all  k. 

(2)  If  a  group  has  more  than  one  A  or  B  representation  they  are  given  subscripts  1  and  2 
according  to  whether  the  character  is  +1  or  -1  in  the  column  representing  a  rotation  or 
improper  rotation  about  an  axis  other  than  the  Principal  axis.  For  example,  in  the  groups 
Dn  a  representation  is  given  a  subscript  1  if  the  character  under  C2  about  the  axis  is  +1  and 
2  if  it  is  -1. 

(3)  Two  dimensional  representations  are  labeled  E. 

(4)  Three  dimensional  representations  are  usually  labeled  T. 


Finally,  the  reducible  representations  are  used  to  obtain  the  irreducible  representations  by  the 
applications  of  the  formula 


(4) 
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3.  VIBRATIONS  OF  CYCLIC  MECHANICAL  SYSTEMS 


Many  vibrating  systems  possess  sufficient  symmetry  to  allow  us  to  use  group  theory  which 
reduces  the  amount  of  work  involved  in  the  calculations  and  also  furnishes  us  with  an  insight 
into  the  nature  of  the  vibrations.  Consider  n  nominally  identical  masses  connected  via  ground 
springs  kt  and  coupling  springs  at  the  edges  of  an  inscribed  n*  order  regular  polygon  in  a  cir¬ 
cle  of  radius  r.  Rotational  symmetry  about  one  n-fold  axis  perpendicular  to  the  plane  of  motion 
and  n,  2-fold  axes  give  2n  number  of  elements  for  the  corresponding  symmetry  group.  In  fact, 

this  is  the  dihedral  group  Dn.  Dn  has  (2  +  —  *  )  conjugate  classes  when  n  is  odd  and  (3  + 

classes  when  n  is  even.  Utilizing  group  character  table  construction  rules  [1,2],  it  can  be  shown 
that  the  number  of  possible  degenerate  eigen  levels  are: 


n  even 


single  degenerate  levels  =  4 
^  double  degenerate  levels  =y  -  1 


n  odd 


single  degenerate  levels  =  2 
11')  double  degenerate  levels  =  ■^r— 


(5) 


(6) 


Results  (i)  and  (ii)  imply  that  a  cyclic  symmetric  system  of  this  type,  at  worst,  can  have  double 
degenerate  vibrational  modes  for  any  finite  n. 

Once  the  system  is  subjected  to  a  random  parameter  perturbation,  some  of  the  symmetry 
may  be  lost  and  consequently,  we  get  a  new  group  which  is  likely  to  be  a  subgroup  of  the  origi- 
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nal  group.  Equations  (5)  and  (6)  can  be  used  to  get  a  qualitative  information  about  the  new 
degenerate  eigenvalues  of  the  perturbed  system.  Since  any  higher  order  cyclic  symmetric  sys¬ 
tem  of  this  type  will  have  at  most  doubly  degenerate  eigenvalues  it  is  sufficient  to  consider  an 
example  with  n  =  3. 

Consider  three  normally  identical  masses  connected  via  both,  ground  springs  kt  and  cou¬ 
pling  springs  kc,  at  the  edges  of  an  inscribed  isosceles  triangle  in  a  circle  of  radius  r.  We  wish 
to  study  the  number  and  degeneracy  of  the  eigenvalues  of  such  systems  under  random  differen¬ 
tial  perturbations  in  the  elements  {kt,  k^  m}. 

As  defined  in  Section  2,  a  symmetry  operation  is  one  which  leaves  the  undistorted  system 
indistinguishable  from  its  previous  orientation.  Such  an  operation  interchanges  equivalent 
masses.  However,  in  the  vibrational  state,  the  system  is  in  a  distorted  configuration  and,  when 
the  symmetry  operation  is  performed  on  the  distorted  mass  the  effect  is  the  same  as  that  obtained 
by  interchanging  displacement  vectors  amongst  equivalent  masses.  Therefore  we  can  define  the 
action  of  a  symmetry  operation  for  each  mass  in  a  distorted  system  to  be  a  displacement  through 
vector  X,  from  its  equilibrium.  When  a  symmetry  operation  is  applied  we  can  assume  that  the 
mass  positions  remain  invariant.  Furthermore  the  symmetry  operation  can  have  no  effect  on  the 
potential  or  kinetic  energy  of  the  system,  or  even  the  angles  between  the  connections.  Conse¬ 
quently  the  quadratic  forms  of  the  kinetic  energy  T  and  the  potential  energy  V  remain  invariant 
under  the  action  of  the  group  transformations.  Group  theory  can  thus  be  used  to  determine  and 
classify  the  normal  modes  of  the  vibrating  system. 

We  begin  with  the  3N  dimensional  representation  of  the  group  of  symmetry  operations  of 
the  undistorted  system.  By  reducing  this  representation  using  the  character  table  for  the 
corresponding  group,  (and  reduction  of  reducible  representations)  we  can  determine  the 
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irreducible  representations  to  which  the  3N  translational,  rotational  and  vibrational  modes 
belong.  Also  we  can  immediately  find  the  degeneracy  of  each  normal  mode.  In  addition,  by 
considering  the  symmetrized  basis  S,  we  can  bring  the  mass  matrix  M  and  the  stiffness  matrix  K 
into  block  form  and  thus  greatly  simplifying  the  solution  of  the  characteristic  or  the  frequency 
equation. 

We  apply  these  group  theoretic  techniques  to  the  problem  shown  in  Figure  1  [4].  This  sys¬ 
tem  belongs  to  the  group  D3  which  contains  the  symmetry  operations  E,  C$,  C^,  C“,  C^,  and 
C2.  Now  by  applying  these  operations  to  the  nine  cartesian  coordinates  xj,  X2,  X3,  y3,  y2,  y3,  z\, 
z2>  z3  we  obtain  their  nine  dimensional  reducible  representations.  This  can  be  accomplished  by 
finding  the  corresponding  matrix  representation,  and  using  the  equation 

X'  =  r(R)X,  (7) 

where  R  is  the  appropriate  symmetry  operation,  and  X  and  X'  denote  the  nine  dimensional  vec¬ 
tors  rerresenting  the  cartesian  coordinates  in  original  and  transformed  planes  respectively.  For 
example,  under  the  operation  of  the  system  configuration  in  Figure  1  is  transformed  to  the 
configuration  in  Figure  2.  This  also  clearly  indicates  the  manner  in  which  the  coordinates 
undergo  rotation.  In  Figure  2,  the  axes  zlt  z2  and  z3  are  pointing  out  of  the  plane  of  the  paper. 
The  new  and  the  old  coordinates  are  related  by  the  relations 

(*t»  yt.zt)  (x’l.  Y\,  z'i)  ■(-  x2  sin  30  -  y2  cos  30,  x2  cos  30 -y2  sin  30,  z2),  (8) 

(x2,  y2,  z2)  -4  (x'2,  /2  z'2)  a  (-  x3  sin  30  -  y3  cos  30,  x3  cos  30  -  y2  sin  30,  z3) ,  (9) 

(x3.  y3.  z3)  —*  (x'3.  /3>  z'l)  =(  -  x3  sin  30  -  y3  cos  30,  x3  cos  30-  y2  sin  30,  z3) .  (10) 

In  matrix  form,  equations  (8)  -  (10)  are  represented  as 


0 

-1/2 

0 

0 

-V3/2 

0 

0 

0 

0 

0 

0 

-1/2 

0 

0 

-V3/2 

0 

0 

0 

-1/2 

0 

0 

-V3/2 

0 

0 

0 

0 

0 

0 

V3/2 

0 

0 

-1/2 

0 

0 

0 

0 

0 

0 

V3/2 

0 

0 

-1/2 

0 

0 

0 

-V3/2 

0 

0 

-1/2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

where  the  coefficient  matrix  is  the  matrix  representation  of  the  group  element  C3.  The  character 
(trace)  of  the  matrix  is  then 

Xc(C])  =  0.  (12) 

Since  Cf  and  C]  belong  to  the  same  %c  (C^)  =  0.  It  is  clear  that  (E)  =  9.  In  a  similar 
manner 

Xc(C5)  =  -i, 

xc  (C5)  =  -  1  ,  (13) 

xc  (C5)  =  - 1 . 

Now  using  the  character  table  for  D3  (Table  1)  and  equation  (4),  we  can  determine  the 
irreducible  components  of  Tc  (R)  as: 


Table  1 
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Here,  Ai,  A2  represent  one  dimensional  representations,  E  represents  two  dimensional 
representations,  and  Tx,  Ty,  Tz  and  Rx,  Ry,  Rz  represent  unit  translational  and  rotational  vectors 
respectively.  Consequently 

aAl  =  j  (9  x  1  +  2  (0  x  1)  +  3  (1  x  -1))  =  1 , 

aAz  =  j  (9  x  1  +  2  (0  x  1)  +  3  (-1  x  -1))  =  2  , 

aAj  =  1  (9  x  2  +  2  (0  x  -1)  +  3  (1  x  0))  =  3  . 

Hence  T0"4*  (R)  =  A\  ©  2A2  ©  3E.  Since  the  system  has  three  masses,  N  =  3  and  there  are  9 
degrees-of-freedom  for  the  system.  Since  3N  -  6  =  3,  we  have  only  three  vibrational  modes,  the 
other  six  corresponding  to  zero  frequency  modes  and  to  pure  rotations  and  pure  translations. 
This  can  be  seen  from  the  character  table, 

rt(R)  =  A2©E, 
r  (R)  =  A2©E. 

Since  we  are  primarily  interested  in  the  vibrational  modes  of  the  system  and 

T c"d  (R)  =  Ai  ©  2A2  ©  3E  , 
r*  (R)  ©  r  (R)  =  2A2  ©  2E, 

we  get  rv"d  (R)  =  Ai  @  E. 

Since  Aj  is  one  dimensional,  the  vibrating  system  has  one  non  degenerate  eigenvalue  and 
since  E  is  2  dimensional  there  is  a  one  degenerate  eigenvalue.  These  results  are  consistent  with 
the  exact  eigenvalues  given  in  appendix  A. 

Using  group  theory  arguments  to  predict  the  number  of  degenerate  eigenvalues  becomes 
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very  useful  as  the  order  of  the  system  increases  since  the  exact  calculation  of  eigenvalues  for 
such  high  order  systems  becomes  increasingly  burdensome. 

4.  THE  PERTURBED  SYSTEM 

Suppose  that  the  Hamiltonian  of  the  unperturbed  system  is  Hq.  Then  Ho  is  invariant  under 
its  symmetry  group  G.  Suppose  further  that  the  system  is  subjected  to  a  perturbation  with  Ham¬ 
iltonian  V.  The  perturbed  Hamiltonian  H  =  Hq  +  V,  will  then  have  a  symmetry  group  which  is 
necessarily  a  subgroup  of  G.  Two  possible  cases  arise. 

CASE  I 

If  the  perturbation  V  has  symmetry  at  least  as  great  as  Ho,  the  group  G  will  still  be  the 
symmetry  group  of  the  total  Hamiltonian  H.  In  this  case  the  possible  types  of  eigenvalues  will 
be  unchanged  by  the  perturbation.  In  fact  no  splitting  of  degenerate  levels  occurs. 

CASE  II 

If  the  perturbation  V  has  symmetry  lower  than  Ho,  the  total  Hamiltonian  H  will  have  a 
symmetry  group  G1  which  is  a  subgroup  of  G.  This  subgroup  G1  is  invariant  under  the  perturba¬ 
tion.  Because  of  the  perturbation,  some  of  the  degenerate  eigenvalues  may  split.  This  can  be 
explained  by  using  group  representation  theory. 

For  a  given  representation  D(G)  of  the  group  G,  we  now  obtain  the  invariant  subgroup  G1 . 
Even  if  D(G)  is  an  irreducible  representation  of  G,  the  representations  of  G1  which  we  derive  in 
this  way  may  be  reducible.  In  other  words,  even  though  we  cannot  find  a  subset  of  the  basis  vec¬ 
tors  of  D(G)  which  is  invariant  under  all  transformations  of  the  group  G,  we  may  be  able  to  find 


-  17- 


a  subspace  which  is  invariant  under  all  transformations  belonging  to  the  eigenvalue  X  from  a 
basis  of  an  irreducible  representation  of  G.  This  representation  may  be  reducible  for  the  sub¬ 
group  G1 .  The  perturbation  V  will  then  split  the  level. 

We  now  apply  the  above  ideas  and  show  the  appropriate  methodology  in  the  context  of  the 
three  mass  system. 

4.1  Ground  Spring  Perturbed  System 

First  we  consider  the  situation  when  one  of  the  ground  springs  is  perturbed.  The  unper¬ 
turbed  system  can  be  represented  by  the  group  D3  =  {E,  C3,  C^,  Ca,  Q,»  Cc }.  Once  the  system 
is  subjected  to  a  ground  spring  perturbation,  for  this  particular  system  G1  =  {E,  Ca }  is  the  invari¬ 
ant  subgroup. 

The  character  table  for  G1  =  {E,  Ca },  is  as  follows: 


E 

Ca 

A' 

1 

1 

A" 

1 

-1  . 

Considering  the  part  of  the  character  table  of  D3  which  refers  to  the  operations  of  the  sub¬ 
group  G1  =  (E,  Ca},  we  have 


E 


E 


2 


Ca 


0. 


Utilizing  (4),  the  irreducible  components  are  then  given  by 
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aA*  =  y(lx2xl  +  lxOxl)=l, 
aA"  =  1  (1  x  2  x  1  +  1  x  0 -1)  =  1  . 

Thus,  the  doubly  degenerate  E  level  of  the  unperturbed  system  splits  into  single  levels  A'  and  A" 
of  G1  under  the  ground  spring  perturbation.  As  a  result,  degenerate  eigenvalue  of  the  perturbed 
system  separates.  Hence,  for  this  particular  system  we  get  three  distinct  eigenvalues.  Coupling 
spring  perturbation  leads  to  a  case  where  there  is  no  invariant  subgroup  left  and  consequently, 
group  theoretical  arguments  do  not  work  for  this  particular  situation.  We  conjecture  that  this  is 
indicative  of  those  cases  where  perturbations  do  not  lead  to  radical  changes  in  the  eigenvector 
directions. 

4.2  The  Mass  Perturbed  System 

A  system  consisting  of  three  particles,  two  with  mass  m  and  the  other  with  mass  M,  is 
illustrated  in  Figure  3.  Considering  rotational  symmetry  of  the  system,  we  can  see  from  Figure  3 
that  this  system  belongs  to  the  group  D2 

D2  =  {E,  C2}. 

The  character  of  each  of  the  elements  of  the  group  can  be  obtained  from  the  reducible  represen¬ 
tation  whose  matrix  representation  is  obtained  by  the  use  of  the  coordinate  transformation 

X1  =  r  (R)  X. 

Performing  these  operations,  and  following  the  steps  along  the  lines  of  work  in  section  3,  we  can 
show  that  the  reducible  characters  of  each  element  are 


Xc  (C2)  =  -  1  , 
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XC(E)  =  9. 

Using  the  character  table  of  D2  and  (4),  the  irreducible  components  of  the  reducible  representa¬ 
tions  can  be  then  determined  as 


d2 

E 

c2 

A 

1 

1 

z 

Rz 

B 

1 

-1 

x,y 

RX,  Ry 

x  1  -  1  x  1 


x  1  +-1  x  - 


Therefore, 

rc™“  (R)  =  4A@5B  . 

In  general,  T  **  (R)  =  T  t"d  (R)  ©  (R)  ©  Tv~d  (R). 


(14) 


By  placing  a  coordinate  system  XYZ  at  the  center  of  mass,  translations  and  rotations  can  be 
represented  as  shown  in  Figure  4.  The  representations  for  C2  and  E  are  given  by 

Xl(C2)  =  ~l  .  X‘(E)  =  3, 

Xr(C2)=-l  ,  Xr(E)  =  3. 

By  the  application  of  the  a,  equations  (4),  we  get 


(i)  Translation 
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ak  =  y[lx3-lxl]=l, 
a^  =  y[lx3-lx-lj  =  2. 

(ii)  Rotation 

ak  =  y  [lx3-lxl]=l  , 
arB  =  y  [1  x3  -  lx— 1]  =  2. 

Therefore, 

(R)  =  A@2B  ,  (15) 

rf"d(R)  =  A@2B.  (16) 

Hence  from  equations  (14)-(16),  we  get  the  result  that 

rVrcd  (R)  =  2A@B.  (17) 

This  shows  that  there  are  three  vibrational  modes  and  each  eigen  level  is  non-degenerate  since  A 
and  B  are  one  dimensional  representations. 

5.  SUMMARY  AND  CONCLUSIONS 

This  work  uses  results  from  group  theory  and  applies  it  to  perturbed  cyclic  symmetric 
vibratory  systems.  It  is  shown  that: 

a.  The  number  and  order  of  degenerate  eigenvalues  in  a  symmetric  system  can  be  predicted 
apriori  by  using  group  theory  without  explicitly  determining  the  eigenvalues. 

b.  Cyclic  symmetric  vibrating  systems  possess  degenerate  eigenvalues  for  n  >  2.  For  strong 
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coupling,  these  eigenvalues  occur  in  doubly  degenerate  pairs  and  in  single  nondegenerate 
levels. 

c.  Random  parameter  perturbations  may  partially  or  totally  destroy  the  symmetry  of  the  sys¬ 
tem.  Accordingly  these  perturbations  lift  some  of  the  degeneracy  of  eigenvalues.  As  a 
result,  eigenvalue  loci  veering  [7]  occurs  when  the  parameters  are  continuously  varied.  This 
may  also  lead  to  a  mode  localization  or  rapid  variation  in  the  eigenfunctions  [7]. 
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Appendix  A 

Equations  of  motion  for  the  system  in  Figure  2  can  be  written  in  the  form: 

Mx  +  Kx  =  0  ,  where 


m  0  0 

M  = 

0  m  0 

,  K= 

-1^  a 

0  0m 

kc  kc  a 

and 


kt 

3  =  2^  +  -^- 


The  eigenvalues  are  determined  by 

det  |K-C02  M|  =0. 

Therefore, 


3  kc 
m 


3  kc 
m 


and 


kt 

mr2 


This  shows  that  the  cyclic  symmetric  system  has  a  double  degenerate  eigenvalue  and  a  sin¬ 
gle  nondegenerate  eigenvalue. 
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Figure  1.  Model  representing  the  unperturbed  three  bladed  disk  assembly  in  its  identity 
orientation. 

Figure  2.  Model  of  figure  1  counter/clockwise  rotated  by  2rt/3  radians  about  the  center  of  the 
system,  showing  change  in  cartesian  coordinate  axes. 

Figure  3.  Mass  perturbed  system  in  its  identity  orientation. 
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coordinate  system. 


Figure  2.  Model  of  figure  1  counter/clockwise  rotated  by  2it/3  radians  about  the  center  of  the 
system,  showing  change  in  cartesian  coordinate  axes. 
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In  recent  years  there  has  been  tremendous  interest  in  the  vibrations  and  structural  dynamics 
community  in  the  phenomenon  of  mode  localization.  This  interest  stems  from  the  recognition 
that  large  systems  composed  of  nominally  identical  subsystems  inevitably  involve  minor  devia¬ 
tions  from  the  idealized  structures  and  these  disorders  or  perturbations  can,  under  appropriate 
conditions,  cause  disproportionately  large  deviations  from  the  predicted  behavior  in  the  nominal 
or  idealized  system  modes.  Important  technical  applications  of  these  include  mistuned  bladed 
disk  assemblies  [1]  and  large  space  structures  [2|. 

It  is  well  understood  by  now  that  the  presence  of  small  irregularities  in  nearly  periodic  struc¬ 
tures  may  inhibit  the  propagation  of  vibration  and  localize  the  vibration  modes.  Depending  on 
the  magnitude  of  perturbations  (disorder)  and  on  the  strength  of  internal  coupling  between  the 
subsystems,  the  mode  shapes  may  undergo  dramatic  changes  to  become  strong’  localized  when 
small  perturbations  are  introduced,  thereby  confining  the  energy  associated  with  a  given  mode  to 
a  small  geometric  region.  This  phenomenon  is  referred  to  as  mode  localization.  Pierre  [3] 
showed  that  strong  mode  localization  and  eigenvalue  curve  veering,  are  two  manifestations  of 
the  same  phenomenon.  Therefore,  the  investigation  of  the  curves  of  the  eigenvalues  or  natural 
frequencies  in  the  neighborhood  of  the  ordered  state  is  sufficient  for  determining  the  occurrence 
of  strong  mode  localization.  That  eigenioci  veering  phenomenon  can  occur  in  disordered  struc¬ 
tures  under  certain  conditions  has  been  explained  qualitatively  using  geometric  arguments  in  [4]; 
where  it  was  also  hinted  that  quantitative  results  should  be  obtainable  by  the  use  of  singular  per¬ 
turbation  analysis. 

For  mistuned  linear  chains,  Pierre  [3]  and.  Pierre  and  Dowell  [5]  showed  that  the  straight  for¬ 
ward  expansion  in  terms  of  mistuning  parameters  breakdown  in  the  case  of  weak  coupling.  This 
arises  because  the  idealized  system  that  is  being  perturbed  has  natural  frequencies  with 
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multiplicity  >  1.  They  then  developed  a  so-called  "modified"  perturbation  technique  which  pro¬ 
vided  a  gccd  approximation  to  the  exact  eigenfrequencies  and  showed  good  agreement  with 
experimental  results. 

In  the  present  note  we  show  that  the  singularity  causing  the  breakdown  of  the  straight  for¬ 
ward  expansion  can  be  analyzed  by  the  '•veil  developed  singular  perturbation  techniques  [6]  and 
an  appropriate  asymptotic  expansion  for  ;he  eigenfrequencies  can  be  constructed  which  provides 
a  correct  qualitative  and  good  quantitative  approximation.  In  order  to  explain  the  ideas  and  iu 
keep  the  algebraic  manipulations  to  a  minimum,  the  attention  is  focused  on  the  now  standard 
example  [3]  of  the  coupled  penduli  shown  in  Figure  1. 

The  basic  idea  of  the  technique  is  the  following:  by  applying  the  regular  perturbation  tech¬ 
nique  to  the  characteristic  equation  F(X,  e,  8)  =  0,  of  the  system,  we  can  obtain  algebraic  expres¬ 
sions  for  the  natural  frequencies  (eigenvalues)  as  a  power  series  in  the  small  parameter  or  pertur¬ 
bation  (say  5).  The  coefficients  of  the  power  series  are  dependent  on  the  second  parameter  e  and 
these  expansions  are  valid  for  all  values  of  e  so  long  as  no  singularities  arise.  Singularities  occur 
for  values  of  E  where  the  eigenfrequencies  lose  their  smoothness  and  it  is  said  that  the  expansion 
is  not  uniformly  valid  for  all  E.  Away  from  the  singular  parameter  (e)  values,  the  straightforward 
expansions  are  good  approximations  and  are  called  the  "outer  expansions”.  The  neighborhood 
of  the  singular  parameter  point  is  then  stretched  or  rescaled  in  terms  of  a  new  parameter  so  as  to 
remove  the  singularity.  The  expansions  in  terms  of  the  new  parameter  is  valid  only  in  the  neigh¬ 
borhood  of  the  singular  point  and  is  called  the  "inner  expansion".  The  inner  and  the  outer  solu¬ 
tions  can  be  matched  where  their  domains  of  validity  overlap  and  then  a  composite  expansion 
can  be  constructed  which  is  valid  uniformly  throughout  the  function  domain  for  all  values  of  the 


parameter  e. 
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Consider  the  system  of  two  weakly  coupled  penduli  system  as  shown  in  Figure  1 .  The  two 
important  parameters  are  the  dimensionless  coupling  between  pendulums  R2  =  (k/m )/( g/1),  and 
the  dimensionless  length  change  Al.  The  corresponding  eigenvalue  problem  generated  by  the 
above  system  is  given  by: 


1+R2  -R2 
-R2  R:+(l  +A/)_1 


O  = 


where  R2  =  co^/g)2  ,  G>£  =  k/m  .  co'  =  g //  . 


( 1 . 1 1 


This  eigenvalue  problem  results  in  the  following  characteristic  equation: 


F0.,e,5)  =  X2  -  (1+2&*-— )X.  +  —  +  6  =  0, 

1 +£  1  +£ 


(1.2) 


where  A /  =  £  and  R2  =  5. 


We  can  express  the  solutions  to  (1.2)  as  regular  functions  of  the  parameters  A/,  6  as  follows: 


A.[  (A/)  —  1  +  6  + 


1  + 


A/ 


52  +  0(53)  sA.f  (A/) 


(1.3) 


52  +  0(53)  =  Xj  (A/) .  (1.4) 

The  expressions  (1.3),  (1.4)  are  the  regular  expansions  of  the  eigenvalue  problem  for  small  cou¬ 
pling.  When  A /  — >  0,  X\  and  X2  become  unbounded  and  the  continuity  of  the  eigenvalues  with 
respect  to  the  perturbation  A /  breaks  down,  as  shown  in  Figure  2.  Each  eigencurve  has  two 
branches,  one  valid  for  Al  >  0  and  the  other  for  A/  <  0.  These  branches  are  indicated  by  the 
superscripts  ’+’  and  which  correspond  to  Al  >  0  and  Al  <  0.  respectively.  Note  that,  since  Al 
and  R2  =  5  have  been  treated  as  two  independent  parameters,  we  have  no  control  over  expres- 


X2(A  /)  = 


1 


l  +  A 1 


+  5- 


1  + 


Al 
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sior.s  (1.3)  and  (1.4)  In  the  limiting  process  when  11  -»  0  and  5  — >  0.  By  forming  the  inner 
expansion"  however,  we  can  find  an  exact  relation  between  these  two  parameters  by  taking  into 
consideration  the  nature  of  the  singularity.  Then  11  and  o  become  dependent  parameters 
Asymptotically  matching  the  inner  and  the  outer  expansions,  then  gives  the  composite  expan¬ 
sions  which  are  valid  throughout  the  region  of  interest. 

For  the  in.n.er  ^pmsion.  u.e  assume  that  the  physical  parameters  R~  =  6  and  V  -  <-: 
related  (dependent)  by  a  set  of  mathematical  parameters:  Ci,  .  .  and  u  by  a  "stretcning 

transformation  of  the  form: 

50 

e  =  e0 +q|ia  +  £ ,  d.5i 

j-- 

where  pi  is  a  new  small  parameter  that  is  defined  by: 

5(|i)  =  (sgn5)|ib  .  (1.6) 

The  positive  constants  a  and  b  are  to  be  determined  by  the  nature  of  the  singularities  of  FT.. 

£,  0)  near  e  =  £o,  where  £q  IS  the  singular  point  of  interest.  Let  the  dependent  variable 

ztjrj  =  A.(e(|i),5(|i))  be  written  as  the  expansion 

z=£zjfiJ.  IT  I 

J=0 

Note  that  for  the  pendulum  problem  Ft  a.  £.  0)  =  0  has  a  singular  point  at  £  =  £,,  =  <>.  The 
expansions  (1.7)  are  called  the  inner  expansions"  and  the  ZjS  are  called  the  inner  coefficients 
Substituting  ( 1.5).  ( 1.6)  and  ( 1.7)  into  ( 1.2).  and  by  simplifying  the  inner  expansions  aith  a=b=2 
one  obtains  the  following  solutions 
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2-1;  +  V^2  +  4 

tt2  ,  q2 

r  ^ 

i-  * 

z.[  ’ll 

L  2  J 

4  +  2 

.  Vq2+4. 

Z2  —  1  *+■ 


4. 

r  x 

2-q  -  V^2  +4 

!l2 

Li2  + 

1  +  q 

2 

„  Vq2  +  4  _ 

(i4  +0(^i5) , 
jj.4  +  0(fi5) , 


A  f  ** 

M  =  e  =  cir  , 


(1.9) 


(1.10) 


and 

R2  =  5  =  (i2  .  d.ll  i 

Keeping  (i  fixed  and  taking  the  limit  |q|  — >  <».  it  is  easy  to  see  that  z\  matches  asymptoti¬ 
cally  with  a. ^  for  q  — ■»  oo  and  with  for  q  — »-oo.  Similarly  the  inner  solution  matches 
asymptotically  ( |q  |  with  X[  and  X|.  Now  combining  the  inner  and  the  outer  expansions 

appropriately,  we  get  the  composite  expansions: 


X,  = 


1  +R2  +(l  +  —  )R4 
Al 


(1  -u(A/))  + 


1  + 


2R2  -  Al  -  R2V(A//R2)2  +4 


Al2  1 


+ 


Al 


R 


W(A//R2)2 


+  4 


— L7  +  R2-(1  +  -L)r4 

1  +  Al  Al 


u(A/)-(l  -u(A/)) 


1  +  R2  +  (1  +  -i-)R4 


-  u(A /) 


l  -  Al  +  Al2  +  R2  -  (1  +  — )R4 

Al 


+  0(R6) 


(1.12) 
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— 


— 5—  +  R2  -(1  +  —  ;R4 
1+A/  A/ 


( 1  -  u(A/)) 


j  2R2-A/+R2V(A//R2)2  +4 


+ 


+ 


A/2  h_ 


A/ 


1 


R2V(A//R2)2  +4  J 


l  +  R“  +  ( 1  +  — iR 
A/ 


4|uiA/>-  !:  l-A/+Ai2+R2-ll  +  —  )R4 


A/ 


(1  -  u( A/)) 


-  [1  +  R2  +  (1  +  -y  iR4  |u(A/)  •+■  Q(  R6  )  . 


where 


u(A/)  = 


1  ,  A/  >  0 


0  ,  A/  <  0 

The  plots  of  eigenfrequencies  ,  a_2  versus  A /  are  given  in  Figure  3  for  both  the  exact  solu¬ 
tions  and  the  solutions  obtained  above  by  the  singular  perturbation  technique.  These  are  in 
excellent  agreement.  Thus,  the  singular  perturbation  technique  leads  to  qualitatively  correct 
asymptotic  approximations  that  are  often  close  to  true  solutions  and  can  be  used  as  a  mathemati¬ 
cal  tool  to  generate  quantitatively  accurate  solutions  for  a  wide  variety  of  linear  and  nonlinear 
structural  dynamics  problems.  The  methodology  is  general  and  systematic  and  when  combined 
with  elementary  singularity  theory,  should  provide  a  powerful  technique  to  study  the  mode  local¬ 
ization  phenomenon  in  any  finite  order  linear  or  cyclic  dynamic  chain. 
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ABSTRACT 

An  investigation  of  the  eigenvalue  loci  veering  and  mode  localization  phenomenon  is 
presented  for  mistuned  structural  systems.  Examples  from  both,  the  weakly  coupled  uniaxial 
component  systems  and  the  cyclic  symmetric  systems,  are  considered.  The  analysis  is  based  on 
the  singular  perturbation  techniques.  It  is  shown  that  uniform  asymptotic  expansions  for  the 
eigenvalues  and  eigenvectors  can  be  constructed  in  terms  of  the  mistuning  parameters  and  these 
solutions  are  in  excellent  agreement  with  the  exact  solutions.  The  asymptotic  expansions  are 
then  used  to  clearly  show  how  singular  behavior  in  the  eigenfunctions  or  modeshapes  leads  to 


mode  localization. 
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1.  INTRODUCTION 

In  recent  years  there  has  been  tremendous  interest  in  the  vibrations  and  structural  dynamics 
community  in  the  phenomenon  of  mode  localization.  This  interest  stems  from  the  recognition 
that  large  systems  composed  of  nominally  identical  subsystems  inevitably  involve  minor  devia¬ 
tions  from  the  idealized  structures  and  these  disorders  or  perturbations  can,  under  appropriate 
conditions,  cause  unexpectedly  large  deviations  from  the  predicted  behavior  in  the  nominal  or 
idealized  system  modes.  Important  technical  applications  where  these  problems  arise  include 
mistuned  bladed  disk  assemblies  [1,2]  and  large  space  structures  [3,4]. 

It  is  well  understood  by  now  that  the  presence  of  small  irregularities  in  nearly  periodic  struc¬ 
tures  may  inhibit  the  propagation  of  vibration  and  localize  the  vibration  modes.  Depending  on 
the  magnitude  of  perturbations  (disorder)  in  the  individual  components  and  the  strength  of  inter¬ 
nal  coupling  between  the  subsystems,  the  mode  shapes  may  undergo  dramatic  changes  and 
become  strongly  localized  when  small  perturbations  are  introduced,  thereby  confining  the  energy 
associated  with  a  given  mode  to  a  small  geometric  region.  This  phenomenon  is  referred  to  as 
mode  localization.  Pierre  [5]  suggested  that  strong  mode  localization  and  eigenvalue  curve  veer¬ 
ing  are  two  manifestations  of  the  same  phenomenon.  Therefore,  the  investigation  of  the  curves 
of  the  eigenvalues  or  natural  frequencies  in  the  neighborhood  of  the  ordered  state  is  sufficient  for 
determining  the  occurrence  of  strong  mode  localization.  That  eigenloci  veering  phenomenon 
can  occur  in  disordered  structures  under  certain  conditions  has  been  explained  qualitatively 
using  geometric  arguments  in  [6];  where  it  was  also  hinted  that  quantitative  results  might  be 
obtainable  by  the  use  of  singular  perturbation  analysis. 


For  mistuned  linear  chains,  Pierre  [5],  and  Pierre  and  Dowell  [7]  showed  that  the  straightfor- 


-  4  - 


ward  expansion  in  terms  of  mistuning  parameters  breaksdown  in  the  case  of  weak  coupling. 
This  arises  because  the  idealized  system  that  is  being  perturbed  has  natural  frequencies  with 
multiplicity  >  1.  They  then  developed  a  so-called  "modified"  perturbation  technique  which  pro¬ 
vided  a  good  approximation  to  the  exact  eigenfrequencies  and  showed  good  agreement  with 
experimental  results.  In  the  case  of  strong  coupling  between  identical  subsystems,  no  such 
difficulty  arises  and  regular  perturbation  expansions  in  terms  of  mistuning  parameters  are  uni¬ 
formly  valid. 

For  systems  with  cyclic  symmetry  or  spatial  periodicity,  however,  mode  localization  can 
arise  in  the  presence  of  perturbations  which  split  the  degenerate  or  coincident  eigenvalues," 
irrespective  of  the  strength  of  internal  coupling  [6].  Using  differential  topological  ideas  it  was 
shown  qualitatively  in  [6]  that  circularly  configured  systems  which  have  cyclic  symmetry  exhibit 
complicated  topological  behavior  even  for  strong  coupling  when  small  perturbations  are 
imposed.  Furthermore,  the  frequency  response  of  a  perturbed  cyclic  system  depends 
significantly  on  the  form  of  the  perturbation.  Such  cyclic  periodic  systems  are  important  to  the 
analysis  of  vibrations  of  bladed  disk  assemblies. 

In  the  present  work  we  show  that  the  singularity  causing  the  breakdown  of  the  straightfor¬ 
ward  expansion  can  be  analyzed  by  the  well  developed  singular  perturbation  techniques  [8]  and 
appropriate  uniform  asymptotic  expansions  for  the  eigenfrequencies  and  eigenvectors  can  be 
constructed  which  provide  a  correct  qualitative  and  good  quantitative  approximation.  Prelim¬ 
inary  results  on  eigenvalue  veering  for  the  now  standard  example  [5]  of  the  coupled  penduli 
shown  in  Figure  1  were  recently  reported  in  a  short  paper  [9].  Here  we  present  complete  details 
of  the  singular  perturbation  analysis  for  the  eigenvalue  problem  of  the  coupled  penduli  system. 
Using  the  uniform  asymptotic  expansions  for  eigenvectors,  the  occurance  of  mode  localization  is 
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then  related  to  sensitivity  with  respect  to  parameter  variations.  We  then  consider  the  simplest  of 
examples  of  systems  with  cyclic  symmetry  consisting  of  three  identical  masses  arranged  in  a 
ring,  interconnected  by  identical  springs  and  having  individual  torsional  stiffnesses.  The  eigen¬ 
value  veering  is  here  shown  to  exist  even  for  the  strong  coupling  case.  Finally,  based  on  the 
solutions  for  the  strong  coupling  case,  behavior  for  the  weak  coupling  limit  is  explored. 

The  basic  idea  of  analysis  by  the  singular  perturbation  technique  is  the  following:  by  apply¬ 
ing  the  regular  perturbation  technique  to  the  eigenvalue  problem  A0  =  X0,  of  the  system,  we  can 
obtain  algebraic  expressions  for  the  eigenvalues  and  eigenfuctions  as  a  power  series  in  the  small 
parameter  or  perturbation  (say  5).  The  coefficients  of  the  power  series  are  dependent  on  a 7 
second  parameter  e  and  these  expansions  are  valid  for  sufficiently  small  5,  for  all  values  of  £,  so 
long  as  no  singularities  arise.  Singularities  occur  for  values  of  E  where  the  eigenffequencies  and 
eigenfunctions  lose  their  smoothness  and  it  is  said  that  the  expansion  is  not  uniformly  valid  for 
all  e.  Away  from  the  singular  values  of  the  parameter  (e),  the  straightforward  expansions  are 
good  approximations  and  are  called  the  "outer  expansions".  The  neighborhood  of  the  singular 
parameter  point  is  then  stretched  or  rescaled  in  terms  of  a  new  parameter  so  as  to  remove  the 
singularity.  The  expansion  in  terms  of  the  new  parameter  is  valid  only  in  the  neighborhood  of 
the  singular  point  and  is  called  the  "inner  expansion".  The  inner  and  the  outer  solutions  can  be 
matched  where  their  domains  of  validity  overlap  and  then  a  composite  expansion  can  be  con¬ 
structed  which  is  valid  uniformly  throughout  the  function  domain  for  all  values  of  the  parameter 
e. 
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2.  THE  COUPLED  PENDULI 
2.1  Singular  Perturbation  Analysis 

Consider  the  system  consisting  of  two  weakly  coupled  penduli  as  shown  in  Figure  1.  The 
two  important  parameters  are  the  dimensionless  length  A/,  and  the  dimensionless  coupling 
between  the  two  pendulums  8  =  R2  =  (k/m)/(g//).  The  dimensionless  parameter  A /  represents  the 
disorder  or  perturbation  in  the  individual  pendulums.  It  is  important  to  point  out  that  the  dimen¬ 
sionless  coupling  between  the  two  pendulums  8  «  1  for  weak  coupling  irrespective  of  A/.  The 
resulting  eigenvalue  problem  in  symmetric  form  generated  by  the  above  system  is  given  by 

A<J>  =  X<t>,  (1) 

where  6  =  R2  =  C0k/C0g  ,  CDk  =  k/m  ,  coj  =  gfl , 

1+R2  -R2 

A  — 

[-R2  R2+(l+A/r1 

For  small  values  of  5,  it  is  natural  to  expand  eigenvalues  and  eigenfuctions  in  the  regular  expan¬ 
sion  as  powers  of  8  regarding  A/  as  a  parameter  in  the  range  of  interest.  Thus  we  write  A,  X,  and 
<()  in  powers  of  8  as 


A  =  Ao  4-  Aj8  +  A2&2  +  0(53) , 

(2) 

X  =  XQ  +  Xj  8  +  X2  8~  +  0(83 ) , 

(3) 

4)  =  4)o  +<th8  +  02S2  +  0(83 ) . 

(4) 

Substituting  (2),  (3)  and  (4)  into  (1),  equating  coefficients  of  each  power  of  8  to  zero,  and  solv¬ 
ing  the  resulting  sequence  of  homogeneous  and  nonhomogeneous  linear  systems  gives  the  fol¬ 
lowing  expansions  for  the  eigenvalues  and  the  corresponding  eigenvectors 


-7- 


X1  =  l+8  +  (l+A/)/Al]  62  +  0(63),  (5) 

S  N 

X2  =  1/(1  +  A/)  +  5  -  (I  +  M)/M  S2  +  0(53) ,  (6) 

k  - 


In  (7)-(8),  C  and  are  arbitrary  constants.  The  expressions  (5)  -  (8)  are  the  regular  (outer) 
expansions  of  the  eigenvalue  problem  for  small  coupling  5  which  depend  on  the  parameter  A 1. 
These  are  valid  for  sufficiently  small  5  for  all  values  of  A/.  When  A/  — » 0,  (5)  -  (8)  become 
unbounded  and  the  continuity  of  eigenvalues  and  eigenvectors  with  respect  to  the  perturbation  A / 
breaks  down.  Thus,  in  the  neighborhood  of  A/  =  0  the  expansions  (5)-(8 )  become  nonuniform 
and  singular  or  non-analytic  points  have  therefore  been  identified.  These  eigenvalues  in  (5)  and 
(6)  are  plotted  in  Figure  2  for  some  small  but  fixed  5  as  a  function  of  A/.  Each  eigencurve  has 
two  branches,  one  valid  for  A/  >  0  and  the  other  for  A1  <  0.  These  branches  are  indicated  by  the 
superscripts  *+’  and  *-’  and  correspond  to  A/  >  0  and  A/  <  0,  respectively.  Note  that  since  A /  and 
8  have  been  treated  as  two  independent  parameters  there  is  no  control  over  expressions  (5)  -  (8) 
in  the  limiting  process  when  A/  -»  0  and  5  -»  0.  By  stretching  the  neighborhood  of  the  singular 
parameter  value  A/  =  0  and  by  taking  into  consideration  the  nature  of  the  singularity  we  can  find 
an  exact  relation  between  these  two  parameters.  Then  Ai  and  5  become  dependent  in  the  neigh¬ 
borhood  of  the  singular  parameter  value,  called  the  "inner  region".  The  solutions  of  the  problem 
in  the  inner  region  are  called  the  inner  expansions.  Asymptotically  matching  the  inner  and  the 
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outer  expansions,  and  combining  them  appropriately  then  gives  the  composite  expansions  which 
are  valid  throughout  the  interval  of  interest  in  A /  for  sufficiently  small  5.  Before  finding  inner 
expansions  for  the  eigenvalues  and  the  eigenfunctions,  we  mass  normalize  eigenfunctions  01  and 
<f)2  to  get  a  unique  set  of  eigenfunctions 

5  +  0(52) ,  (9; 


5  +  0(52) . 


(10>- 


In  the  inner  expansion,  we  assume  that  the  physical  parameters  R2  =  8  and  A/  =  e  are  related 
(dependent)  by  a  set  of  mathematical  parameters  £j,  ^  ^3,  .  .  .  and  p  through  a  "stretching" 


transformation  of  the  form 


£  =  e0+^a  +  X^(pay ,  (11) 

po 

where  p  is  a  new  small  parameter  that  is  defined  by: 

5(p)  =  (sgn5)pb  .  (12) 

For  fixed  p,  the  quantity  ^  serves  as  the  internal  variable.  The  positive  constants  a  and  b  are  to 
be  determined  by  the  nature  of  the  singularities  of  the  characteristic  equation  F(X,e,8)  =  0  of  (1) 
near  e  =  e0,  where  e0  is  the  singular  point  of  interest.  Let  the  eigenfunctions  Q  and  the  eigen¬ 
values  X  be  written,  in  the  inner  region,  as  the  expansions 
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<J)(e,5)  =  z(e(ji),5<4))  =  z(n)  =  £  ,  (13) 

j=o 

X(e,5)  =  A.(e(n),8(n))  =  O.  (4)  =  £  0^  .  (14) 

j=o 

Note  that  for  the  pendulum  problem  F(k,£,0)  =  0  has  a  singular  point  at  e  =  £c  =  0.  The  expan¬ 
sions  (13)  and  (14)  are  the  "inner  expansions"  and  the  Zj’s  and  ftj’s  are  called  the  "inner 
coefficients".  Substituting  (11),  (12),  (13),  and  (14)  into  (1)  simplifying  the  inner  expansions 
with  a  =  b  =  1,  and  solving  the  sequence  of  eigenvalue  problems  obtained  by  equating  each 
power  of  p.  to  zero,  we  obtain  the  following  inner  eigenvalues  Q1,^2  and  inner  eigenfunctions 
z1  and  z2 


n'  =  i+  „+  77^^77  n2+o(n3), 

-  >  ^ 


H2+0(ii3), 


r~ki 

z1  =  j  +  1  +  k2  H  +  0(42) , 

-*  0 


r  -I  b0^2m2 

■)  ,  m  - —  - 

z  =  b0  i  +  1  +  m2  |4  +  0(p.2) , 


where  A/  =  £  =  £4  ,  5  =  R2  =  ji  , 


V^2+4 


(19) 
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(20) 


Here  the  eigenvectors  z1  and  z2  have  been  mass  normalized.  Keeping  p.  fixed  and  taking  the 
limit  |^|  —>  oo,  it  is  easy  to  see  that  Q1  matches  asymptotically  with  X\  for  ^  —>  °°,  and  with  A.2 
for  \  ->  -oo.  Similarly  the  inner  solution  Q2  matches  asymptotically  ( |^|  -+  «>)  with  Xl  and  A.2. 
In  fact,  it  can  be  easily  shown  that  the  outer  eigenvectors  <J>„  and  0^  match  with  the  inner  eigen¬ 
vectors  z1  and  z2  in  exactly  the  way  the  eigenvalue  branches  match.  The  composite  expansions 
for  the  eigenvalues  and  the  eigenfunctions  are  now  obtained  by  combining  the  inner  and  the 
outer  expansions,  and  subtracting  the  common  part  of  the  two  expansions.  The  resulting  expres¬ 
sions  for  the  eigenvalues  and  the  eigenvectors  are 


A  1 

^comp 


r  y. 

=  T-^rr  +  5-  1  +  T7  52  (1  -u(A/))  +  1  +  6+  1  +  4t  52  u(A/) 
1  +  A  /  A  /  A/ 
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— (1  — u(A/)) (1  -A/ +A/2  +  8-(l+l/A/)82)-(l  +8+(l  +  l/A/)82)u(A/)  +  0(83) , 


^LP=  1+5+  1  +  — 


Jr  52  (1  -u(A/))+  T^-rr  +  5-  1  +  -Jr  S2  u(A /) 

1  +  A/  A/ 

J  J 


r  , _  ^  (A//5)2 

j1+  zHmz$^I±  s+  ,  f(A//5)+^inf  5 
1  J  l  2  J 


-( 1  -  u(A/))(  1  +  5  +  ( 1  +  -!-)52)  -  u(A/)[  1  -A/  +  A/2  +  8-(l  +  -5-)82  ]  +  0(63 ) , 

A  /  A/ 


-1  +  8/2  A/ 

0comP  =  z1  +  uCA/)^  +  (l-uCA/))^  -  u(A/)  g/A/  _  83 

2  A/ 3 


-fl-u(A/) 


A/  1  +  ■ 


c2 

5(1  +  — — -T-) 
2  A/ 2 


5  +0(52), 


which  simplifies  to 


Lkl  k2A/2/82 

4*  comp  —  aoj  ,  +  1+k2  8  +  0(82) , 
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and 


4> 


2 

comp 


m2A/2/52 
1  +  m2 


+  0(82) , 


(26) 


where  k,  m,  a<j,  and  bQ  are  already  defined  in  equations  (19),  (20),  (21)  and  (22)  respectively  and, 


u(A /)  = 


fl  ,  A/  >  0 
jo  ,  A/  <  0  ' 


The  composite  solutions  (23)  -  (26)  of  the  eigenvalue  problem  (1)  are  the  asymptotic  approxima-  - 
tions  for  small  coupling  5,  and  are  uniformly  valid  for  all  mistunings  A/. 

The  eigenfrequencies  A.c0mp  and  X^mp,  obtained  by  the  singular  perturbation  analysis,  are 
plotted  in  Figure  3  as  a  function  of  the  parameter  A/.  Since  the  coupled  penduli  system  is  sim¬ 
ple,  exact  expressions  for  the  eigenfrequencies  are  easily  obtained  and  they  are  also  plotted  in 
the  figure.  The  exact  solutions  for  the  eigenvalue  problem  are  given  in  the  Appendix.  Clearly, 
there  is  an  excellent  agreement  in  the  exact  frequencies  and  their  asymptotic  approximations. 
The  asymptotic  solutions  also  clearly  display  the  veering  phenomenon. 


We  now  study  the  behavior  of  eigenvectors  for  the  case  of  weak  coupling  when  A /  -  0. 


2.2  Eigenvector  Rotations  and  the  Sensitivity  Function 

In  our  earlier  work  [6]  with  mistuned  cyclic  systems,  it  was  suggested  that  localization  of 
modes  can  be  investigated  by  considering  the  sensitivity  function,  and  the  rotations  of  eigenvec¬ 
tors  under  variations  of  parameters.  The  sensitivity  function  of  eigenvectors  or  eigenvector  sen- 


-  13- 


sitivity,  in  short,  is  defined  as 

1 1  Su  1 1  =  *\/tr(SuSu)  (Frobenius  norm) , 

where  Su  =  uj^Au,  Su  denotes  the  complex  conjugate  of  transpose  of  Su,  uQ  is  the  modal  matrix 
of  eigenvectors  for  zero  coupling  (5  =  0),  ‘tr’  denotes  the  trace  of  the  matrix,  and  Au  =  u  -  uQ 
where  u  is  the  modal  matrix  for  nonzero  coupling  (5  *  0).  The  eigenvector  sensitivity  evaluated 
for  the  pendulum  problem  turns  out  to  be 

I  I  Su  |  |  =  V(“Pl  1  +  P21  )2  +  (Pi  1  +  P21  )2  +  (-Pl2  +  P22)2  +  (Pl2  +  P22)2  ,  (27) 

v2 

where 


!  aok  2(A1)2 

Pi  1  =  —aok  +  — —  +  — - r— 

V2  6(1 +  k2) 


1  b0m(A /) 

P12  =  b0m - ~  + 


V5"  6(1  +  m2) 


P2I  =  3o  ~  -p-  , 

v2 

P22=b„-^=. 

V2 

We  can  also  define  the  angles  between  the  eigenvectors  $1  and  02  for  the  unperturbed  (5  =  0) 
system,  and  the  eigenvectors  fylomp  and  <t>comp  for  the  perturbed  system.  These  angles  are  given 
by 


<<j>0  .  4* comp  >  1  -  Pi 

II  4>o  II II  ^comp  II  ’ 


(28) 
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COS  02  = 


<$0  »  $comp  ->  _  1^*2 


llflll  llOimpll 


■Jlyfl  +  P< 


where 


Pi 


kW_ 

8(l+k2)  ’ 


P2  =  m  + 


m 2  (A/)2 
6(l+m2)  ' 


Plots  of  the  sensitivity  function  of  eigenvectors  and  the  cosine  of  the  angle  between  nominal 
(5  =  0)  and  perturbed  eigenvectors,  as  a  function  of  the  mistuning  A/,  are  given  in  Figures  4  and 
5,  for  both  the  exact  solutions  and  the  asymptotic  approximations  obtained  above  by  the  singular 
perturbation  technique.  These  solutions  are  in  excellent  agreement.  Figures  4  and  5  clearly 
show  and  confirm  the  expectation  that  the  eigenvectors  for  the  weakly  coupled  system  undergo 
rapid  changes  in  the  vicinity  of  the  singular  point.  Furthermore,  either  of  the  two  criterion  can 
be  effectively  used  as  a  quantitative  measure  and  indicator  of  the  mode  localization 
phenomenon. 


3  CYCLIC  SYSTEMS 

Consider  three  identical  particles,  each  of  mass  m  arranged  in  a  ring  and  interconnected  by 
identical  springs  of  stiffness  kc.  Assume  that  all  the  masses  are  hinged  to  the  ground  by  torsional 
springs  of  stiffness  kt  and  that  the  radius  of  the  ring  is  r,  as  shown  in  Figure  6.  As  a  perturbed 
system  we  consider  the  case  when  two  of  the  torsional  springs  are  perturbed  by  £i  and  e2.  The 
eigenvalue  problem  corresponding  to  this  system  is  given  by 
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a  +  £i  -oil  -to l 

-co l  a+£2  -<Oc  <j>  =  X<j>,  (29) 

-al  -Oil  a 

2kc  kt  ,  2 

where  a  = - +  — r-  and  0i£  =  —  • 

m  mr2  m 

As  was  the  case  for  the  example  of  linear  chains,  the  interest  here  is  in  the  development  of 
asymptotic  expansions  for  the  eigenvalues  and  eigenvectors  in  terms  of  the  perturbation  parame¬ 
ters  Ei  and  £2.  As  we  shall  see  the  unperturbed  cyclic  system  (ei  =  £2  =  0),  has  a  double  eigen¬ 
value  and  an  isolated  eigenvalue,  and  thus,  introduction  of  perturbations  is  expected  to  split  the 
degenerate  eigenvalue  pair.  Expecting  that  the  eigenvector  behavior  will  be  governed  by  the 
eigenvalue  behavior,  similar  to  the  case  of  the  coupled  penduli,  we  restrict  the  developments  to 
expansions  for  only  the  eigenvalues  of  the  perturbed  three  particle  system. 

The  eigenvalue  problem  (29)  results  in  the  characteristic  equation 

F(X,£!  ,£2) =  — (3a  +  £i  +  £2)X2  +  (3a2  +  2a£2  +2a£2  +£j  £2  —  3c0c)X. 

+  2c0c +o)c(3a  +  £i +£2)-(a2 +  a£1)(a  +  £2)j  =0.  (30) 

First  consider  the  unperturbed  cyclic  system  the  eigenvalues  or  roots  of  (30)  are  given  by 

r-1  r-2  r-3  2 

X=A=a  +  o>c,A.=a  -  2o)c  . 

The  corresponding  roots  of  (31)  are 

—1  —2  2  _  2 

X  =X  =<0c  ,  X  =  -2co c  • 

Thus,  there  is  a  coincident  pair  of  eigenvalues  and  one  isolated  eigenvalue.  So  long  as 
oil  ~  0(1),  the  two  distinct  eigenvalue  are  well  separated. 
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To  study  the  perturbed  problem  we  first  introduce  the  coordinate  transformation  x  =  X  -  a,  so 
that  equation  (30)  results  in 

X3  -  (ei+e2)x2  +  (£i£2  ~  3c o*)x  +  +  e+e2)  =  0  (31) 

We  can  express  the  solutions  to  (31)  as  regular  functions  of  the  parameters  ej,  e2,  and  co3  by 

oo 

writing  x(ei-e»)  =  SXj(£i)ei  as  a  power  series  in  e2.  Substituting  the  resulting  expression  in 
j=o 

(31)  and  proceeding  in  the  usual  manner,  the  expansions  for  the  three  roots  of  (31)  turn  out  to  be 


i  i  £2  1  C0c~ ,,  , 

x1  =co3+  +  el+0(ei), 

2  8  a)‘£! 

k.  * 

(32) 

X2  =  X()2  +  Xl2£2  +  X22£2  +  0(£^)  , 

(33) 

X3  =  X03  +  Xl3£2  +  X23£2  +  0(£2)  » 

(34) 

where 

£i  -  w  1+  Pi  ^  _  X02-X02£l  ~  ©c 

X02  >  Xl2  -5  2t  t  4  ’ 

2  3xo2-2Xo2£i  -  3g)c 

Xl2(£l“3X02)  +  Xl2(2X02_el) 

112  3X  02-2X02£1  -  3o)c 

E^col-Pj  ^  _  X03-X03£l  -  ©c 

X03  —  *  X 13  —  2  -i  -5  4’ 

2  3x  53  -2X03^1  -  3 coj 

_  X 13 (el  — 3X02)  +  Xl3(2X03-£l) 

X23  3x&3-2xo3£i  -  3c0c 

and 
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Pl  =  ^9(0%  +  2ci)f£1  +  El  . 

The  expressions  (32)  -  (34)  are  the  straightforward  expansions  of  the  eigenvalue  problem  for 
small  £2-  When  £i  -»0,  and  X2  become  unbounded  and  the  continuity  of  the  eigenvalues 
with  respect  to  the  parameter  or  perturbation  £\  breaks  down.  The  third  eigenvalue  x3  always 
remain  bounded  and  continuity  is  preserved  for  all  values  of  so  long  as  the  interconnecting  or 
coupling  spring  constant  is  0(1).  These  eigenvalues  (32)-(34)  are  shown  in  Figure  7.  Clearly 
£]  =  0,  £?  =  0  is  a  singular  point  of  the  expansions  and  (32)-(34)  are  the  outer  expansions  valid 
for  small  £2  away  from  £1  =  0. 

Inner  expansions,  which  are  valid  in  the  neighborhood  of  singular  point  or  the  parameter 
values  where  the  outer  expansions  breakdown  are  now  obtained.  The  expansion  process  for  the 
cyclic  system  is  very  similar  to  the  one  presented  in  section  2  for  the  linear  chain  system.  Thus, 
the  physical  parameter  perturbations  £t  and  £2  are  related  to  a  set  of  parameters  ^ ,  £,2>  ^3 >  and  |i 
via  a  "stretching"  transformation  of  the  form 

£  =  ^a  +  £^(4ay  ,  (35) 

j=o 

where  u  is  a  new  small  parameters  that  is  defined  by 

5((i)  =  (sgn  6)nb  .  (36) 

The  positive  constants  a  and  b  are  to  be  determined  by  the  nature  of  the  singularities  of 

F(X,£i  ,£2)  =  0  near  £j  =0,  the  singular  point  of  interest.  Let  the  dependent  variable  or  the 

eigenvalue  £2(p.)  =  Q(Ei(ji),  £2(4))  be  written  as  the  expansion 

ft  =  £nj  M.J  • 


(37) 
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Note  that  for  the  cyclic  system  F(X,£i,£2  =  0)  has  a  singular  point  at  £t  =  e2  =  0.  The  expan¬ 
sions  (37)  are  called  the  "inner  expressions”  and  the  ftjs  are  called  the  "inner  coefficients".  Sub¬ 
stituting  (35),  (36)  and  (37)  into  (30),  simplifying  the  inner  expansions  with  a=b=l,  and  perform¬ 
ing  the  perturbation  analysis,  the  following  roots  of  (30)  are  obtained 

,  (£+l)Qn  - Qii-^Qn 

ft  =  co2  +  ftn  P  +  - 7  /  ==  -  P2  +  0(p3),  (38) 

2a)?V^+  1 

»  - 


ft2  =  co2  +  012  p  + 


fll2+4fli2~(^+l)Qi2 

20)c2^2-^+1 


+  0(p3). 


ft3  =  -2c o2  + 


(S2-^+l)p2+0(p3), 


where 


(40) 


«,  -fr  .  *-n.  ^  n12,  . 


Keeping  p  fixed  and  taking  the  limit  |  — > «»,  it  is  easy  to  see  that  ft1  matches  asymptotically 

with  xl  for  ^  -»  00  and  with  X-  for  £  -» -00.  Similarly  the  inner  solution  Q2  matches  asymptoti¬ 
cally  (|^|  -»  00)  with  xl  and  %-•  As  expected,  the  third  root  ft3  automatically  matches  with  x3 
as  no  singular  behavior  is  displayed  in  this  case.  Now  combining  the  inner  and  the  outer  expan¬ 
sions  appropriately,  we  get  the  composite  expansions 

xLnp  =  X2(u(ei))  +  XI(1-u(ei))  +  ft1  -  u(et )  (common  pans  of  x2  ,  ft1) 


-  ( l-u(£t ))  (common  pans  of  x1  ,  ft1 ) , 
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Xremp  =  X1^!))  +  Xz(l-u(£i))  +  &  -  u(£!)  (common  pans  of  xl  .  ***) 


1  r»2s 


-  (l-utejlXcommon  pans  of  %2,Q2) , 


Xcomp  ~  X  > 


where  u(£i )  = 


1  ,  £i  ^0 

0  ,  et  <0  ’ 


Transforming  these  expansions  back  into  the  original  coordinates  gives  the  following  composite 
expansions  for  eigenvalues 


^-comp  —  a+a)c  f^ll  ^2 


(1  +  ei/E2)^  1 1  ~  ^11  ~  (£i/e2)Qii 

+  /  — . - . .  e2  +  °(el) » 

2coaV(e1/e2)-(ei/e2)+  1 


^•comp  —  a+Ct)c  +  Qi2  £2  + 


Q12  +  Ei/e2Qi2  -  (6j/£2  +  l)fti2 

20 )cV(£l/£2)2  ~(el/e2)+  1 


el  +  0(el). 


^-comp  =  X03  +  (Xi3)£2  +  X23(ei)  +  O(eJ) , 


where 


f2n  — 


(ej/£2  +  1)  +  \(£i/£2)  -(Ei/£2)+1 


= 


(£t/£2  +  1)  ~  V(£i/e2)2  -  (Ci/e2)  +  1 
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=  X03  ~X03et  ~  mc 

3x^3  “  2X03^1  ~  3t»)2 

£i  -  ©2  -  Pi 

X03  - - 2  ’ 

Pt  =  V9C0c  +  200^!  +£?  - 

The  plots  of  eigenfrequencies  Xc0mp,  Xc0mp  and  Xcomp  as  a  function  of  £i  are  given  in  Figure  8a 
and  8b.  The  third  eigenvalue  X^mp  always  behaves  as  a  regular  function  as  is  clear  from  Figure 
8a.  These  results  show  that  for  the  cyclic  system  with  strong  coupling  the  effect  of  perturbation 
is  identical  to  that  for  the  pendulum  problem  with  weak  coupling.  That  is,  the  curve  veering  of 
the  eigenvalues  X*omp  and  X^omp  is  almost  the  same  as  the  curve  veering  of  the  weakly  coupled 
pendulum.  Consequently,  it  is  expected  that  the  mode  localization  of  the  eigenvectors  within  the 
strongly  coupled  cyclic  system  be  similar  to  that  in  the  weakly  coupled  pendulum  problem. 
Therefore,  the  results  should  be  obtained  for  the  singular  behavior  of  the  eigenvectors  of  the 
weakly  coupled  pendulum  problem  in  section  2  should  be  valid  quatitatively  for  the  cyclic  sys¬ 
tem. 

The  coupling  constant  Iq.  plays  a  very  important  role  in  the  eigenvalue  veering  behavior.  For 
the  cyclic  system,  kc  =  0  leads  to  three  coincident  eigenvalues  as  opposed  to  the  case  of  strong 
coupling  when  only  a  double  eigenvalue  appears.  The  composite  expansion*  obtained  earlier 
were  determined  under  the  assumption  that  —  0(1).  We  now  use  these  expansions  and 
explore  the  behavior  of  eigenvalues  in  the  limiting  ease  of  kc  — *  0.  Note  that  the  coupling  con¬ 
stant  ©c  appcsu-  in  the  denominator  of  the  asymptotic  expansions  (41)-(43).  As  ©c  is  reduced 
the  isolated  eigenvalue  Xc0mp  rnoves  closer  to  and  X^,,^  variation  of  X.  s  with  £2 
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becomes  very  rapid,  as  is  evident  from  Figures  9  and  10.  Figure  11  shows  the  three  composite 
eigenvalue  expansions  for  small  £.\  and  £ 2  as  a  function  of  kc  (or  C0c )  and  the  expansions  clearly 
breakdown  in  the  vicinity  of  kc  =0.  Thus,  the  composite  expansions  (41)-(43)  are  behaving  as 
outer  expansions  in  the  weak  coupling  limit.  It  should  be  possible  to  now  construct  uniformly 
valid  expansions  for  eigenvalues  as  a  function  of  the  coupling  constant  Iq.  by  one  more  use  of  the 
singular  perturbation  technique  whereby  the  neighborhood  of  lq  =  0  is  stretched  and  an  inner 
expansion  is  developed.  The  localization  and  veering  behavior  of  thus  obtained  composite 
eigenvalue  expansions,  which  will  be  valid  for  small  enough  £2  for  all  £1  and  lq,  is  expected  to 
be  much  more  interesting  and  is  being  presently  studied. 

4.  SUMMARY  AND  CONCLUSIONS 

Singular  perturbation  technique  has  been  applied  to  two  parameter  eigenvalue  problems  to 
obtain  uniformly  valid  algebraic  expansions  for  the  eigenvalues  and  the  eigenvectors  for  two 
example  systems.  Utilizing  these  expansions,  eigenloci  veering  and  the  mode  localization 
phenomenon  have  been  studied.  A  sensitivity  function  and  the  rotation  of  eigenvectors  have 
been  introduced  as  criteria  to  visualize  mode  localization  phenomenon  in  the  vicinity  of  singular 
points.  One  example,  that  of  the  two  weakly  coupled  penduli,  represents  systems  of  the  linear 
chain  type  with  only  one  weak  coupling  spring.  The  example  of  three  mass  particles  belongs  to 
strongly  coupled  systems  with  cyclic  symmetry.  For  the  coupled  penduli  system,  the  eigenvalue 
and  eigenvector  expansions  are  found  to  be  in  excellend  agreement  with  the  exact  results. 

It  is  shown  that  eigenvalue  curve  veering  occurs  both  in  the  weakly  coupled  penduli  and  the 
strongly  coupled  cyclic  system.  The  effects  of  mistuning  perturbations  which  split  a  pair  of 
coincident  eigenvalues  is  identical  in  both  the  cases.  The  eigenvector  sensitivity  function  and 
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the  angle  of  rotation  of  eigenvectors  are  shown  to  be  two  equally  good  candidates  for  visualizing 
mode  localization  phenomenon  near  singular  points.  The  composite  expansions  for  the  per¬ 
turbed  cyclic  system,  which  are  uniformly  valid  in  the  case  of  strong  coupling,  are  shown  to 
breakdown  in  the  limiting  case  when  the  coupling  stiffness  goes  to  zero.  This  clearly  is  related 
to  the  fact  that  all  the  eigenvalues  for  the  cyclic  system  in  the  weak  coupling  limit  are  identical. 
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APPENDIX 

It  can  be  easily  shown  that  the  two  free  vibration  natural  frequencies  to2  and  the  correspond¬ 
ing  mass  normalized  eigenvectors  x1  for  the  coupled  pendula  problem,  obtained  from  the  exact 
solution  of  the  eigenvalue  problem,  are  given  by 


2 cot  ?  —  2R2  +  1  + - + 

1,2  1+A/ 


4R4  +  1 - —  + 
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x  =  a 
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The  sensitivity  of  eigenvectors  is  then  given  by 


ISJII  =  — j=-  'V(— <ii i -+<121  )2  +  (qn+^i)2  +  +  (qi2+ti22)2  » 

v2 


where 


l  1 

Qn  -  Ot  Xn  +  — p- . 

v2 


2  1 


qi2  -Pxn  - 


Vi" 


q:i  =  <x- 


VT 


q:i  =  P  - 


Vi" 


The  angles  between  the  nominal  and  the  perturbed  eigenvectors  are  then 
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Figure  7.  Outer  expansions  for  the  eigenvalues  of  the  perturbed  cyclic  system  indicating  the 
region  of  singular  behavior;  k<.  =  2,  kt  =  1,  e2  =  0.1,  r  =  1,  m  =  1. 
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Figure  8a.  Eigenvalues  from  composite  expansions  for  the  perturbed  cy< 
kt  =  1,  £2  =  0.1,  m  =  1,  r  =  1. 
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I.  Introduction 

An  important  subject  in  the  dynamics  and  control  of  structural  systems  is  the 
behavior  of  structures  under  transient  or  steady  state  excitations.  In  this  work,  we 
examine  the  stability  of  the  geometric  form  of  the  spatial  configuration  of  struc¬ 
tural  systems  when  the  structural  parameters  are  subject  to  small  perturbations, 
and  the  implications  of  this  instability  for  frequency  response.  We  show  that  cir¬ 
cularly  configured  systems  which  nominally  have  cyclic  symmetry  exhibit  com¬ 
plicated  topological  behavior  when  small  perturbations  are  impressed  on  them. 
We  further  show  that  the  frequency  response  of  a  perturbed  cyclic  system 
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depends  considerably  on  the  form  of  perturbation.  On  the  other  hand,  a  rectilinear 
configuration  of  nearly  identical  subsystems  does  not  ehibit  modal  instability. 
Usually,  both  kinds  of  systems  are  implicitly  assumed  to  undergo  similar  qualita¬ 
tive  behavior  under  a  small  perturbation  whereas,  in  fact,  the  cyclic  configuration 
exhibits  a  very  stange  behavior,  [1], 

The  distinction  between  the  behavior  of  cyclic  and  rectilinear  configurations 
under  a  perturbation  is  important  because  many  engineering  structures  are  com¬ 
posed  of  identical  substructures  which  are  replicated  either  in  a  uni-axial  chain, 
or  in  a  closed  cyclic  formation  where  modal  control  is  of  interest.  Examples  of  the 
former  case  of  periodicity  occur  in  structures  such  as  space  platforms  and  bridges, 
which  have  an  obvious  periodicity  of  the  uni-axial  kind.  An  example  of  cyclic 
periodic  systems  is  a  turbine  rotor,  which  consists  of  a  set  of  nominally  identical 
blades  mounted  on  a  central  hub,  and  often  referred  to  as  a  “bladed  disk  assem¬ 
bly”  [2].  When  all  the  blades  are  truly  identical,  then  the  system  is  referred  to  in 
the  literature  as  a  tuned  bladed  disk  assembly.  Practical  realities  of  manufacturing 
processes  preclude  the  existence  of  exact  uniformity  among  all  the  blades.  When 
residual  differences  from  one  blade  to  another — no  matter  how  small — are 
accounted  for  in  the  theoretical  model,  the  assembly  is  then  termed  a  mistuned 
bladed  disk. 

Our  primary  focus  in  this  investigation  is  on  bladed  disk  assemblies.  However, 
since  we  approach  the  problem  from  a  generalized  viewpoint,  the  conclusions  to 
be  drawn  will  be  of  relevance  to  other  periodic  systems.  Therefore,  in  the  sequel, 
we  borrow  the  ‘tuned’  and  ‘mistuned’  terminology  from  the  bladed  disk  litera¬ 
ture,  and  apply  it  to  repetitive  systems  having  cyclic  or  uniaxial  periodicity.  Thus, 
in  a  tuned  periodic  system,  the  nominal  periodicity  is  preserved,  whereas  it  is  des¬ 
troyed  in  a  mistuned  system. 

If  we  examine  the  system  matrices  of  the  linear  and  cyclic  chains,  we  observe 
a  fundamental  difference  in  forms.  The  dynamical  matrix  of  the  linear  chain  is 
usually  banded.  Banded  matrices  are  frequently  encountered  in  structural  dynam¬ 
ics.  A  special  form  of  banded  matrices  that  is  of  interest  here  is  the  tri-diagonal 
form  a,j=Q,\i-j\>  1.  On  the  other  hand,  the  system  matrix  of  a  cyclic  chain  has  a 
circulani  submatrix,  or  is  entirely  circulant  or  block  circulant  [3],  Circulant 
matrices  usually  arise  in  the  study  of  circular  systems.  They  have  interesting  pro¬ 
perties  that  set  them  apart  from  matrices  of  other  forms  [4],  We  note  that  all 


MODAL  STABILITY  OF  IMPERFECT  CYCLIC  SYSTEMS 


139 


circulants  commute  under  multiplication,  and  are  diagonalizable  by  the  fourier 
matrix.  One  of  the  most  important  consequences  of  the  foregoing  is  that  the  cyclic 
chain  has  a  series  of  degenerate  eigenvalues,  whereas  the  eigenvalues  of  the  uni¬ 
axial  chain  are  all  simple. 

We  know  that  a  tuned  circulant  matrix,  having  a  multitude  of  degenerate 
eigenvalues,  lies  on  a  bifurcation  set  [5].  Thus,  the  reduction  of  such  matrices  to 
Jordan  normal  form  is  an  unstable  operation  [6].  Consequently,  if  a  non-singular 
deformation  due  to  rrustumng  is  applied  to  a  circulant  matrix,  then  some  of  the 
eigenvectors  will  undergo  rapid  re-alignment,  if  the  mistiming  leads  to  a  crossing 
of  the  bifurcation  set.  If  however,  no  crossing  of  the  bifurcation  set  takes  place, 
then  the  tuned  system’s  eigenvectors  will  be  very  stable,  preserving  their  align¬ 
ment  under  small  perturbations.  In  contrast,  the  eigenvectors  of  a  tuned  banded 
matrix,  being  analytically  dependent  on  parameters,  are  not  generally  disoriented 
by  mistuning  until  the  eigenvalues  are  pathologically  close  [7]. 

If  one  examines  the  literature  in  structural  dynamics,  it  is  observed  that  some 
unusual  behavior  has  been  reported  in  the  study  of  perturbed  cyclic  systems.  This 
has  been  the  case  in  various  studies  of  rings  [8],  circular  saws  [9],  and  other 
cyclic  structures  [10],  But  that  such  anomalous  behavior  is  due  to  a  “geometric 
instability”  inherent  in  the  cyclicity  of  the  tuned  system  has  not  been  previously 
established  in  the  literature,  to  our  knowledge.  Indeed,  it  is  often  assumed  (see, 
for  instance,  [11])  that  the  linear  and  cyclic  chains  would  undergo  the  same  quali¬ 
tative  behavior  under  slight  parameter  perturbations  so  that  small  order  perturba¬ 
tions  of  the  system  matrix  will  lead  to  no  more  than  small  order  differences  in  the 
system  response  relative  to  the  unperturbed  case,  if  the  system  has  “strong  cou¬ 
pling". 

In  this  paper,  we  show  that  such  an  assumption  regarding  qualitative  behavior 
does  not  actually  hold  in  the  case  of  cyclic  systems;  that  cyclic  systems  exhibit  a 
peculiarity  of  their  own  under  parameter  perturbation;  that,  although  a  certain 
amount  of  mistuning  may  produce  little  difference  relative  to  the  tuned  datum  in 
one  case,  a  considerable  change  could  be  induced  if  a  slightly  different  kind  of 
mistuning  is  applied  to  the  same  cyclic  system;  that  such  apparently  erratic 
behavior  arises  in  cyclic  system,  even  when  the  system  has  “strong”  coupling.  In 
carrying  out  this  work,  we  borrow  from  certain  developments  in  differential  topol¬ 
ogy  specifically,  from  Arnold’s  monumental  work  in  singularity  theory  [6,  12-16], 
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II.  Topological  Dynamics  of  Quadratic  Systems 

In  mistuned  dynamical  systems,  a  major  concern  is  to  understand  which 
specific  kinds  of  mistuning  parameters,  or  combinations  thereof,  lead  to  unaccept¬ 
ably  high  amplitude  ratios.  In  this  section,  we  give  an  indication  of  the  taxonomy 
of  the  different  consequences  of  mistuning  in  the  hope  of  isolating  those  that  lead 
to  high  ratios. 

Consider  a  mechanical  system  under  small  oscillations  with  kinetic  and  poten¬ 
tial  energies  given  by: 

7  = /)  jc  *  .V/x  >  0.  ('  =  Zix*  Kx  >  0;  x,  x*0.  (2.1) 

Under  the  influence  of  a  forcing  function  /  (r),  (2.1)  produces  the  following  equa¬ 
tions  of  motion  by  application  of  Lagrange’s  formula: 

Vfx  +  Kx  =  /;  x./e  <T  (2.2) 

where  M  and  K  are  symmetric  •'.d  positive  definite.  A  theorem  in  linear 
algebra  shows  that  there  exists  some  non-singular  transformation  matrix  P  such 
that: 

PT\tP  =  /,  and  PtKP  =  A  (2.3) 

where  A  is  a  diagonal  matrix  of  eigenvalues  whose  elements  satisfy  the  equation: 

det(M  -  \K)  =  0  (2.4) 

Consequently,  by  putting 

x  =  Pq ,  (2.5) 

substituting  for  q  in  (2.1),  and  premultiplying  every  term  of  the  resultant  equation 
by  PT ,  we  obtain  a  new  equation  set: 

q  +  Aq  =  /',  (2.6) 

where  /'  =  PTf.  Hence: 

<7.  +  \q>  =/',  for  i  =  1,  2,  •  •  • ,  n.  (2.1) 

Systems  which  can  be  reduced  to  the  above  form  are  called  quadratic  systems. 
They  are  called  quadratic  cyclic  systems  if,  in  addition,  V/  and  K  are  cyclic  or 
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circulant  matrices.  Our  basic  aim  is  to  determine  the  nature  of  the  changes  in  the 
dynamical  properties  of  a  quadratic  system  of  a  given  order,  under  random 
differencial  perturbations  in  M  and/or  K.  Central  to  this  investigation  are  the  topo¬ 
logical  concepts  of  structural  stability  and  genericity. 

Let  N  be  a  set  with  a  topology  and  an  equivalence  relation  e.  An  element  x  s  N 
is  stable  (relative  to  e)  if  the  e-equivalence  class  of  x  contains  a  neighborhood  of 
x. 

A  property  P  of  elements  of  N  is  generic  if  the  set  of  all  x  e  w  satisfying  P 
contains  a  subset  A  which  is  a  countable  intersection  of  open  dense  sets  [17]. 

Genericity  is  important  n  our  context  because  a  generic  system  will  in  effect 
display  a  “typical”  behavior.  More  concretely  if  a  given  generic  system  gives  a 
certain  frequency  response,  all  systems  produced  by  differential  parameter  pertur¬ 
bations  about  the  nominal  system  will  also  produce  frequency  response  curves 
that  are  not  only  slight  perturbations  of  the  original  nominal  response  but  also 
geometrically  (isomorphic)  equivalent  to  it.  Such  systems  are  called  versa!  defor¬ 
mations  of  the  nominal  system  [14],  A  versal  deformation  of  a  system  is  a  normal 
form  to  which  it  is  possible  to  reduce  not  only  a  suitable  representation  of  a  nomi¬ 
nal  system,  but  also  the  representation  of  all  nearby  systems  such  that  the  reduc¬ 
tion  transformation  depends  smoothly  on  parameters.  The  key  to  establishing  ver- 
sality,  and  hence  genericity,  is  the  topological  concept  of  transversality. 

Let  iVc.W  be  a  smooth  submanifold  of  the  manifold  Vf.  Consider  a  smooth 
mapping  f-T  ->  M  of  the  parameter  space  F  into  M\  and  let  p  be  a  point  in  F  such 
that  /  (p.)  e  N. 

The  mapping  /  is  transversal  to  N  at  p  if  the  tangent  space  to  M  at  /  (p)  is  the 
sum: 


TMfiiD  =/•  TF^  +  TNf(^ 

Consequently,  two  manifolds  intersect  transversally  if  either  they  do  not  intersect 
at  all  or  intersect  properly  such  that  perturbations  of  the  manifolds  will  neither 
remove  the  intersection  nor  alter  the  type  of  intersection. 

Lemma  2.1,  see  ref  [14], 

A  deformation  f(p)  is  versal  if  and  only  if  the  mapping  f:T  — >  M  is  transversal 
to  the  orbit  off  at  p  =  0. 
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The  abov*  result  is  crucially  important  because: 

(i)  It  classifies  from  the  set  of  all  perturbations  of  a  given  nominal  system, 
those  that  do  not  lead  to  radically  different  dynamical  properties  from  the 
nominal. 

(ii)  It  separates  the  “good"  from  the  “bad”  perturbations  and  hence  enables  us 
to  concentrate  our  study  on  the  bad  perturbations.  Let  Q  denote  the  family 
of  all  real  quadratic  systems  in  R\  The  set  Q  has  the  structure  of  a  vector 
space  of  dimension  Zi[n[n  1 1).  It  can  be  shown  that  Q  also  has  the  struc¬ 
ture  of  a  differentiable  manifold  {13]. 

Let  Qv  denote  the  set  of  quadratic  systems  having  v2  eigenvalues  of 
multiplicity  2,  V3  eigenvalues  of  multiplicity  3  etc.  Qv  is  called  the  degen¬ 
erate  subfamily  of  Q. 


Theorem  2.1,  see  ref  [13]. 

The  transformation  f:F  — >  Q  is  transversal  to  Qv. 

Consequently,  a  generic  family  of  quadratic  sys,..ns  of  a  given  order  is  given 
by  a  transformation,  /,  of  the  space  of  parameters  T  into  the  space  of  all  quadratic 
systems  Q,  such  that  /  is  transversal  to  the  space  of  all  degenerate  quadratic  sys¬ 
tems  Qv. 

Hence  Qv  is  the  degenerate  (bad)  set  and  QlQv  :s  the  generic  set.  Observe 
that  Q  iQ  v  and  Q  v  are  transversal.  Consequendy,  the  fundamental  group  of  the 
space  of  generic  real  quadratic  systems  is  isomorphic  to  the  manifold  of  systems 
without  degenerate  eigenvalues. 

The  above  discussion  leads  inevitably  to  the  following  conclusions: 

(i)  Radical  changes  in  the  dynamical  properties  of  a  nominal  system  occurs 
under  perturbations,  when  the  perturbations  take  the  system  across  the 
boundary  from  Q/Qv  to  Qv  and  vice-versa. 

(ii)  Qv  is  a  smooth  semi- algebraic  submanifold  of  Q ,  and  can  therefore  be 
stratified  into  distinct  fiber  bundles  [14].  By  a  bundle,  we  mean  the  set  of 
all  systems  which  differ  only  by  the  exact  values  of  their  eigenvalues;  but 
for  which  the  number  of  distinct  eigenvalues  as  well  as  the  respective 
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orders  of  the  degenerate  eigenvalues  are  the  same.  Within  the  degenerate 
set,  Qv,  the  crossing  from  one  bundle  to  another  can  also  lead  to  radical 
dynamical  changes.  Each  bundle  is  represented  by  a  specific  Jordan  block 
of  a  certain  order.  Note  that  each  bundle  is  also  transversal  to  Q. 


Theorem  2.2,  ref  [14], 

Qv  is  a  finite  union  of  smooth  sub-manifolds  with  codimension  satisfying 
Codim  Qv  >  2. 

Theorem  2.2  has  the  following  implications: 

(i)  Q'Qv  is  topologically  path  connected.  This  means  that  by  smooth  parame¬ 
ter  variations,  provided  that  the  number  of  variable  parameters  is  less  than 
the  codimension  of  Q  v,  it  is  possible  to  smoothly  pass  from  one  member  of 
Q/Qv  to  another  without  reaching  any  singularity;  that  is,  without 
encountering  any  member  of  Qv.  Such  parameter  variations  will  typically 
not  lead  to  radical  dynamical  changes  in  QiQv. 

(ii)  Because  codim  Qv  >  2.  it  follows  that  a  generic  one-parameter  family  of 
quadratic  systems  cannot  contain  any  degenerate  subfamilies.  Therefore 
under  one-parameter  deformations  of  a  generic  family,  some  eigenvalue 
pairs  may  approach  each  other  but  cannot  be  coincident  (i.e.  cannot  col¬ 
lide).  After  approaching  each  other,  they  must  veer  away  rapidly,  giving 
rise  to  the  so-called  eigenvalue  loci-veering  phenomenon  [18],  under  one- 
parameter  deformations  of  generic  families.  This  offers  a  theoretical  expla¬ 
nation  for  the  eigenloci  veering  phenomenon  which  has  been  observed  in 
perturbed  periodic  systems  without  a  corresponding  phenomenological 
base  [18,  19].  Furthermore,  this  phenomenon  holds  provided  the  system 
has  a  quadratic  structure,  irrespective  of  whether  the  model  arose  from  a 
continuous  or  discrete  structural  system  [20]. 

This  rapid  eigenloci  veering  can,  under  the  right  conditions,  produce  the  mode 
localization  phenomenon  [18].  Since  the  dynamical  properties  of  any  linear 
constant-coefficient  system  are  totally  determined  by  its  eigen-structure  (eigen¬ 
values  and  eigenvectors),  and  since  the  eigenvalues  are  continuous  functions  of 
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the  matrix  elements,  it  follows  that  radical  changes  in  the  dynamical  properties  of 
a  given  system  under  differential  parameter  perturbations  ensue  principally  from 
a  large  disorientation  between  the  eigenvectors  of  the  tuned  (unperturbed)  and 
mistimed  (perturbed)  systems.  We  study,  in  Section  IV,  the  variation  of  eigenvec¬ 
tors  of  generic  famines  under  differential  random  parameter  oerturbations. 


III.  Bounds  on  Amplitude  Ratios 

Consider,  again,  the  equation  set  for  the  dynamics  of  quadranc  systems: 

V/.t’o  -*■  Kxq  =  /,  (3.1) 

where  M  and  K  are  positive  definite  matrices.  For  tuned  cyclic  systems,  M  and  K 
have  the  additional  property  of  being  circulant.  Taking  the  Laplace  transform  of 
(3.2)  under  zero  initial  conditions,  gives: 

(Ms1  +  K)Xq(s)  =  F(s),  (3.2) 


or 


A(s) -X0(r)  =  F(s)  (3.3) 

where  A  =Ms 2  +  K.  Suppressing  s  in  all  subsequent  calculations  leads  to: 

X0  =  A~lF.  (3.4) 

The  positive  definite  nature  of  M  and  K  guarantees  that  A  exists  for  all  s  on  the 
Nyquist  contour.  Under  normal  operations  of  the  system,  suppose  A  vanes  to 
A  +  AA  :=  Ae.  Let  3fo  then  change  to  Xo  +  :=  Xe.  Then,  for  the  same  excita¬ 

tion  force  as  in  the  tuned  state, 

Xt  =  (A  +  AM)'1  F.  (3.5) 

The  physical  nature  of  the  system  guarantees  that  A  +  will  always  remain 
symmetric  but  not  necessarily  circulant  since  a  true  mistumng  destroys  cyclicity. 
Equation  (3.5)  can  be  rewritten  as: 

Xt  =  (A  +  AA)_1  •  F  =  (/  +  A-]±AV]  ■  A~lF.  (3.6) 


Substituting  (3.4)  into  (3.6)  gives: 
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Xe  —  (I  +  A  *  A )  ■  Xq. 

Normally  AA  will  be  a  differential  perturbation  of  A,  so  that: 

p(A"lAA)  <  1, 

where  p(  )  is  the  spectral  radius  of  (•).  Hence 

(/  A'1  AA)_i  =  f  (-l/fA^AA/. 

<  =o 

Substituting  (3.8)  into  (  3.7)  gives: 

X'  =  Z  AA/  Xr>. 

K  =0 

Taking  norms  in  (,3.9)  gives: 

:|*«»  =  "  S  (-'/^-‘AA)*^!! 
4=0 


(3.7) 


(3.8) 


(3.9) 


s  £  n  a'!aa  r  ■  iiA’o  ii. 

4=0 


(3.10) 


Let  I A  ‘AAII  =  r.  Because  AA  is  a  differential  perturbation  of  A,  it  follows 
that  r  <  1.  Hence: 


11  £  ItXoll  £  rk  =  IIXqII 

4=0 


1  +  r  +  r 


II X0  II 

1  -r 


since  r  <  1. 


(3.11) 


Or: 


l*il  1.  _  I 
ILXqII  ~  l  -r  ~  J_  • 


(3.12) 


Write 
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A  =  D  +  C  =D(/  +  D-1C)  (3.13) 

where  ^  is  a  diagonal  matrix  of  the  uncoupled  system  dynamic  matrix  and  C  is 
the  relative  coupling  dynamic  matrix,  such  that  the  minimum  eigenvalue  of 
£>~l  C  at  any  frequency  gives  the  coupling  index  of  the  system  at  that  frequency 
[21].  If  the  norms  in  (3.12)  are  H“-  norms,  then,  over  the  frequency  interval  fi: 


ilXe((t)jil. 

ess.sup  — - - 

uieCl  IX.,  (CO)  11, 


£  ess.sup  ( 

0J6O 


qmixAA(co) 
atmnA  (co) 


where  amiuc)  and  correspond  to  maximum  and  minimum  singular  values 

of  (•)  respectively.  Note  that  all  the  matrices  and  vectors  considered  above  are 
functions  of  frequency  s-iu >. 

Because  .A  is  symmetric  it  follows  from  (3.13)  that: 

c„un^)  =  <Jmin(0[/>Z?-lC]) 

=  '  \runU  +  D  lC), 


=  (3.15) 

by  the  eigenvalue  shift  theorem,  where  [s  the  minimum  eigenvalue  of  D.  Let 
,  ,,  ,  ,  ilXe(co)  II « 


tfnun  =  a,  and  X^iD  ‘C)  =  Xq.  At  any  frequer-'v  tu,  let 
Then  (3.14)  reduces  to: 


!X0(co)  II, 


=  n(co). 


8e  =  ess.sup  n«o)  < 

wen 


,  A.A  ( CO)  I 

I  -  ess.sup  - ; - : 

wen  a(co)(l  +  A0(to.))| 


where  [2  is  a  frequency  interval  of  interest.  In  some  cases  it  is  possible  to  define 
f2  by  the  semi-open  interval  Q  =  [  0  =»).  Here  Xo  is  called  the  coupling  index  of 
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the  system.  The  system  is  decoupled  when  Xo  =  0.  It  is  weakly  coupled  if 
Xo  <  1,  and  is  strongly  coupled  if  Xo  5  1.  In  general,  0  S  Xo  2  =».  Observe  that 
Xo (k,  co)  is  a  function  of  both  the  structural  coupling  k,  and  frequency  co.  Inequal¬ 
ity  (3.16)  leads  to  the  following  conclusions: 

(i)  The  mistuned  to  tuned  amplitude  ratio  is  determined  by  the  maximum  peak 
of  the  mistunmg  strength  aTilxAA(c)),  the  minimum  strength  of  the  weakest 
link  in  the  system  auo),  and  the  minimum  peak  of  the  coupling  index 
(strength)  Xo(co). 

(ii)  A  variation  in  rigidity  (coupling)  affects  the  ratio  of  (3.16)  monotonically 
for  fixed  values  of  a~,„tAA)  and  a.  This  is  because  at  any  given  fre¬ 
quency,  Xo  vanes  continuously  and  monotonically  as  the  coupling  is  varied 
[13]. 

(iii)  A  reduction  in  a  caused  by  a  reduction  of  mass  of  the  blades,  and/or  more 
flexible  blades,  increases  the  ratio  (3.16)  monotonically.  More  specifically, 
at  any  frequency  when  Xo  ->  0,  from  (3.16): 


Se  <,  ess. sup 
<ue  Q 


_ _ 

a  (co)  -  om„AA(co) 


>  1,  for  omiXAA(co)  >  0,  VcoeQ. 


Hence  under  weak  coupling  across  the  frequency  interval,  the  amplitude 
ratio  depends  entirely  on  the  relationship  between  the  frequency  response 
of  the  mistuning  strength  and  that  of  the  strength  of  the  weakest  blade  in 
the  assembly.  Under  these  conditions,  the  maximum  amplitude  ratio  will 
arise  from  the  blade  with  the  worst  mistune  [22]. 

IV.  Eigenvector  Rotations 

In  section  II,  we  showed  that  generic  systems  Q.QV  will  typically  have  dis¬ 
tinct  eigenvalues,  while  degenerate  systems  Qv  will  typically  have  repeated 
eigenvalues.  To  study  eigenvector  perturbations  for  generic  systems,  regular 
analytical  methods  will  work,  while  for  eigenvector  variations  in  the  system  Qv 
we  require  singular  perturbations  [23].  Let  A  e  C'1"’  be  the  dynamic  matrix 
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arising  from  any  system  Q,  s  QIQV.  Let  T  represent  the  parameter  space  and  let 
pe  T  be  a  p-dimensional  parameter  vector.  If  Codim  Qv  >  r,  then  for  any 
p  e  T  e  Rp,  where  p<r,  differential  parameter  variations  in  A(p)  will  not  lead  to 
eigenvalue  degeneracies.  Thus,  if  the  eigenvalues  of  A(p),  given  by 
X2(p),  •  ■  •  ,A.a(p),  are  distinct  when  p=0  they  will  continue  to  remain  dis¬ 
tinct  when  p  is  small,  by  continuity  arguments. 

Let 


A  (5p)  =  .4(0)  +  5A, 


(4.1) 


where: 

5A  =  Sp*  •  -  A(p)  |  4=0  (4.2) 

54  can  be  expanded  in  Taylor  senes  form  as:  54  =  AA  +  higher  terms  dependent 
on  p.  To  a  first  order  approximation  we  can  write  the  perturbed  matrix  as: 

4  =  Aq  +  AA.  (4.3) 


Write 


An  =  U  A  U~ 


where  U  is  the  modal  matrix  of  Aq,  and  V'*  =  C/  1  where: 

U  -  [Ml,  U  2,  ■  •  •  .  Un] 
and 


i/*  r  •  • 

V  =  fvi  v2  ,  •  •  •  ,J 


(4.4) 


with 


A  =  diag(Xi,  X2,  •  •  • ,  X„). 

where  (•)*  is  the  complex  conjugate  transpose  of  (•)• 

Since  Aq  is  also  generic,  we  can  write  the  perturbed  modal  expression  as: 

A0  +  M=[U +AU}[A  +  Aj\)[U +  AU]~X  (4.5) 

where  A U  is  the  perturbation  in  U  resulting  from  A A  while  AA  is  the 
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corresponding  perturbation  in  A  resulting  from  A4.  Under  eigenvector  normali¬ 
zation,  II  u,  II  =  1,  and  II  u,  +  Am,  !l  =  1.  Equation  (4.5)  can  be  solved  as: 

AqU  +  A AU  +  AqAU  —  U A  +■  U AA  +  AU A  (4.6) 

where  we  neglect  second  order  terms  like  AUAA  and  AAAU  [24].  As  a  measure 
of  the  eigenvector  variations,  we  begin  by  writing  A u,  as  a  linear  combination  of 
all  the  eigenvectors  since  the  eigenvectors  span  the  whole  n-dimensional  space. 
Thus: 


Or: 


n 

Au,  =  £  Z„  u, 


j=  i 


(4.7) 


AU  =  UL.  (4.8) 

Now  solving  for  AA  in  (4.6)  gives: 

AA  =  U~l  AqU  +  U~l  A AU  +  U~x  A0  AU  -  A  -  U~l  AU  A.  (4.9) 
Observe  that  U~l  A qU  -  A  =  0,  so  that 

AA  =  U~x  AAU  +  AL  -LA.  (4.10) 

Notice  that  the  diagonal  elements  of  (AL  -  LA)  are  zero.  Hence: 

AX,  =  [  U~lAAU]u  =  v*  AAUi. 

To  solve  for  AU,  we  need  L.  The  off-diagonal  elements  of  L  are  given  by  Skelton 
[24] 

Iji  =  (X;  -  X,r‘  v*  A AUj  for  i*j,  ij  =1,2  •  •  •  n  , 


lij  =  (X,  -  X;)_l  v*  A Aui  fori*/,  i.j  =  1,  2  ••  •  n. 

To  determine  observe  that  the  constraint  equation  II  u,  +  Au,  II  =  i  contains 
Zu.  Thus 

Hu,  +  Am,  II  =  {<  ut  +  A uit  u ,  +  Au,  >  }*  =  1.  (4.12) 
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Or: 


u'  u,  +  2 u*  An,  +  A u]  A u,  -  1. 


(4.13) 


But  ut  Ui  =  1,  so  chat 


Therefore: 


2 n,  An,  -t-  An,  Aw,  =  0. 


AM;  —  ^  l  tlUj  —  2  ^jl  UJ  l ii 

;  =  i  l  =  l 


(4.14) 


=  X,  +  /„  M, 


where: 


Thus: 


Letting: 


and 


x, 


=  Z  ur 


)  «i 

i  *  t 


lu  +  (2  +  2u’x,  )/,-,  +  (2n'r,  +  jc-jc;)  =  0 


z.  =  «,*Xi 


(4.15) 


(4.16) 


y<  =x‘x,, 

gives  (on  accepting  the  positive  solution  of  the  quadratic): 

hi  =  -(1  +  2.)  +  (1  +  zf -yi)'A.  (4.17) 

Since  the  eigenvectors  u,  and  «,•  +  Am,  can  be  normalized  to  unity  and  since  each 
vector  is  represented  by  a  magnitude  m,  and  an  angle  0,,  the  natural  measure  of 
modal  variations  is  0,  since  m,sl  after  normalization.  Knowing  all  the  elements 
of  L,  we  can  now  determine  0,  as: 
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<uit  Ui  +  A Ui>  =  II  u,  II  II m,  +  i\Ui  II  cos  9 (4.18) 

But  II  Ui  II  =  II  u,  +  A ^  II  =  1. 

Hence: 


cos  0,  =  <Ui,  ut>  +  <Ui,  A ut> 


=  1  +  u'(Xi  +  ltl  u ,) 


1  f  u'x,  +  /„ 


=  (1  -f  zf  -y,/1  ,  0  S  9 ,■  S  rt/2.  (4.19) 

Consequently  for  the  occurrence  of  no  vector  rotation  under  parameter  variations, 
we  require: 

z?  -  y,  =  0  (4.20) 


Or: 


*? - x'xi  =  JC* (MjU*  - /)jt,  =  0  (4.21) 

This  implies  jc,  belongs  to  the  null  space  of  («,u*  -  /),  that  is: 

(utu'  -/)  £  Iji  Uj  =  0  (4.22) 

j  =  i 
y 

where: 


/;»  =  (A*  -  A,)'1  vj  A/4  ^  i*j 

The  nearer  the  expression  (4.20)  is  to  zero,  the  less  the  corresponding  eigenvector 
rotation  under  the  given  perturbation.  Let 


“i  =  Z  (/<  My  (w.«i*  -/)  Z  '  =  L  2  •  •  •  n. 

;  =  *  y  =  • 

;  *<  j *  i 

Then  max  (a, }  gives  the  eigenvector  with  maximum  rotation. 


The  conclusions  are  the  following: 
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a)  If  the  separation  between  the  eigenvalues  is  very  large,  (i.e.  (X,  -X;)  is 
very  large  for  all  i,j),  then  /;1  is  relatively  small  and  eigenvector  rotation 
will  be  correspondingly  small. 

b)  If  v*  AAu,  =  0,  then  eigenvector  rotation  will  also  be  relatively  small,  pro¬ 
vided  Iji  *  00  ■ 

For  example,  if  A0  is  Hermitian  as  is  the  case  in  all  quadratic  systems,  and 
AA  =  a/,  ae  C,  then 

v’  A Auj  s  0,  V  i.j  =  1,2,  •  ■  •  n. 

Thus,  identical  increases  or  decreases  in  the  diagonal  elements  of  a  quadratic  sys¬ 
tem  will  not  produce  unexpected  amplitude  excursions  [25]  because  it  cannot 
produce  eigenvalue  splittings  in  formerly  degenerate  families.  Therefore,  such 
perturbation  cannot  take  a  system  either  across  the  boundary  of  the  bifurcation  set 
or  across  different  bundles  of  Qv.  Geometrically,  this  implies  that  degenerate 
eigenvalues  in  systems  belonging  to  a  bundle  in  Qv  cannot  be  lifted  by  perturba¬ 
tions  that  leave  the  perturbed  system  in  the  same  bundle  of  Qv.  Indeed,  define 
the  eigenvector  sensitivity  matrix  of  a  quadratic  system  as 

5  =  AULT1  =  L, 

from  eqn.  (4.8).  Defining  the  eigenvector  sensitivity  metric  measure  by 

Sp  =  IIS  II  J  =  £  l\ 

!  =  1 
i  =  1 

where  Sp  is  the  Frobenius  norm  of  S  shows  that  the  maximum  eigenvector  sensi¬ 
tivity  is  obtained  at  the  positions  of  minimum  eigenvalue  separation,  which  is  not 
difficult  to  compute.  Alternatively,  (Sp)mtx  also  corresponds  to  the  position  of 
maximum  angular  rotation  between  the  tuned  and  mistuned  system  eigenvectors. 
This  condition  is  evidenced  by  strong  eigenloci  deformations. 

If  A(a>)  is  a  frequency  response  matrix  arising  from  a  generic  system,  the 
eigenvalues  Xi(CD)  and  eigenvectors  M;(o>)  are  also  continuous  functions  of  fre¬ 
quency.  We  can  therefore  plot  the  frequency  response  functions  Sp (co)  to  deter¬ 
mine  the  fre  -'iencies  at  which  maximum  deformations  take  place. 
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V.  Examples 

To  illustrate  the  theory  so  far  developed  we  consider  two  examples.  The  first 
is  an  interconnected  linear  chain  of  oscillators.  This  has  been  studied  by  Arnold 
[13]  and  more  recendy  by  Pierre  [181. 

Example  1;  Mode  Localization  in  Generic  Periodic  Systems. 

Considt:.  a  coupled  pendulum,  as  shown  in  Fig  1,  with  identical  masses  but  of 
different  lengths  l\  and  Z2,  where  /;  is  a  perturbation  of  l\,  i.e.,  Z2  =  (i\  +  A/i).  If 
we  put  /i  =  /,  then  the  kinetic  energy  is  given  by 

T  =  7  m[/20f  +  (/  +  Si)2  Q\]  (5.1) 

while  the  potential  energy  is  given  by 

d2  q2 

U  =  ml-~  +■  mil  +  &l)~  +  “ (0i  -  ©2)2-  (5-2) 


Fig.  I.  Two  coupled  oscillators. 


Under  unit  gravitational  force,  application  of  Lagrange’s  equations  results  in  the 
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equation  of  motion: 


M 


K 


=  0 


where 


M  = 


ml 2  0 

0  m  {l  +  Ai)2 


and 


K  = 


ml  +  k  -k 

-k  m{l  +  Ai)2  +  k 


The  dynamic  matrix  for  tne  above  system  is  given  by 
A  (to)  = 


[ml  +  k  -  co 2ml2  -k 

-k  m{l  +  &l)  +  k-di^il  +  Alj2 


Rewrite  A  (co)  as: 


A  (CO) 


*  fr  ?]• 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


The  characterisdc  equation  of  A  (co)  is  given  by 

X2  -  (a  +b)X  +  {ab  -  k2)  =  0 

Both  M  and  K  are  symmetric  and  positive  definite.  The  eigenvalues  of  A  (  )  are: 


*1.2  -  , 


(5.8) 


Note  that  A.],2  cannot  be  degenerate.  Thus  under  one-parameter  deformations,  the 
eigenvalues  can  deform  but  cannot  collide. 

Indeed, 


3^ 

da 


1 

2 


a  -  b 

3 ka-b )2  +  4 k2  . 


(5.9) 


and 
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Hence: 


dh  _  1 

da  2 


dX[  dX2 
da  da 


1  - 


a-b 


V(a  -  b)2  +  4 k2 


1 

2 


when  a=b. 


(5.10) 


(5.11) 


The  distance  between  the  eigenvalues  is  given  by: 

=|*i  -X2|=  'i(a-b)2  +  4*2  (5.12) 


which  assumes  its  minimum  value  of  2k  when  a-b  or  when 


dXi/da 


=  Vi  This 


represents  the  tuned  state  of  the  linear  chain.  For  a  fixed  mistiming  value  (a-b), 
d\  depends  essentially  on  k.  If  (a-b)  is  small,  it  is  clear  that  Sp  -> »  as  k  -*  0. 
The  modal  matrix  of  the  chain  is  given  by 


i/  = 


-i 


-1 


(a-b)  -  V(a  -  b )2  +■  41:2  (a-b)  +  -b)2  +  4  k2 

2k  2k 


(5.13) 


O';  serve  that  u‘  Uj  =0,  V  k,a,b.  Under  tuned  conditions,  a=b,  then 

u.  ■  [:!  i 


(5.14) 


However,  consider  the  very  interesting  situation  when  the  mismning  to  coupling 
ratio  is  rather  large.  That  is  to  say: 

k 

Then  (a  -  b)2  3*  k2,  and  k2  becomes  negligible  in  the  eigenvector  expressions. 
Expanding  the  term  under  the  radicals  and  neglecting  second  and  higher  order 
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terms  gives: 

[{a  -b)2  +  k2)Vl  =  (a  -b)  1  + 

!  k  l2 

=  (a  - b)[l  +  -r  - -  +  (5.15) 

2  a  -  b 

In  this  case  the  modal  matrix  reduces  to: 

-1  -1 

k  a  -  b 
a  -  b  k 

An  energy  exchange  now  takes  place.  The  second  component  of  the  first  mode 
becomes  vanishingly  small  while  the  corresponding  component  of  the  second 
mode  becomes  extremely  large.  This  is  an  extreme  case  of  classical  vibration 
absorber,  and  is  the  mode  localization  phenomenon.  We  therefore  conclude  that 
mode  localization  (or  extreme  energy  exchange)  will  occur  in  a  generic  system 
under  one-parameter  deformations  if  the  following  conditions  are  satisfied: 

•  at  any  frequency  co  where  the  system  is  almost  decoupled,  i.e.,  X(co)0  -*  0. 
(Note  that  A.(cu)o  ->  0  as  *->0). 

(1  ~ 

•  when  the  mistuning  to  coupling  ratio  — - —  »•  1. 


(5.16) 


At  the  localization  stage  the  eigenvalue  and  eigenvector  sensitivities  take  on  their 
maximum  values,  i.e.  both  II  AAA'1  II  \  and  II  At/  U~l  Ilf  have  their  maximum 
values.  Localized  modes  always  produce: 


5e  = 


\\xe  II «, 
lijCQ  H- 


(5.17) 


Example  2:  Cyclic  Systems. 

Consider  three  identical  masses,  m,  arranged  in  a  ring  structure  and  intercon¬ 
nected  by  identical  springs  kc.  Assume  that  all  the  masses  are  hinged  to  the 
ground  by  torsional  springs  of  strength  ktl  and  that  the  radius  of  the  ring  is  r;  as 
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shown  in  Fig.  2. 

The  basic  equations  of  motion  of  this  “ring”  is 

Mx  +  Kx  =  f  (5.18) 

where 


k, 

2kc  +  — 

r 

~kc 

-*c 

m 

0 

01  1 

I, 

0 

m 

ii 

5* 

O 

-kc 

2kc  *  ~ 

—ic 

0 

0 

m\ 

r 

c 

-*C 

-  kc 

k, 

2  kc  +  — 

r 

Fig.  2.  Models  of  (a)  the  cyclic  chain,  (b)  the  linear  chain 
with  three  degress  of  freedom. 
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Using  group  theoretic  arguments  [13],  we  can  easily  deduce  that  the  above  sys¬ 
tem  has  degenerate  eigenvalues  occurring  as  doublets,  by  cyclicity  of  the 
corresponding  system  matrices.  Consequently,  every  quadratic  cyclic  system 
Qc  cfiv  Furthermore  all  perturbations  of  the  above  system  preserving  the 
cyclic  structure,  leaves  the  modal  geometry  invariant  [3,  25].  Indeed  the  eigen¬ 
values  of  the  above  system  are  given  as: 


X-—  X  -  —  +  —  X)  -  —  +  — 

1  mr '  mr  m  mr  m 


(5.20) 


Write  the  dynamic  matrix  of  this  system  as: 


Ao(cu)  = 


’  a 

-b 

rb 


-b 

a 

-b 


-b 

-b 

a 


where 


a  =  2kc  +  k,  -  0)2m,  b  =  kc. 


We  now  consider  a  diagonal  perturbation  of  the  form: 

E  =  diag  («i,  ei,  e3). 


Then 


A*(o))  =  A0(co)  +  E. 


(5.21) 


(5.22) 


(5.23) 


This  would  correspond  to  the  realistic  situation  where  there  are  slight  changes  in 
the  values  of  the  ground  spring  kt,  depending  for  example  on  how  the  blades  are 
coupled  to  the  disk  in  bladed  disk  assemblies  [22].  The  major  difference  between 
the  behavior  of  (degenerate)  cyclic  systems  and  generic  systems  are  the  follow¬ 
ing: 

(i)  For  generic  systems,  all  the  eigenvalues  and  the  distance  between  adjacent 
pairs  increases  as  the  coupling  kc  increases.  Consequently  the  probability 
of  mode  localization  decreases  as  kc  and  hence  X(d))0  increases.  On  the 
other  hand,  perturbations  which  split  the  degenerate  eigenvalue  of  cyclic 
systems  turn  them  into  generic  systems  with  pathologically  close  eigen¬ 
values  [7].  Hence  for  previously  cyclic  systems  whose  eigenvalues  bifur¬ 
cate  under  perturbations,  Sp  is  very  large.  Therefore  such  systems  are 
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susceptible  to  mode  localization,  independent  of  the  values  of  the  coupling 
strength  kc .  Recall  that 

SF  =  II L  II 

where 

ltJ  -  (Xj  ~  A-,  J_l  v*  A AUj,  i*j,  ij  =1,2,  •  ■  ■  n. 

(ii)  Consequently  the  only  way  to  avoid  large  values  of  SF  in  such  a  situation  is 
if  and  only  if  :l  v*  A Aw,  II  s  0  or  in  the  neighborhood  of  zero.  Perturbations 
that  induce  this  condition  are  precisely  those  that  will  not  induce  radical 
dynamical  changes  in  mistuned  cyclic  systems.  It  was  already  shown  that 
if  AA  =  a/  ,  then  ;i  v*  AAu,  II  =  0 

(iii)  Of  the  remaining  possible  perturbations  those  that  have 
II  vf  A Au,  II  =  e  <<  1  will  produce  minimum  dynamical  changes.  All  oth¬ 
ers  for  which  ;l  v*  A Aw,  II  is  not  small  will  give  susceptibility  to  mode  local¬ 
ization,  no  matter  how  strong  the  interblade  coupling. 

The  following  numerical  example  amplifies  the  above  observations.  We  con¬ 
sider  the  case  of  the  so-called  ‘strong  coupling’,  using  the  following  values: 
kc=  9.5,  k,=  1,  a  =  20,  b  =  9.5,  =0,  ei =-0.1,  ex  =0.1.  Clearly,  the  ratio  of  mis- 

tuning  to  coupling  strength  is  very  small.  Now,  in  order  to  compute  the  frequency 
response  curves,  we  need  some  damping  to  obtain  finite  amplitudes  at  resonance. 
Assume  hysteretic  damping  of  0.01  for  all  cases.  Without  loss  of  generality,  the 
response  to  be  computed  is  uie  direct  recept ance,  i.e.  the  response  of  each  node  to 
individual  excitation.  We  turn  the  ring  into  a  linear  chain  by  putting 
b  =  k i3  =£31  =0  in  equation  (5.21).  Then  Ag  becomes  a  tridiagonal  banded 
matrix. 

The  frequency  response  of  the  tuned  and  mistuned  systems  of  the  linear  chain 
are  shown  in  Fig  3.  The  illustration  is  windowed  around  one  of  the  resonant  fre¬ 
quencies  of  the  coupled  system.  Nonce  that,  at  the  tuned  state,  ,.ie  amplitudes  of 
nodes  1  and  3  are  equal  on  account  of  symmetry,  while  that  of  node  2  is  double 
that  magnitude. 

Because  the  system  is  now  generic,  and  therefore  exhibits  modal  stability,  all 
nodes  have  almost  the  same  response  patterns  and  magnitudes  as  in  the  tuned 
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system.  This  is  also  the  case  when  we  change  the  sign  of  e2,  from  -0.1  to  0.1. 

When  wt  repeat  exactly  the  same  procedure  for  the  circulant  system,  a  very 
dl'Jerent  picture  is  obtained.  Fig  4  shows  the  response  of  individual  nodes  com¬ 
pared  with  the  tuned  case.  This  case  corresponds  to  a  2-parameter  perturbation, 
with  e  i  =  0. 1 ,  e  2  =  -0.1,  =  0. 


(. a )  node  ' 


(0)  node  2 


(c )  node  3 


Fig,  3.  Effect  of  mistuning  on  th „  response  curves  of  the  linear  chain.  Note  the 
preservation  of  the  shape  of  the  curves  around  resonance,  and  the  minimal 
difference  in  the  peak  amplitudes  of  the  tuned  and  mistuned  systems  -  -  -  -  tuned 
systems;  _ mistuned  system). 
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Ul  u*  III 

frequency 


Fig.  4.  Effect  of  two  parameter  mistiming  on  the  response  curve  of  the  cyclic 
chain.  Note  the  severe  reduction  in  the  amplitude  at  node  3,  which  is  only  50%  of 
the  tuned  system  (  -  -  -  -  tuned  systems;  _ mis  tuned  system). 


Notice  that  the  node  with  zero  mistuning  (mode  3)  now  has  a  reduction  in 
amplitude  of  almost  50%.  This  extremely  unequal  amplitude  distortion  (Fig  4)  is 
the  case  no  matter  how  small  the  magnitude  of  the  perturbation  is,  so  long  as  we 
keep  the  form  of  mistuning,  and  the  mistuning  does  not  actually  vanish. 

If  we  now  change  the  mistuning  matrix  in  a  very  small  way,  by  making 
e 2=0.1,  we  obtain  the  response  curves  in  Fig  5.  We  now  notice  a  substantial 
difference  in  the  geometry  of  the  curves  in  Fig  5,  compared  to  those  in  Fig  4. 
Thus,  a  very  small  change  in  the  perturbation  matrix,  now  results  in  a  consider¬ 
able  difference  in  the  vibration  response  at  the  individual  nodes.  The  question  of 
which  node  will  be  most  responding,  or  the  one  having  the  least  amplitude,  is  now 
not  as  easy  as  one  would  have  expected.  In  Fig  4,  it  is  node  3,  while  it  is  node  2  in 
Fig  5.  In  fact,  the  amplitude  of  node  3  has  been  increased  by  about  100%  from 
Fig  4  to  Fig  5,  merely  by  changing  only  one  entry  in  the  system  matrix  from  19.9 
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to  20.1,  a  change  of  less  than  \%  ! 

The  foregoing  examples,  based  on  a  simple  3  degrees  of  freedom  model  of  a 
circular  ring  or  disk  only,  illustrates  the  instability  induced  by  cyclicity.  It  is 
clear  that  the  qualitative  conclusions  to  be  drawn  from  Fig  4  are  inconsistent  with 
those  from  Fig  5,  although  the  difference  between  the  two  mistuned  matrices  is 
very  small  indeed.  We  emphasize  that  these  results,  obtained  for  just  a  cyclic 
chain,  are  not  necessarily  applicable  to  bladed  disks  in  all  generality,  especially 
chose  models  in  which  cyclicicy  is  ignored.  However,  when  bladed  disk  systems 
are  well-modeled  to  include  the  effects  of  blade  coupling,  blade  or  isk  mistuning 
and  cyclicity,  similar  distortions  in  the  geometry  of  the  frequency  response  curves 
can  result.  The  subject  is  currently  under  investigation  by  us. 


(a )  nod*  1 


(&)  node 2  node 3 


Fig.  5.  Effect  of  one-parameter  mistuning  on  the  response  curve  of  the  cyclic 

chain.  Note  the  symmetrical  unfolding  of  the  degenerate  singularity  ( - tuned 

systems;  _ mistuned  system). 
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VI.  Conclusions 

(i)  For  generic  systems,  to  which  linear  periodic  chains  of  oscillators  belong, 
differential  parameter  perturbations  are  significant  for  the  system  dynamics 
only  under  weak  coupling  conditions  when  the  mistuning  to  coupling  ratio 
exceeds  unity  (Example  1).  Under  all  other  conditions  that  do  not  induce 
eigenvalue  degeneracy;  small  magnitudes  of  mistuning,  or  the  type  of  mis¬ 
tuning,  is  irrelevant  to  system  dynamics. 

(ii)  For  degenerate  systems  to  which  a  tuned  cyclic  system  with  circulant 
dynamic  matrices  belongs,  it  is  not  just  the  mistuning  to  coupling  ratio 
which  is  significant  in  the  determination  of  the  perturbed  system  dynamics. 
The  type  of  mistuning  assumes  a  far  greater  importance  than  the  mistuning 
to  coupling  ratio.  All  types  of  mistuning  that  move  the  system  either  across 
the  boundary  of  the  bifurcation  set,  or  from  one  fiber  bundle  of  the  degen¬ 
erate  set  to  another  within  Qv  will  lead  to  topological  catastrophes  [15]. 
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APPENDIX  4 


The  Stability  of  Frequency  Response  Curves 


SUMMARY 


In  this  Appendix,  we  highlight  one  of  the  results  obtained  so  far,  namely:  that  the 
modes  of  vibration  of  cyclic  structures  are  unstable  under  arbitrarily  small 
perturbation.  This  instability  is  not  the  usual  ill-conditioned  problem  of 
numerical  analysis.  It  has  nothing  to  do  with  the  computational  algorithm.  The 
eigenvector  instability  results  because  the  perfect  cyclic  system  is  "physically  ill- 
conditioned",  since  a  very  small  perturbation  changes  its  dynamics 
characteristics  dramatically.  This  is  significant  because  many  aerospace 
structures  have  circular  profiles.  The  implication  of  eigenvector  instability  for 
modal  control,  forced  response  amplitudes,  sensitivity  analysis,  etc,  therefore 
needs  further  investigation. 


Eigenvector  Stability,  Forced  Response,  and  Turbine  Blade  Failure 

The  structural  integrity  of  turbine  blades  used  in  jet  propulsion  systems  is 
sometimes  compromised  by  the  rare,  but  very  dangerous,  failure  of  some  "rogue 
blades".  This  problem  has  been  addressed  by  different  investigators  of  the 
mistiming  problem.  However,  they  often  obtained  conflicting  results. 

This  is  because  the  response  obtained  from  each  rotor  studied  by  each  author 
depends  on  the  eigenvectors  of  the  rotor  system  matrix  A.  In  general,  the  matrix 
A  will  be  different  for  each  model  used  by  each  author,  although  the  difference 
may  be  very  small.  In  fact,  mistiming  is  usually  small. 

However,  the  problem  created  by  mistiming  is  not  always  small.  Thus,  although 
each  A  in  the  family  is  differentiably  dependent  on  the  mistiming  parameter  e  in 
the  neighborhood  of  the  origin  of  £,  the  corresponding  eigenvectors  is  not. 


-2- 


Consequently  two  almost  similar  rotors  may  produce  dramatically  different 
vibration  responses,  if  their  respective  system  matrices  are  different 
perturbations  of  the  same  nominal  matrix.  An  effective  demonstration  of 
unstable  frequency  response  curves  in  a  simple  3  degree  of  freedom  cyclic 
system  is  given  in  the  following  examples.  First,  we  examine  the  instability  of 
eigenvectors,  then  the  instability  of  frequency  response  curves. 


Numerical  Examples  Illustrating  Eigenvector  Instability 


Example  1 

At  least  three  coordinates  are  required  to  define  a  cyclic  system  uniquely. 
Therefore,  we  consider  the  simplest  possible  example:  a  3x3  circulant  matrix 
with  real  elements,  a,  b  €  R. 


a  0  = 


a  b  b' 
b  a  b 
b  b  a 


(A4.1) 


Using  the  following  perturbation  matrices,  where  £  e  R  is  a  very  small  parame¬ 
ter,  we  can  generate  two  matrices  A\  =a0  +  E\  and  A2  =  a0  +  E2  that  are  very 
close,  and  such  that  these  depend  smoothly  on  e,  and  as  e  0,  A  \  ->  a  0  <-  A  2- 
Thus,  if 


then 


and 


’  £ 

0  O' 

'  £ 

0 

O' 

0 

e  0 

e2  = 

0 

-£ 

0 

.0 

0  0, 

0 

0 

0 

(A4.2,A4.2) 


A\  - 


'  a  +  e 
b 
b 


b  b 
a  +  £  b 
b  a 


(A4.4) 


A  2 


a  +  £  b  b 
b  a-£  b 
b  b  a 


(A4.5) 
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Note  that  I  I  A\  I  I  =  I  I  A2  I  I ,  where  !  !  •  II  is  some  norm. 

Consider,  for  an  illustration,  a  situation  when  a  =  100,  b  =  45,  e  =  0.1.  We  can  com¬ 
pute  the  eigenvalues  of  A j  and  A  2  respectively  as 

A!  =  diag  (10.0667,  145.0333,  145.1000)  (A4.6) 


and 


A2  =  diag  (10.0000,  144.9423,  145.0578)  (A4.7) 


Notice  that  the  eigenvalues  of  the  two  matrices  are  very  close.  If  we  now  com¬ 
pute  the  corresponding  eigenvectors,  we  get 


U 1  = 


'  .9993  -.5004 
.9993  -.5004 
1.0000  1.0000 


1.0000' 

-1.0000 

0 


(A4.8) 


U2  = 


.9985 

1.0000 

.9993 


-  .2681 
1.0000 
-  .7328 


1.0000' 

-  .2678 

-  .7313 


(A4.9) 


We  now  notice  a  significant  difference  between  the  eigenvectors  at  modes  2  and  3 
of  matrices  A\  and  A2  respectively.  For  example,  there  is  no  node  (a  point  where 
displacement  is  zero)  in  the  third  mode  of  U2/  whereas  there  exists  such  a  node 
mUv 


Example  2 


In  the  second  example,  we  consider  the  following  circulant;  its  elements  are  com¬ 
plex  but  the  matrix  is  not  symmetric.  It  may  be  regarded  as  a  deformation  of  a 
symmetric  circulant. 


^1 


‘  200  +  i(-10) 
-95  + 1  (5) 
-95  +  i(  —  5) 


-95  +  t’(  —  5) 
200  +  i(  - 10) 
-95  +  i(5) 


-95  +  1(5)  ' 
-9  5  + 1(  -  5) 
200  +  /'( -  10) 


(A4.10) 


We  test  for  modal  stability  by  computing  the  eigenvalues 


and  the  eigenvectors 


'1  +  i(0)  -  .5  +  it 866)  -  .5  +  i(  -  .866)' 

Ji,=  1  + 1(0)  1  +  i(0)  1  +  i(0)  (A4.12) 

1  +  i(0)  -  .5  +  i(  -  .866)  -  .5  +  i(.866) 


Now,  we  apply  a  very  small  perturbation  to  the  matrix  At  to  get: 


'201  +  i(  - 10)  -  95  +  i(  -  5)  -  95  +  i(5)  ' 

A2  =  -95  +  i(5)  199  +  i(  - 10)  -95  +  i(-5)  (A4.13) 

[-95  +  /  ( —  5)  -  95  +  i(5)  200  +  i(  - 10) 


It  is  clear  that  the  matrices  A\  and  A2  are  'close',  since  I  I  £  I  I  =0,  where 


'  - 1  +  i(0)  0  +  i(0)  0  +  i(0) " 

E  =  Ai  -  A2  ~  0  +  i(0)  l  +  i(0)  0  +  i(0)  (A4.14) 

0  +  i(0)  0  + 1(0)  0  +  i(0) 


The  computed  eigenvalues  of  A  2  are: 


A2  =  diag  9.9977  +  i(  - 10),  286.3218  +  i(  - 10),  303.6807  +  i(  - 10)  KA4.15) 


Now,  notice  what  happens  to  the  third  eigenvector  of  A  i  (eq.  (A4.12)),  as  a  very 
small  change  is  made  using  £,  eq.  (A4.13),  to  transform  it  to  A  2.  The  eigenvector 
matrix  of  A  2  is: 


’  0.993  +  i(O.OOO)  -  .474  +  i(0.820)  -  0.556  +  i(  -  .867)' 

U2  =  1.000  +  1(0.000)  1.000  +  i( -0.00)  -  .444  +  i(-.861)  (A4.16) 

0.997  +  i(-0.000)  -  .531  +  i (0.817)  1.000  +  i(0.000) 


Again,  it  should  be  noted  that  a  very  small  change  in  the  matrix  A  ^  induces  a 
significant  qualitative  difference  in  the  eigenvector  at  certain  modes  of  A]  (eq. 
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A4.10)  compared  with  the  corresponding  eigenvectors  of  A  2,  (eq.  A4.16). 

From  Modal  Analysis,  sometimes  known  as  eigenfunction  expansion,  we  know 
that  the  forced  response  amplitudes  are  related  to  eigenvectors.  Thus,  if  the 
eigenvectors  are  unstable  under  arbitrary  perturbation,  then,  the  forced  response 
curves  will  also  be  unstable  under  arbitrary  perturbation.  This  is  illustrated  in 
Figs  A4.1  to  A4.2  below. 


node  2  node  3 


Fin  44-1  Effect  of  one-paramctcr  mistuning  on  the  response  curve  of  ilic  cyclic  mem- 
braneT" Note  the  symmetrical  unfolding  of  the  degenerate  singularity. 

( - tuned  system;  _ mtstuned  system) 


FigA4-2  Effect  of  two-parameter  mistiming  on  the  response  curve  of  the  cyclic  mem¬ 
brane.  Note  the  severe  reduction  in  the  amplitude  at  node  3,  which  is  only  50%  of  the 
tuned  system 


( - tuned  system; 


mistuned  system) 
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IEEE  Mexicon  1976,  Conference.  Mexico  City  1976,  Paper  1.9. 

3.  Nwokah,  O.D.I.,  Progress  in  the  design  of  multivariable  control  systems,  Proc.  4th 
IFAC  Symposium  on  Computer  aided  design  of  multivariable  technological  system, 
Fredericton,  Canada,  1977,  Paper  3.1. 

4.  Nwokah,  O.D.I.,  Hadamard  transfer  matrices  and  diagonal  dominance.  Proc.  IEEE 
Mexicon  1977  Conference.  Mexico  City  1977,  Paper  4.7. 

5.  Nwokah,  O.D.I.,  Synthesis  of  feedback  systems  for  specified  time  domain 
insensitivity  to  interaction  induced  plant  ignorance.  Proc.  6th  IFAC  Symposium  on 
computer  aided  design  of  multivariable  technological  systems,  Purdue,  Indiana, 
U.S.A.  1982. 

6.  Nwokah,  O.D.I.^  The  robust  decentralized  stabilization  of  composite  feedback 
systems,  Proc.  4  IFAC  Symposium  on  Large  Scale  Systems,  Zurich,  Switzerland, 
1986. 

7.  Nwokah,  O.D.I.,  On  decentralized  control  of  complex  feedback  systems, 
Proceedings  American  Control  Conference,  Seattle,  WA,  Vol.  1,  pp.  326-331, 1986. 

8.  Nwokah.,O.D.I.,  The  quantitative  design  of  robust  multivariable  control  systems, 
Proc.  25m  IEEE  Conference  on  Decision  and  Control,  Athens,  Greece,  pp.  16-24, 
1986. 

9.  Nwokah,  O.D.I.,  Pseudo-derivative  feedback  control.  Proceedings  American 
Control  Conference,  Minneapolis,  MN,  pp.  1811-1814,  1987. 

10.  Nwokah,  O.D.I.,  Strong  Robustness  in  Uncertain  Multivariable  Systems,  Proc.  27th 
IEEE  Conference  on  Decision  and  Control,  Austin,  Texas,  December  1988. 

11.  Nwokah,  O.D.I.,  Strong  diagonal  dominance  and  closed  loop  performance,  Proc. 
ASME  Winter  Annual  Meeting,  Chicago,  Illinois,  December  1988. 

12.  Nwokah,  O.D.I. ,  Strong  robustness  in  uncertain  multivariable  systems,  Proc.  12th 
IMACS  World  Congress  on  Scientific  Computation,  Paris,  1988. 
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13.  Nwokah,  O.D.I.,  Thompson,  D.F.,  Algebraic  and  topological  aspects  of  quantitative 
feedback  theory,  Proc.  American  Control  Conferences,  Pittsburgh,  PA,  June  21-23, 
1989. 

14.  Afolabi,  D.,  Nwokah,  O.D.I.,  Fci  lurbations  of  periodic  systems,  Proc.  21st  Midwest 
Mechanics  Conference,  Michigan  Tech.  University,  July  1989. 

15.  Afolabi,  D.,  Nwokah,  O.D.I.,  The  frequency  response  of  mistuned  cyclic  systems, 
Proc.  ASME  Vibration  Conference,  Montreal,  PQ,  Canada,  September  17-19, 1989. 

16.  Burhoe,  J.C.A.,  Nwokah,  O.D.I.,  Multivariable  control  of  a  biaxial  machine  tool, 
ASME  Winter  Annual  Meeting,  San  Francisco,  December  1989. 

17.  Perez,  R.A.,  Nwokah,  O.D.I.,  The  integrated  control  of  a  propulsion-air  frame 
system,  ASME  Winter  Annual  Meeting,  San  Francisco,  December  1989. 

18.  Thompson,  D.F.,  Nwokah,  O.D.I.,  Stability  and  optimal  design  in  quantitative 
feedback  theory,  ASME  Winter  Annual  Meeting,  San  Francisco,  December  1989. 

19.  Nwokah,  O.D.I.,  Thompson,  D.F.,  Robust  low  order  controller  design  in  quantitative 
feedback  theory,  Proc.  28  CDC,  Tampa,  FL,  December  13-15, 1989. 

20.  Perez,  R.  A.,  Nwokah,  O.D.I.,  Multivariable  control  of  an  integrated  propulsion- 
airframe  system,  3rd  USC  Workshop  on  Control  Mechanics,  January  1990. 

21.  Happawana,  G.S.,  Nwokah,  O.D.I.,  Bajaj,  A.K.,  Afolabi,  D.,  On  the  dynamics  of 
perturbed  cyclic  symmetric  systems,  3rd  USC  Workshop  on  Control  Mechanics, 
January  1990. 

22.  Thompson,  D.F.,  Nwokah,  O.D.I.,  Frequency  response  specifications  and  sensitivity 
functions  in  quantitative  feedback  theory,  Proc.  ACC,  San  Diego,  599-604, 1990. 

23.  Thompson,  D.F.,  Nwokah,  O.D.I.,  Optimal  loop  synthesis  in  quantitative  feedback 
theory,  Proc.  ACC,  San  Diego,  626-631, 1990. 

24.  Afolabi,  D.,  Nwokah,  O.D.I.  Effects  of  mild  perturbations  on  the  dynamics  of 
structures  with  circulant  matrices,  Proc.  1990  AIAA/ASME/ASCE/AHS,  Structures, 
Structural  Dynamics  and  Materials  Conference,  Long  Beach,  CA,  1990. 

25.  Happawana,  G.S.,  Afolabi,  D.,  Nwokah,  O.D.I. ,  Bajaj,  A.K.,  On  the  dynamics  of 
perturbed  symmetric  systems,  Proc.  11th  U.S.  National  Congress  of  Applied 
Mechanics,  Tucson,  AZ,  May  21-25, 1990. 

26.  Nwokah,  O.D.I.,  Perez,  R.A.,  On  multivariable  stability  in  the  gain  spaced,  Proc. 
EEEE  Conference  on  Decision  and  Control,  Honolulu,  Hawaii,  December  1990. 

27.  Nwokah,  O.D.I.,  Thompson,  D.F.,  Perez,  R.A.,  On  the  existence  of  QFT  controllers, 
ASME-WAM,  Dallas,  TX,  1990. 

28.  Nwokah,  O.D.I.,  Jayasuriya,  S„  Chait  Y.,  Constrained  robust  disturbance 
accommodation  by  quantitative  feedback  theory  (QFT),  4th  USC  Workshop  on 
Control  Mechanics,  University  of  Southern  California,  January  22-24,  1991. 
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D.  Technical  Reports: 

1.  Nwokah,  O.D.I.,  On  design  of  multivariable  control  systems.  Control  Systems 
Centre  report  No.  429,  UMIST,  Manchester,  England,  September  1978. 

2.  Nwokah,  O.D.I.,  Weak  diagonal  dominance  and  the  Nyquist  array,  Research  report 
IAS  (1980),  Manchester  Polytechnic,  Manchester,  England. 

3.  Nwokah,  O.D.I.,  Strongly  restricted  interaction  and  minimum  variance  controllers 
for  multivariable  systems.  Research  report  IAS,  (1980)  Manchester  Polytechnic, 
Manchester,  England. 

4.  Nwokah,  O.D.I.,  The  assignment  of  loop  interaction  in  linear  multivariable  systems. 
Research  report  IAS  (1980),  Manchester,  England. 

5.  Nwokah,  O.D.I.,  Multivariable  describing  functions  and  the  Nyquist  array,  Research 
report  IAS  (1979),  Manchester  Polytechnic,  Manchester,  England. 

6.  Nwokah,  O.D.I.,  The  integrity  and  robustness  of  systems  with  integral  control. 
Control  Systems  Technical  Note  #30,  Dept  of  Electrical  Engr.,  University  of  Notre 
Dame,  Oct.  1984. 

7.  Nwokah,  O.D.I.,  The  robust  decentralized  stabilization  of  composite  dynamical 
systems.  Control  Systems  Technical  Note  #31,  Dept,  of  Electrical  Engr.,  University 
of  Notre  Dame,  Oct.  1984. 

8.  Nwokah,  O.D.I.,  The  robust  multi-level  hierarchical  control  of  composite  feedback 
systems.  Control  Systems  Technical  Note  #32,  Dept  of  Electrical  Engr.,  University 
of  Notre  Dame,  Oct.  1984. 

9.  Nwokah,  O.D.I.,  On  non-singular  value  based  design  of  controllers  for  robust 
stability.  Control  Systems,  Technical  Note  #33,  Dept  of  Electrical  Engr.,  University 
of  Notre  Dame,  Dec.  1984. 

10.  Nwokah,  O.D.I.,  Sensitivity  and  robustness  optimization  of  linear  multivariable 
systems.  Control  Systems  Technical  Note  #35,  Dept,  of  Electrical  Engr.,  University 
of  Notre  Dame,  March  1985. 

E.  Some  Invited  Presentations. 

1.  The  Design  of  Linear  Multivariable  Systems  -  University  of  London  (Imperial 
College),  Dept,  of  Electrical  Engineering,  June  18, 1975. 

2.  Rapid  Changes  of  State  in  the  Nigerian  Socio-economic  System:  Applications  of 
Singularity  Theory,  University  Engineering  Lecture  Series,  University  of  Nigeria, 
March  1976. 

3.  Progress  in  the  Design  of  Multivariable  Systems,  Control  Systems  Center, 
University  of  Manchester  Institute  of  Science  and  Technology,  July  1977. 

4.  Applications  of  Catastrophe  Theory  to  Modelling  of  the  Nigerian  Socio-economic 
System.  Nigerian  Institute  for  Social  and  Economic  Research,  University  of  Ibadan, 
April  1977. 

5.  Interaction  Indices  for  Multivariable  Systems,  Joint  Control  Seminar  Series, 
Sheffield  University/Sheffield  Polytechnic,  Dept,  of  Electrical  Engineering, 
Sheffield,  England,  November  1979. 
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6.  Multivariable  Control  of  Aircraft  Fatigue  Testing  Rigs,  Department  of  Electrical 
and  Electronic  Engineering,  University  of  Surrey,  England,  March  1980. 

7.  A  Control  Theory  Component  to  the  German  Academic  Exchange  Program 
(DAAD)  on  Industrial  Mathematics,  Department  of  Mathematics,  University  of 
Ibadan,  Nigeria,  February  1981. 

8.  Frequency  response  methods  for  the  control  of  uncertainty  in  multivariable  plants. 
Case  Western  Reserve  University,  Department  of  Electrical  Engineering  and 
Applied  Physics,  Cleveland,  Ohio,  April  1985. 

9.  Polynomial  based  feedback  design  and  Bode’s  frequency  response  conditions. 
Department  of  Electrical  Engineering,  Seattle  University,  Seattle,  Washington, 
March  1985. 

10.  Are  Engineering  Control  and  H“  Optimization  Compatible,  Department  of 
Electrical  Engineering,  University  of  Notre  Dame,  Notre  Dame,  Indiana,  February 
1985. 

11.  M-matrices  and  Quantitative  Robustness  Measures  in  Multivariable  Uncertain 
Systems,  Harris  Corporation,  Government  Aerospace  Systems  Division,  Melbourne, 
Florida,  August  1987. 

12.  Practical  Aspects  of  Sensitivity  Reduction  in  Engineering  Control  Systems,  Purdue 
University  at  Indianapolis,  School  of  Engineering  and  Technology,  November  1987. 

13.  The  Dynamics  of  Mistuned  Cyclic  Systems,  Electric  Power  Research  Institute,  Palo 
Alto,  CA,  October  13, 1988. 

14.  Integrated  Dynamics  and  Control  of  Propulsion-Airframe  System,  Duke  University, 
Durham,  NC,  Dept,  of  Mechanical  Engineering  and  Materials  Science,  April  5, 
1989. 

15.  Application  of  QFT  to  Engine- Airframe  Integration,  University  of  Illnois  at 
Urbana-Champaign,  Dept,  of  Mechanical  and  Industrial  Engineering,  April  18, 
1989. 

16.  Model  Reference  QFT  Application  to  Ship  Propulsion  Systems,  Texas  A&M 
Univesity,  Dept,  of  Mechanical  Engineering,  March  22,  1991. 

F.  Thesis  Supervision.  (Major  Advisor) 


(i)  Ph.D. 

David  Thompson: 

Ronald  Perez: 

Germu  Happawana: 

Sobita  Samanarayake: 


Optimal  Loop  Shapes  in  Quantitative  Feedback  Theory; 
August  1990. 

Integrated  Propulsion-Airframe  Dynamics  and  Control; 
August  1990. 

Dynamics  and  Control  of  Vibration  in  Bladed  Disk 
Assemblies  (Degree  Expected  December  1992) 

The  Algebraic  and  Topological  Dynamics  of  Periodic 
Systems  (Degree  Expected  December  1993) 
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Chin-Homg  Yau: 

Generic  Modelling  and  Control  of  Integrated  Aero 
Engine/Air  Frame  Systems  (Degree  Expected  December 
1992) 

Richard  Nordgren 

Robust  Multivariable-QFT  Rotorcract  Control  (Degree 
Expected  December  1993) 

(ii) 

MS  ME  (Thesis) 

Alley  Butler: 

Analysis  of  the  Control  of  Kamyr  Digester  Based  on  the 
Purdue  Kamyr  Digester  Model;  August  1988. 

John  Bur  hoe: 

Coordinated  Control  of  Contour  and  Tracking  Accuracy 
in  X-Y  Tables;  December  1988. 

Ali  Shahin: 

Modeling  and  Control  of  the  PD  514-1  STOVL  Engine 
(December  1990) 

Joseph  E.  Gallagher 

Multivariable  -  QFT  Control  of  the  Allison  PD  514-1 
STOVL  Engine  (December  1991) 

(iii) 

MSME  (Projects) 

Robert  Matzen: 

Experimental  Modelling  and  Analysis  of  Contour 
Tracking  X-Y  Tables;  May  1988. 

Nader  Motamedi: 

Vibration  Dynamics  of  Bladed  Disk  Assemblies: 

May  1988. 

Christine  Treml: 

Computational  Issues  in  Signal  Processing 
(Supervised  by  P.  Sherman),  May  1988. 

Richard  Walker: 

Controller  Design  for  Coordinated  Control  Applications; 
December  1988. 

Douglas  Starfield: 

Expert  Systems  Control  Software  Development," 

August  1989. 

Way  Ion  Wang 

Active  Control  of  Mechanical  Impedance 

December  1989. 

(iv) 

MSME  (By  Course  Work  Only) 

Richard  Strehler: 

May  1987. 

(v) 

MSEE  (Thesis) 

Gloria  Nnebedum: 

Modelling  and  Control  of  the  Fermentation  Process  in  Garri 
Production,  University  of  Nigeria,  December  1983. 

Alex  Anyaegbunam: 

Microprocessor  Control  of  the  Garri  Frying  Machine, 
University  of  Nigeria,  July  1984. 

Clair  Courtine: 

Nyquist  Array  Based  Methods  for  the  Control  of 
Multivariable  Systems,  Purdue  University  (Co-advisor:  J. 
Chiasson),  May  1988. 
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G.  Teaching  Experience:  September  1975  to  Present. 

Courses  Taught: 

(i)  Undergraduate:  Control  Engineering,  Engineering  Laboratories,  Engineering 
Mathematics,  Modelling  of  Mechanical  and  Electromechanical  Systems,  Electronic 
Circuits,  Digital  Logic  Design,  Electronic  Circuit  Design, 

(ii)  Graduate:  Linear  System  Theory,  Theory  and  Design  of  Control  Systems, 
Multivariable  Control  System  Design. 


(i)  University  of  Nigeria 


ECE512: 

Engr.  301: 
Engr.  401: 
ECE  412: 
EE  312: 

ECE  612: 


Control  Engineering,  1975,  1976, 1977,  1978,  1979. 

(Core  Undergraduate) 

Engineering  Mathematics  1,  1975  —  1984.  (Core  Undergraduate) 
Engineering  Mathematics  2, 1975  —  1984.  (Core  Undergraduate) 
Instrumentation,  1980  -  1984.  (Core  Undergraduate) 

Modelling  of  Mechanical  and  Electromechanical  Systems. 

(Core  Undergraduate) 

Linear  System  Theory.  (Core  Graduate) 


(ii)  Manchester  Polytechnic 

EE  214:  Electronics,  Fall  1979,  Spring  1980.  (Core  Undergraduate) 

EE  314:  Automatic  Control,  Fall  1979,  Spring  1980.  (Core  Undergraduate) 


(iii)  University  of  Notre  Dame 

EE  21 1:  Digital  Logic  Design,  Fall  1984,  Spring  1985.  (Core  Undergraduate) 

EE  31 1:  Electronic  Circuits  Design,  Spring  1985.  (Core  Undergraduate) 


(iv)  Purdue  University 

ME  375:  Modelling  and  Analysis  of  Physical  Systems,  Spring  1987  Fall  1988. 

(Core  Undergraduate) 

ME  463:  Mechanical  Engineering  Design,  Fall,  1988,  Spring  1989 

(Core  Undergraduate) 

ME  475:  Automatic  Control,  Fall  1985,  Spring  1986,  Fall  1986. 

ME  475L:  Mechanical  Engineering  Control  Lab,  Fall  1987,  Spring,  1988,  Spring  1989 

(Core  Undergraduate) 

ME  575:  Theory  and  Design  of  Control  Systems,  Fall  1987. 

(Core  Undergraduate) 

ME  597N:  Multivariable  Control  Systems  Design,  Fall  1986. 

(Graduate) 

ME  697N:  Multivariable  Control  Systems  Design,  Spring  1988. 

(Graduate) 
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H.  New  Courses  Developed 

Reorganized  and  updated  ME  475  laboratory  experiments  with  a  grant 
from  AT&T  Foundation. 

Contributed  to  the  reorganization  of  the  ME  475  course  structure 
and  material. 

Introduced  and  developed  the  new  graduate  course:  ME  697N. 

Multivariable  Control  System  Design 

This  is  the  only  frequency  domain  control  offering  at  Purdue 

I.  Works  under  Review 

1.  Perez,  R.A.,  Nwokah,  O.D.I.,  Thompson,  D.F.,  Almost  decoupling  by  quantitative 
feedback  theory,  Submitted  to  ASME  Journal  Dynamic  Systems  and  Control. 

2.  Nwokah,  O.D.I.,  Jayasuriya,  S.,  Chait,  Y,  Parametric  robust  control  by  quantitative 
feedback  theory.  Submitted  to  AIAA  Journal  Guidance  and  Control. 

3.  Jayasuriya,  S.,  Nwokah,  O.D.I.,  Yaniv,  O.,  The  Benchmark  problem  solution  by 
quantitative  feedback  theory.  (Invited  1991  ACC  Paper)  Also  submitted  for  a 
special  issue  of  AIAA  Journal  Guidance  and  Control. 

4.  Happawana,  G.S.,  Bajaj,  A.K.,  Nwokah,  O.D.I.,  A  singular  perturbation  analysis  of 
eigenvalue  veering  and  mode  localization  in  perturbed  linear  chain  and  cyclic 
systems.  Submitted  to  Journal  of  Sound  and  Vibrations 

J.  Professional  Activities 

1.  Organizer  and  Chairman:  Session  on  Quantitative  Feedback  Theory,  ASME-WAM, 
San  Francisco,  December  1989. 

2.  Organizer  and  Co-chair:  Session  on  QFT  Approach  to  Robust  Control,  ACC,  San 
Diego,  CA,  May  1990. 

3.  Organizer  and  Chair:  Two  Sessions  on  Recent  Developments  in  QFT;  ASME- 
WAM,  Dallas,  TX,  December  1990. 

4.  Chairman:  Session  on  Linear  Stability,  IEEE-CDC,  Honolulu,  Hawaii,  December 

1990. 

5.  Organizer  and  Co-chair:  Session  Parametric  Robust  Control  by  Quantitative 
Feedback  Theory,  ACC,  Boston,  MA,  June  1991. 

6.  Organizer  and  Chair:  2  Sessions  on  QFT,  ASME-WAM,  Atlanta,  GA,  December 

1991. 
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Personal  Data: 
Address: 

Office  Phone: 
Birthdate: 

Education: 

1976-1981 

1974-1976 

1968-1973 

Work  Experience: 
Aug.  1985-present 

June-July  1988 

Feb. -April  1988 


ANIL  KUMAR  BAJAJ 
(updated  January  91) 


School  of  Mechanical  Engineering 

Purdue  University 

West  Lafayette,  IN  47907 

(3 17)  494-6896  Home  Phone:  (317)  463-5487 

April  25,  1951,  Allahabad,  India 


Ph.D.  in  Mechanics,  University  of  Minnesota,  Minneapolis, 
Minnesota.  (Advisor:  Professor  P.  R.  Sethna) 

M.Tech.  in  Mechanical  Engineering,  Indian  Institute  of 

Technology,  Kanpur,  India. 

B.Tech.  in  Mechanical  Engineering,  Indian  Institute  of 

Technology,  Kharagpur,  India. 


Associate  Professor,  School  of  Mechanical 
Engineering,  Purdue  University, 

West  Lafayette,  IN  47907 


Visiting  Associate  Professor,  Department  of  Aerospace 
Engineering  and  Mechanics,  University  of  Minnesota, 
Minneapolis,  MN  55455 


Visiting  Faculty,  Department  of  Mechanical  Engineering 
Indian  Institute  of  Technology,  Kanpur,  India 
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Jan.  198 1-July  1985  Assistant  Professor,  School  of  Mechanical 

Engineering,  Purdue  University, 

West  Lafayette,  IN  47907 


June-August  1981  Visiting  Research  Associate,  Department  of 

Aerospace  Engineering  and  Mechanics, 
University  of  Minnesota, 

Minneapolis,  MN  55455 

Sept.  1976-Dec.  1980  Teaching  and  Research  Assistant, 

Department  of  Aerospace  Engineering 
and  Mechanics,  University  of  Minnesota, 
Minneapolis,  MN  55455 


Aug.  1974-June  1976  Research  Assistant,  Department  of 

Mechanical  Engineering,  Indian  Institute 
of  Technology,  Kanpur,  India 


Honors,  Activities,  Organizations: 


American  Academy  of  Mechanics  -  Member. 

AS  ME  -  Associate  Member 

Sigma  Gamma  Tau  (Aero.  Eng.  Honorary  Society)  -  Member 
B.Tech  -  Institute  Silver  Medal 
Sigma  Xi  -  Member 

Society  of  Engineering  Science  -  Member 

Society  of  Industrial  and  Applied  Mathematics  (SIAM)  -  Member 


Society  Activities: 


Member,  ASME  Technical  Committee  on  Dynamics  of  Structures  and  Systems,  Applied 
Mechanics  Division.  1989-present 
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Recent  Teaching  Experience: 
Purdue  University: 


ME  270  -  Basic  Mechanics  I  (Statics)  -  S81,F82JF83,F84,F85,SS86,F86,F87,F88,F89,F90 

ME  2701  -  Basic  Mechanics  I  (Independent)  -  S82 

ME  274  -  Basic  Mechanics  II  (Dynamics)  -  F81,SS84,S86,SS87,S89 

ME  562  -  Advanced  Dynamics  -  S82,S83,S84,S86,S89,S90  (TV),S91 

ME  563  -  Mechanical  Vibrations  -  F83,S85,SS86,S87,F87 

ME  580  -  Nonlinear  Systems  -  F81,F82,F84,F85,F86,F88,F89,F90 

ME  597T  -  Analysis  and  Design  of  Manipulators  -  S83,S84,S85 

ME  665  -  Dynamic  Stability  of  Elastic  Systems  -  S83 

ME  697B  -  Introduction  to  Bifurcations  and  Chaos  -  S87,S90 


Other  Teaching  Related  Activities: 

1 .  Revised  graduate  courses  ME  562,  ME  580. 

2.  Completely  revised  and  reorganized  ME  665. 

3.  Introduced  new  course:  Introduction  to  Bifurcations  and  Chaos. 


University  and  Depanmental  Activities 

1.  Graduate  School  Fellowship  Committee:  1983,  1984,  1985,  1986,  1987. 

2.  Mechanical  Engineering  Graduate  Committee:  81-82,  82-83,  83-84,  84-85,  85-86. 

3.  Mechanical  Engineering  Curriculum  Committee:  86-87,  87-88. 

4.  Dean’s  Ad  Hoc  Committee  on  Undergraduate  Students  distribution  in  the  Schools  of 
Engineering:  87. 

5.  Schools  of  Engineering  Academic  Personnel  Grievance  Committee:  1988-89,  1989-1990. 

6.  Departmental  Task  Force  on  Secretarial  Services:  1989  Summer. 


Refereeing  Activities: 


ASME  Journal  of  Applied  Mechanics 
SIAM  Journal  on  Applied  Mathematics 

ASME  Journal  of  Dynamic  Systems,  Measurements  and  Control 

ASME  Journal  of  Mechanisms,  Transmissions,  and  Automation  in  Design 

Journal  of  Sound  and  Vibration 

Dynamics  and  Stability  of  Systems 

IEEE  Transactions  on  Circuits  and  Systems 

Journal  of  Fluids  and  Structures 
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Intemaiional  Journal  of  Nonlinear  Mechanics 
Journal  of  Fluid  Mechanics 

ASME  Journal  of  Vibration,  Acoustics,  Stress,  and  Reliability  in  Design 
Physica  D, 

Nonlinear  Dynamics 


Reviewing  Activities: 


National  Science  Foundation,  Division  ot  Mechanics,  Structures  and  Materials. 
Army  Research  Office,  Division  of  Engineering  Sciences. 


Consulting: 


Consultant  to  Rubbermaid  Incorporated,  Reynolds,  IN,  in  structural  re-design  of  a 
molded,  polyester  resin  part  (1989). 

Consultant  to  Arrowhead  Plastic  Engineering,  Inc.,  Muncie,  IN,  in  Failure  Diagnosis  and 
Redesign  of  a  Plastic  Fertilizer  Storage  Box  (1989). 

Consultant  to  O’Reilly,  Cunningham,  Norton  &  Mancini,  Law  Offices,  Wheaton,  IL,  on 
behalf  of  LMP  Steel,  Defendant,  in  a  lawsuit  involving  stacking  of  bundles  of  steel  bars 
constrained  by  crimped  spring  steel  bands  (1989  -  present). 
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Research  and  Other  Grants: 

1.  XL  Summer  Grant  -  Purdue  University  -  1982  -  $3,300. 

2.  Award  from  National  Science  Foundation,  Division  of  Mechanical  Engineering  and 
Applied  Mechanics  for  the  period  September  1982  -  February  1985  for  project  entitled, 
"Bifurcations  in  Rotationally  Symmetric  and  Perturbed  Mechanical  Systems"  -  $49,991. 

3.  XL  Summer  Grant  -  Purdue  University  -  1983  -  $3,500. 

4.  a.  David  Ross  Grant  -  for  the  project  "Dynamic  Analysis  of  Robotic  Manipulators  With 
Parametric  and  External  Excitation",  -  1984  (with  C.  M.  Krousgrill)  -  $6,600. 

b.  Renewal  of  above,  1985  $6,600. 

5.  PRF  International  Travel  Grant  -  Purdue  University  -  1984  -  $800. 

6.  Travel  Grant-in-aid-U.S.N.C.  on  Theoretical  and  Applied  Mechanics,  National  Research 
Council  -  1984  -  $700. 

7.  AMP  Inc.  "Study  of  Parametric  Instability  of  High  Speed  Elastic  Mechanisms",  September 
1984  -  August  1985  (with  A.  Midha)  -  $12,000. 

8.  a.  David  R^ss  Grant  -  for  the  project  Qualitative  Studies  in  the  Nonlinear  Dynamics  of 
Compliant  Offshore  Struc  tres"  -  1985  (with  C.  M.  Krousgrill)  -  $6,600. 

b.  Renewal  of  above,  1986  -  $6,600. 

9  Govt,  of  Indonesia  grant  for  the  project  "M.S.  -  Mechanical  Engineering  -  Andi 
Mahyuddin"  August  1986  -  July  1988  (with  Mr.  A.  Mahyuddin)  -  $13,458. 

10.  Travel  Grant  -  U.  S.  National  Committee  on  Theoretical  and  Applied  Mechanics,  National 
Research  Council  -  1988  -  $1,025. 

11.  Award  from  AFOSR,  Structural  Mechanics  Program,  Bolling  AFB,  Washington,  D.C.  for 
November  1988  -  October  1990  for  project  "Vibration  Dynamics  and  Control  of  Bladed 
Disk  Assemblies"  -  $124,049  (Co-PI  with  Professor  O.  D.  Nwokah). 

12.  a.  David  Ross  Grant  -  for  the  project  "Analytical  and  Experimental  Investigations  of 
Chaotic  Response  in  Structures  with  Nonlinear  Modal  Interactions"  -  1989  (with  C.  M. 
Krousgrill)  -  $7,890. 

b.  Renewal  of  above,  1990  -  $9,000. 

13.  Dean’s  Club  "Creative  Undergraduate  Instructional  Project”  grant  for  the  project  - 
Development  of  Graphical  Display  and  Simulation  Software  for  the  Enhancement  of 
Visualization  in  Undergraduate  Mechamcs  Education  -  1990  (with  C.  M.  Krousgrill  and 
D.  C.  Anderson)  $9,000. 

14.  Hewlett  -  Packard  Equipment  Grant  (through  the  Head  of  the  School)  for  instructional 
laooratory  in  dynamics  of  mechanical  systems  -  1990  (with  C.  M.  Krousgrill  and  A. 
Midha)  $39,735. 


-48- 


15.  Award  from  Engineering  Sciences  Division,  Army  Research  Office,  Dept,  of  the  Army  for 
the  period  November  1,  1990  -  October  31,  1993  for  project  "Modal  Interactions  and 
Complex  Responses  in  Weakly  Nonlinear  Multi-Degree-of-Freedom  Mechanical 
Systems"  -  $290,000  (PI  with  Prof.  Davies  as  Co-PI). 
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M.S.  Theses  Directed: 

1.  S.  Tousi,  School  of  Mechanical  Engineering:  Period-Doubling  Bifurcations  and 
Modulated  Motions  in  Forced  Mechanical  Systems  (1984). 

2.  A.  I.  Mahyuddin,  School  of  Mechanical  Engineering:  A  Comparative  Parametric  Stability 
Study  for  Flexible  Cam-Follower  Systems  due  to  Varying  Cam  Profiles  (Dec.  1988)  (  Co- 
Advisor  with  Prof.  A.  Midha). 


M.S.  Students  Advised  (Non-thesis  Option) 

1.  Thomas  Berglund,  Tod  Dalrymple,  Carol  Dittrich,  David  Ehinger,  Leslie  Grundman, 
Gerald  Merrill,  Scott  Pantaleo,  Timothy  Rorick,  David  Smith,  David  W.  Turner 


Ph.D.  Theses  Directed: 

1.  M.  Ilkhani-Pour,  School  of  Mechanical  Engineering:  Effect  of  a  Viscoelastic  Boundary 
Inset  on  Vibrational  Characteristics  of  a  Beam  (Dec.  1984)  (Co-Advisor  with  Prof.  E.  C. 
Ting,  CE). 

2.  D.  A.  Streit,  School  of  Mechanical  Engineering:  Dynamics  and  Stability  of  Robotic 
Manipulators  with  Parametric  Excitation  (Aug.  86)  (Co-Advisor  with  Prof.  C.  M. 
Krousgrill). 

3.  Y.  M.  Huang,  School  of  Mechanical  Engineering:  Nonlinear  Dynamics  of  Low  Order 
Models  of  Offshore  Structures  (Dec.  87)  (Co- Advisor  with  Prof.  C.  M.  Krousgrill). 

4.  S.  Tousi,  School  of  Mechanical  Engineering:  Directional  Stability  of  Road  Vehicles  (Dec. 
88)  (Co-Advisor  with  Prof.  W.  Soedel). 

5.  J.  M.  Johnson,  School  of  Mechanical  Engineering:  On  the  Vibration  of  Stretched  Strings 
(May  1989). 
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Joumal  Publications: 

1.  Bajaj,  A.  K.  and  Garg,  V.  K.,  "Linear  Stability  of  Jet  Rows,"  J.  Appl.  Mech.,  Vol.  44, 

Sept.  1977,  pp.  378-384.  ~ 

2.  Sethna,  P.  R.  and  Bajaj,  A.  K.,  "Bifurcations  in  Dynamical  Systems  with  Internal 
Resonance,"  J.  Appl.  Mech.,  Vol.  45,  Dec.  1978,  pp.  895-902. 

3.  Lundgren,  T.  S.,  Sethna,  P.  R.,  and  Bajaj,  A.  K.,  "Stability  Boundaries  for  Flow  Induced 
Motions  of  Tubes  with  an  Inclined  Terminal  Nozzle,”  J.  Sound  and  Vib.,  Vol.  64,  1979, 
pp.  553-571. 

4.  Bajaj,  A.  K.,  Sethna,  P.  R.,  and  Lundgren,  T.  S.,  "Hopf  Bifurcation  Phenomena  in  Tubes 
Carrying  a  Huid,"  SIAM  J.  Appl.  Math.,  Vol.  39,  1980,  pp.  213-230. 

5.  Bajaj,  A.  K.,  "Bifurcating  Periodic  Solutions  in  Rotationally  Symmetric  Systems," 
SIAM  J.  Appl.  Math.,  Vol.  42,  Oct.  1982,  pp.  1078-1098. 

6.  Bajaj,  A.  K.  and  Sethna,  P.  R.,  "Bifurcations  in  Three-Dimensional  Motions  of  Articulated 
Tubes,  Part  1:  Linear  Systems  and  Symmetry,"  J.  Appl.  Mech.,  Vol.  49,  Sept.  1982,  pp. 
606-611. 

7.  Bajaj,  A.  K.  and  Sethna,  P.  R.,  "Bifurcations  in  Three-Dimensional  Motions  of  Articulated 
Tubes,  Part  2:  Nonlinear  Analysis,"  J.  Appl,  Mech.,  Vol.  49,  Sept.  1982,  pp.  612-618. 

8.  Bajaj,  A.  K.  and  Sethna,  P.  R.,  "Row-Induced  Bifurcations  to  Three-Dimensional 

Oscillatory  Motions  in  Continuous  Tubes,"  SIAM  J.  Appl.  Math,  Vol.  44,  Apr.  1984,  pp. 
270-286.  "  ~~ 

9.  Bajaj,  A.  K.,  "Interactions  Between  Self  and  Parametrically  Excited  Motions  in 
Articulated  Tubes,"  J.  Appl.  Mech.,  Vol.  51,  June  1984,  pp.  423-429. 

10.  Sethna,  P.  R.,  Meyer,  K.  R.,  and  Bajaj,  A.  K.,  "On  the  Method  of  Averaging,  Integral 
Manifolds  and  Systems  with  Symmetry,"  SIAM  J.  Appl.  Math.,  Vol.  45,  June  1985,  pp. 
343-359. 

11.  Tousi,  S.  and  Bajaj,  A.  K.,  "Period-Doubling  Bifurcations  and  Modulated  Motions  in 
Forced  Mechanical  Systems,"  J.  Appl.  Mech.,  Vol.  107,  June  1985,  pp.  446-452. 

12.  Bajaj,  A.  K.,  "Resonant  Parametric  Perturbations  of  the  Hopf  Bifurcation,"  J. 
Math.  Anal.  &  Appl,,  Vol.  1 15,  April  1986,  pp.  214-224. 

13.  Bajaj,  A.  K„  "Bifurcations  in  a  Parametrically  Excited  Nonlinear  Oscillator,"  Int. 
J.  Nonlinear  Mech.,  Vol.  22,  1987,  pp.  47-59. 

14.  Streit,  D.  A.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "A  Preliminary  Study  of  the  Dynamic 

Stability  of  Flexible  Manipulators  Performing  Repetitive  Tasks",  ASME 

J.  Dynamic  Syst.,  Meas.  and  Cont.,  Vol.  108,  Sept.  1986,  pp.  206-214. 


15.  Bajaj,  A.  K.,  "Nonlinear  Dynamics  of  Tubes  Carrying  a  Pulsatile  Row", 
Dynamics  and  Stability  of  Systems,  Vol.  2,  1987,  pp.  19-41. 
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16.  Farhang,  K.,  Midha,  A.  and  Bajaj,  A.  K.,  "A  Higher-Order  Analysis  of  Basic  Linkages  for 
Harmonic  Motion  Generation",  ASME  J.  Mechanisms,  Transmissions,  and  Automation 
in  Pes.,  Vol.  109,  Sept.  1987,  pp.  301-307. 

17.  Farhang,  K.,  Midha,  A.  and  Bajaj,  A.  K.,  "Synthesis  of  Harmonic  Motion  Generating 
Linkages,  Part  I:  Function  Generation",  ASME  J.  Mechanisms,  Transmissions,  and 
Automation  in  Des.,  Vol.  1 10,  March  1988,  pp.  16-21. 

18.  Heiman,  M.  S.,  Sherman,  P.  J.  and  Bajaj,  A.  K.,  "On  the  Dynamics  and  Stability  of  an 
Inclined  Impact  Pair",  J.  Sound  and  Vibration,  Vol.  1 14,  No.  3,  1987,  pp.  535-547. 

19.  Streit,  D.  A.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Nonlinear  Response  of  Flexible  Robotic 

Manipulators  Performing  Repetitive  Tasks”,  ASME  J.  Dynamic  Syst.,  Meas.  and  Contr., 
Vol.  1 1 1,  No.  3,  1989,  pp.  470-480.  ~~  ~ "  ""  —  “ 

20.  Huang,  Y.  M.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Nonlinear  Response  of  a  Dynamic 

System  due  to  Oscillatory  Flow",  ASME  J.  Offshore  Mech.  and  Arct.  Engr.,  Vol.  109,  No. 
4,  1987,  pp.  345-356.  '  ~ 

21.  Streit,  D.  A.,  Bajaj,  A.  K.,  and  Krousgrill,  C.  M.,  "Combination  Parametric  Resonance 
Leading  to  Periodic  and  Chaotic  Response  in  Two  Degree-of-Freedom  Systems  with 
Quadratic  Nonlinearities",  J.  Sound  and  Vibration,  Vol.  124,  No.  2,  July  1988,  pp.  297- 
314. 

22.  Tousi,  S.,  Bajaj,  A.  K.  and  Soedel,  W.,  "On  the  Stability  of  a  Flexible  Vehicle  Controlled 
by  a  Human  Pilot",  Vehicle  System  Dynamics,  Vol.  17, 1988,  pp.  37-56. 

23.  Heiman,  M.  S.,  Bajaj,  A.  K.  and  Sherman,  P.  J.,  "Periodic  Motions  and  Bifurcations  in 
Dynamics  of  an  Inclined  Impact  Pair",  J.  Sound  and  Vibration,  Vol.  124,  No.  1,  July  1988, 
pp. 55-78. 

24.  Johnson,  J.  M.  and  Bajaj,  A.  K.,  "Amplitude  Modulated  and  Chaotic  Dynamics  in 
Resonant  Motion  of  Strings",  J.  Sound  and  Vibration,  Vol.  128,  January  1989,  pp.  87-107. 

25.  Huang,  Y.  M.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Dynamic  Behavior  of  Offshore 
Structures  with  Bilinear  Stiffness",  J.  Fluids  and  Structures,  Vol.  3,  July  1989,  pp.  405- 
422. 

26.  Bajaj,  A.  K.  and  Johnson,  J.  M.,  "Asymptotic  Techniques  and  Complex  Dynamics  in 
Weakly  Nonlinear  Forced  Mechanical  Systems",  Int.  J.  Non-Linear  Mech,  Vol.  25,  No. 
2/3,  1990,  pp.  21 1-226. 

27.  Huang,  Y.  M.,  Bajaj,  A.  K.  and  Krousgrill,  C.  M.,  "Complex  In-Line  and  Whirling 

Response  of  Structures  to  Oscillatory  Flow”,  J.  Sound  and  Vibration,  Vol.  136,  No.  3, 
1990,  pp.  491-505.  “ 

28.  Tousi,  S.,  Bajaj,  A.  K.  and  Soedel,  W.,  "Finite  Disturbance  Closed-Loop  Directional 
Stability  of  Vehicles  with  Human  Pilot  Considering  Nonlinear  Cornering  Behavior", 
Vehicle  System  Dynamics  (to  appear). 
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29.  Bajaj,  A.  K.  and  lousi,  S.,  "Torus  Doublings  and  Chaotic  Amplitude  Modulations  in  a 

Two  Degree-of-Freedom  Resonantly  Forced  Mechanical  System",  Int.  J.  Non-Lin.  Mech., 
Vo.  25,  No.  6,  1990,  pp.  625-642.  '  " 

30.  Tousi,  S.,  Bajaj,  A.  K.  and  Soedel,  W.,  "Closed-Loop  Directional  Stability  of  Car-Trailer 
Combinations",  Vehicle  System  Dynamics  (submitted). 

31.  Mahyuddin,  A.  I.,  Midha,  A.  and  Bajaj,  A.  K.,  "On  Methods  for  Evaluation  of  Parametric 
Stability  and  Response  of  Flexible  Cam-Follower  Systems",  ASME  J.  Mech.  Design 
(submitted). 

32.  Restuccio,  J.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Nonlinear  Nonplanar  Dynamics  of  a 
Parametrically  Excited  Inextensional  Elastic  Beam",  Nonlinear  Dynamics  (to  appear). 

33.  Bajaj,  A.  K.  and  Sethna,  P.  R.,  "Effect  of  Symmetry-Breaking  Perturbations  on  Flow- 
Induced  Oscillations  in  Tubes",  J.  of  Fluids  and  Structures  (submitted). 

34.  Happawana,  G.,  Bajaj,  A.  K.  and  Nwokah,  O.D.I.,  "A  Singular  Perturbation  Perspective  on 
Mode  Localization",  J.  Sound  and  Vibration  (to  appear). 

35.  Bajaj,  A.  K.  and  Johnson,  J.  M.,  "On  the  Amplitude  Dynamics  and  "Crisis"  in  Resonant 
Motion  of  Stretched  Strings",  Proceedings  of  the  Royal  Society  of  London  (submitted). 

36.  Happawana,  G.,  Bajaj,  A.  K.  and  Nwokah,  O.D.I.,  "A  Singular  Perturbation  Analysis  of 
Eigenvalue  Veering  and  Mode  Localization  in  Perturbed  Linear  Chain  and  Cyclic 
Systems",  J.  Sound  and  Vibration  (submitted). 


Conference  Publications: 

1.  Sethna,  P.  R.  and  Bajaj,  A.  K.,  "Bifurcation  in  Nonlinear  Oscillatory  Systems,"  Proc.  of 
VHIth  Int.  Conference  on  Nonlinear  Oscillations,  Prague,  1978,  pp.  621-626. 

2.  A.  K.  Bajaj,  "Primary  Parametric  Resonance  in  a  Self-Excited  Oscillator,"  Proceedings  of 
the  XVIHth  Midwestern  Mechanics  Conference,  The  University  of  Iowa,  1983. 

3.  A.  K.  Bajaj,  "Some  Interactions  in  Flow-Induced  Motions  of  Tubes,"  Proceedings  of  the 
ASCE,  Engineering  Mechanics  Specialty  Conference,  Purdue  University,  1983. 

4.  Streit,  D.  A.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Dynamic  Stability  of  Flexible 
Manipulators  Performing  Repetitive  Tasks",  in  Robotics  and  Manufacturing  Automation, 
PED-Vol.  15,  eds.  M.  Donath,  M.  Leu,  ASME,  New  York,  1985. 

5.  Heiman,  M.  S.,  Sherman,  P.  J.  and  Bajaj,  A.  K.,  "On  the  Dynamics  and  Stability  of  an 
Inclined  Impact  Pair",  Proceeding  of  the  27th  AIAA/ASME/ASCE/AHS  Structures, 
Structural  Dynamics  and  Materials  Conference,  San  Antonio,  TX,  May  19-21,  1986,  pp. 
585-591.  Paper  No.  AIAA-86-0998-CP. 

6.  Farhang,  K.,  Midha,  A.  and  Bajaj,  A.  K.,  "A  Higher-Order  Analysis  of  Basic  Linkages  for 
Harmonic  Motion  Generation",  ASME  Paper  86-DET-70,  19th  Biennial  ASME 
Mechanisms  Conference,  Columbus,  Ohio,  October  1986. 
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7.  Bajaj,  A.  K.  and  Krousgrill,  C.  M.,  "Complex  Nonlinear  Dynamics  of  a  Rotating 
Pendulum",  Proceedings  of  the  20th  Midwestern  Mechanics  Conference,  Purdue 
University,  Aug.  31  -  Sept.  2,  1987,  pp.  832-839. 

8.  Johnson,  J.  M.  and  Bajaj,  A.  K.,  "Chaotic  Amplitude  Modulations  in  Forced  Vibrations  of 
Strings",  Proceedings  of  the  20th  Midwestern  Mechanics  Conference,  Purdue  University, 
Aug.  31  -  Sept.  2,  1987,  pp.  851-856. 

9.  Tousi,  S.,  Bajaj,  A.  K.  and  Soedel,  W.,  "On  Stability  of  Flexible  Vehicle  Controlled  by  a 
Human  Pilot",  Proceedings  of  the  20th  Midwestern  Mechanics  Conference,  Purdue 
University,  Aug.  31  -  Sept  2,  1987,  pp.  1354-1359. 

10.  Huang,  Y.  M.,  Krousgrill,  C.  M.,  and  Bajaj,  A.  K.,  "Subharmonic  Resonances  in  a  Bilinear 
Dynamic  System  due  to  Oscillatory  Flow",  Proceedings  of  the  20th  Midwestern 
Mechanics  Conference,  Purdue  University,  Aug.  31  -  Sept.  2,  1987,  pp.  637-642. 

11.  Bajaj,  A.  K.,  "Bifurcations  in  Flow-Induced  Oscillation  in  Tubes  Carrying  a  Fluid", 
Dynamical  Systems  Approaches  to  Nonlinear  Problems  in  Systems  and  Circuits,  eds.  Fathi 
M.  A.  Salam  and  Mark  L.  Levi,  SIAM,  Philadelphia,  1988. 

12.  Farhang,  K.,  Midha,  A.  and  Bajaj,  A.  K.,  "Synthesis  of  Harmonic  Motion  Generating 
Linkages,  Part  I:  Function  Generation",  Advances  in  Design  Automation  -  1987,  DE-Vol. 
10-2,  ed.  S.S.  Rao,  ASME,  New  York,  1987. 

13.  Streit,  D.  A.,  Bajaj,  A.  K.  and  Krousgrill,  C.  M.,  "Parametric  Instabilities  and  Chaotic 
Amplitude  Modulations  in  a  System  with  Two  Degrees-of-Freedom",  Proceedings  of  3rd 
International  Conf.  on  Recent  Advances  in  Structural  Dynamics,  University  of 
Southampton,  Southampton,  July  1988,  pp.  809-818. 

14.  Huang,  Y.  M.,  Krousgrill,  C.  M.  and  Bajaj,  A.  K.,  "Dynamic  Behavior  of  Offshore 
Structures  with  Bilinear  Stiffness",  1988  International  Symposium  on  Flow-Induced 
Vibration  and  Noise,  Vol.  7  -  Nonlinear  Interaction  Effects  and  Chaotic  Motions,  eds:  M. 
M.  Reischman,  M.  P.  Paidoussis,  R.  J.  Hansen,  ASME,  NY  1988. 

15.  Zadoks,  R.  I.  and  Bajaj,  A.  K.,  "Mel'nikov’s  Criterion  and  Chaos  in  a  Model  of  an  Arch", 
Proceedings  of  the  First  Pan-American  Congress  of  Applied  Mechanics,  Rio  de  Janeiro, 
Brazil,  January  3-6,  1989. 

16.  Mahyuddin,  A.  I.,  Midha,  A.  and  Bajaj,  A.  K.,  "On  Methods  for  Evaluation  of  Parametric 
Stability  and  Response  of  Flexible  Cam-Follower  Systems”,  ASME  1990  Mechanisms 
Conference,  Chicago,  IL,  Sept.  15-18,  1990. 

17.  Chang,  S.  I.,  Bajaj,  A.  K.  and  Krousgrill,  C.  M.,  "Amplitude  Modulated  Dynamics  in 
Harmonically  Excited  Non-Linear  Oscillations  of  Rectangular  Plates  with  Internal 
Resonance",  Accepted  for  4th  Int.  Conference  on  Recent  Advances  in  Structural 
Dynamics,  Inst,  of  Sound  and  Vibration  Research,  Southampton,  England,  July  1991. 

18.  Happawana,  G.,  Bajaj,  A.  K.  and  Nwokah,  O.D.I.,  "Mode  Localization  in  Structural 

Systems:  A  Singular  Perturbation  Analysis",  Submitted  for  32nd 

AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural  Dynamics  and  Materials 
Conference,  Baltimore,  Maryland,  April  8-10,  1991. 
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19.  Bajaj,  A.  K.,  "Examples  of  Crisis  Phenomenon  in  Structural  Dynamics",  Accepted  for 
Proceedings  of  the  1990  Conference  on  Bifurcations  and  Chaos;  Theory,  Algorithms  and 
Applications,  Wiirzburg,  W.  Germany. 

20.  Happawana,  G.,  Bajaj,  A.  K.  and  Nwokah,  O.D.I.,  "A  Singular  Perturbation  Analysis  of 
Eigenvalue  Veering  and  Mode  Localization  in  Perturbed  Linear  Chain  and  Cyclic 
Systems",  Submitted  for  the  13th  Biennial  ASME  Conference  on  Mechanical  Vibration 
and  Noise,  Miami,  Florida,  September  22-25,  1991. 

21.  Happawana,  G.,  Nwokah,  O.D.I.  and  Bajaj,  A.  K.,  "On  the  Dynamics  of  Perturbed 
Symmetric  Systems",  Submitted  for  the  13th  Biennial  ASME  Conference  on  Mechanical 
Vibration  and  Noise,  Miami,  Florida,  September  22-25,  1991. 

22.  Mohamad,  A.A.  and  Bajaj,  A.  K.,  "Dynamics  of  a  Spinning  Toroidal  Thermosyphon", 
Accepted  for  the  1991  ASME/AlChE  National  Heat  Transfer  Conference,  Minneapolis, 
Minnesota,  July  28-31,  1991. 
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Conference  Presentations: 

1.  "Bifurcations  in  Dynamical  Systems  with  Internal  Resonance,"  P.  R.  Sethna,  A.  K.  Bajaj, 
ASME  Winter  Annual  Meeting,  San  Francisco,  CA,  December  10-15,  1978. 

2.  "Bifurcations  in  Three-Dimensional  Motions  of  Articulated  Tubes,"  A.  K.  Bajaj,  P.  R. 
Sethna,  EXth  U.S.  National  Congress  of  Applied  Mechanics,  Cornell  University,  Ithaca, 
NY,  June  21-25,  1982. 

3.  "Bifurcations  in  a  Parametrically  Excited  Nonlinear  Oscillator",  A.  K.  Bajaj,  16th  ICTAM 
held  at  Lyngby,  Denmark,  August  19-25,  1984. 

4.  "Period  Doubling  Bifurcations  and  Modulated  Motions  in  Forced  Mechanical  Systems",  S. 
Tousi,  A.  K.  Bajaj,  21st  Annual  Meeting  of  the  Society  of  Engineering  Science,  VPISU, 
Blacksburg,  Virginia,  October  15-17,  1984. 

5.  "Nonlinear  Dynamics  of  Tubes  Carrying  a  Pulsatile  Flow",  A.  K.  Bajaj,  21st  Annual 
Meeting  of  the  Society  of  Engineering  Science,  VPISU,  Blacksburg,  Virginia,  October 
15-17,  1984. 

6.  "Parametric  Excitation,  Period  Doubling  and  Chaos  in  Two-Degree-of-Freedom  Flexible 
Mechanical  Systems",  D.  A.  Streit,  A.  K.  Bajaj,  C.  M.  Krousgrill,  10th  U.S.  National 
Congress  of  Applied  Mechanics,  The  University  of  Texas  at  Austin,  TX,  June  16-20,  1986. 

7.  "Nonlinear  Dynamics,  Resonances  and  Chaotic  Motions  in  Harmonically  Excited 
Mechanical  Systems",  A.  K.  Bajaj,  23rd  Annual  Meeting  of  the  Society  of  Engineering 
Science,  SUNY,  Buffalo,  N.Y.,  August  25-27,  1986. 

8.  "Bifurcations  in  Flow  Induced  Oscillations  in  Tubes  Carrying  a  Fluid",  A.  K.  Bajaj, 
Engineering  Foundation  Conference  on  "Qualitative  Methods  for  the  Analysis  of 
Nonlinear  Dynamics",  Henniker,  N.H.,  June  8-13,  1986. 

9.  "Stability  and  Nonlinear  Dynamics  of  Flexible  Manipulators  Performing  Repetitive 
Tasks",  A.  K.  Bajaj,  C.  M.  Krousgnii,  D.  A.  Streit,  AFOSR/ARO  Conference  on  Non- 
Linear  Vibrations  Stability,  and  Dynamics  of  Structures  and  Mechanisms,  VPISU, 
Blacksburg,  March  23-25,  1987. 

10.  "Nonlinear  Dynamic  Response  of  a  Structure  Due  to  Oscillatory  Flow”,  Y.  M.  Huang,  C. 
M.  Krousgrill  and  A.  K.  Bajaj,  1 1th  Canadian  Congress  of  Applied  Mechanics,  Edmonton, 
May  31  -  June  4,  1987. 

11.  "Periodic  Motions  and  Bifurcations  in  Dynamics  of  an  Inclined  Impact  Pair",  M.  S. 
Heiman,  P.  J.  Sherman  and  A.  K.  Bajaj,  11th  Canadian  Congress  of  Applied  Mechanics, 
Edmonton,  May  31  -  June  4,  1987. 

12.  "Flow-Induced,  Bi-Planar  Motion  of  Marine  Structures",  Y.  M.  Huang,  A.  K.  Bajaj,  C.  M. 
Krousgrill,  Applied  Mechanics  and  Engineering  Science  Conference,  Berkeley.  June  20- 
22,  1988. 

13.  "Amplitude  Modulated  Chaos  in  Harmonically  Excited  Mechanical  Systems",  A.  K.  Bajaj, 
J.  M.  Johnson,  SIAM  1988  Annual  Meeting,  Minneapolis,  Minnesota,  July  11-15,  1988. 
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14.  "Chaotic  Amplitude  Dynamics  of  an  Autoparametric  Two-Degree-of-Freedom  System", 
A.  K.  Bajaj,  J.  M.  Johnson,  Accepted  for  IUTAM  Symposium  on  Nonlinear  Dynamics  in 
Engineering  systems.  University  of  Stuttgart,  Stuttgart,  FRG,  August  21-25,  1989. 

15.  "Nonlinear  Dynamics  of  a  Parametrically  Excited  Inextensional  Elastic  Beam",  J.  M. 
Restuccio,  C.  M.  Krousgrill  and  A.  K.  Bajaj,  SIAM  Conference  on  Dynamical  Systems, 
Orlando,  FL„  May  7-10,  1990. 

16.  "On  the  Dynamics  of  Perturbed  Symmetric  Systems",  G.  Happawana,  O.D.I.  Nwokah,  D. 
Afolabi  and  A.  K.  Bajaj,  11th  U.S.  National  Congress  of  Applied  Mechanics,  Univ.  of 
Arizona,  Tucson,  May  21-25,  1990. 


Invited  Conference  Presentations 

1.  "Asymptotic  Techniques  and  Chaos  in  Weakly  Nonlinear  Forced  Mechanical  Systems",  A. 
K.  Bajaj,  Non-Linear  Vibrations,  Stability,  and  Dynamics  of  Structures  and  Mechanisms 
Conference,  VPISU,  Blacksburg,  June  1-3,  1988. 

2.  "Torus  Doubling  and  Chaotic  Amplitude  Modulations  in  a  Forced  Mechanical  System",  A. 
K.  Bajaj  and  S.  Tousi,  Inter.  Conf.  on  Dynamical  Systems,  Control  Theory  and 
Applications,  Wright  State  University,  Dayton,  June  15-17,  1989. 

3.  "Examples  of  Crisis  Phenomena  in  Structural  Dynamics",  A.  K.  Bajaj,  Joint  ASCE/ASME 
Mechanics  Conference,  San  Diego,  July  10-12,  1989. 

4.  "On  The  Complex  Whirling  Dynamics  of  Strings",  A.  K.  Bajaj,  Third  Conference  on 
Nonlinear  Vibrations,  Stability,  and  Dynamics  of  Structures  and  Mechanisms,  VPISU, 
Blacksburg,  June  25-27,  1990. 

5.  "On  the  Complex  Whirling  Dynamics  of  Strings",  A.  K.  Bajaj,  Conference  on  Bifurcation 
and  Chaos:  Analysis,  Algorithms  and  Applications,  Wurzburg,  West  Germany,  August  20- 
24,  1990. 


Invited  Talks: 

1.  "Bifurcations  in  Systems  with  Symmetry,"  given  in  Control  Science  and  Dynamical  System 
Center,  University  of  Minnesota,  Minneapolis,  MN,  January  29,  1981. 

2.  "Resonant  Parametric  Perturbations  of  the  Hopf  Bifurcation,"  given  at  the  University  of 
Missouri-Rolla,  Rolla,  MO  65401,  on  February  25,  1983,  jointly  sponsored  by  the 
Departments  of  Engineering  Mechanics  and  Mechanical  Engineering. 

3.  "Some  Examples  of  Complex  Dynamics  in  Mechanical  Systems",  given  in  Department  of 
Mechanical  Engineering,  I.I.T.  Kanpur,  India,  on  August  5,  1985. 

4.  "Internal  Resonance  and  Chaos  in  Externally  Excited  Mechanical  Systems"  given  in 
Department  of  Aeronautical  and  Astronautical  Engineering,  University  of  Illinois  at 
Urbana-Champaign,  November  24,  1986. 
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5.  "Chaotic  Amplitude  Modulations  in  Nonlinear  Systems  -  Examples  from  Mechanical 
Oscillators",  given  in  School  of  Chemical  Engineering,  Purdue  University,  West  Lafayette, 
October  1,  1987. 

6.  "Chaotic  Amplitude  Modulations  in  Harmonically  Excitated  Weakly  Nonlinear  Systems", 
given  at  Indian  Institute  of  Technology,  Kanpur,  India,  March  30,  1988. 

7.  "Introduction  to  Bifurcations  and  Chaos",  given  at  Indian  Institute  of  Technology, 
Kharagpur,  India,  April  4,  1988. 

8.  "Chaotic  Dynamics  of  Impacting  Systems",  given  at  Indian  Institute  of  Technology,  Kanpur, 
India,  April  7,  1988. 

9.  On  the  Complex  Resonant  Motions  of  Structures  and  Whirling  Dynamics  of  Strings"  given 
at  Department  of  Mechanical  Engineering,  University  of  Pittsburgh,  Pittsburgh,  October  23, 
1990. 
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Book  Reviews: 


Non-Linear  Oscillations,  by  P.  Hagedom,  Oxford  Univ.  Press,  Mechanism  and  Machine  Theory 
Vol.  20,  No.  3,  1985,  pp.  243. 


Others: 

1.  Attended  Conference  on  "New  Approaches  to  Nonlinear  Problems  in  Dynamics,”  sponsored 
by  Engineering  Foundation  at  Asilomar  Conference  Grounds,  Pacific  Grove,  CA,  December 
9-14,  1979. 

2.  Attended  Workshop  on  "NSF  Research  Needs  in  Theoretical  Foundation  of  Dynamics", 
State  University  of  New  York  at  Buffalo,  August  24,  1986. 

3.  Session  Chairman,  "Chaotic  Motions”  in  20th  Midwestern  Mechanics  Conference,  Purdue 
University,  Aug.  31-Sept.  1,  1987. 

4.  Member,  Organizing  Committee,  20th  Midwestern  Mechanics  Conference,  Purdue 
University,  Aug.  31  -  Sept.  1,  1987. 

5.  Participated  in  Presenting  four-day  short  course  "Applied  Digital  Signal  Processing",  Purdue 
University,  October  1987. 

6.  Attended  Short  Course  "Nonlinear  Dynamics,  Chaos,  and  Bifurcation",  SIAM  Annual 
Meeting,  Minneapolis,  July  10,  1988. 

7.  Participated  in  Presenting  five-day  short  course"  Applied  Digital  Signal  Processing",  Purdue 
University,  June  26-30,  1989. 

8.  Participated  in  Presenting  five-day  short  course  "Applied  Digital  Signal  Processing",  Purdue 
University,  June  18-22,  1990. 

9.  Session  Chairman,  "Methods  I"  at  the  "3rd  Conference  on  Nonlinear  Vibrations,  Stability, 
and  Dynamics  of  Structures  and  Mechanisms",  VPISU,  Blacksburg,  June  25-27,  1990. 


Other  Publications: 


Midha,  A.  annd  Bajaj,  A.  K.,  "The  Convolution  Integral  Solution  of  Linear  Spring-Mass-Damper 
Vibration  System  -  A  Versatile  Tool",  Mech.  Eng’g.  News,  Vol,  23,  No.  3,  1986,  pp.  5-15. 
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General  Record 


1.  Education 

PhD  1983  (Mechanical  Engineering)  Imperial  College,  London,  England. 

MSc  1978  (Mechanical  Engineering)  Imperial  College,  London,  England. 

BSc  1976  (Mechanical  Engineering)  Thames  Polytechnic,  London,  England. 

2.  Academic  Appointments 

Visiting  Scientist,  Summer  1990 

Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts. 

Visiting  Scientist,  Summer  1989 

Institute  for  Computational  Mechanics  in  Propulsion, 

NASA  Lewis  Research  Center,  Cleveland,  Ohio. 

Associate  Professor  of  Mechanical  Engineering,  July 1989-present 
Indiana  University-Purdue  University  at  Indianapolis. 

Assistant  Professor  of  Mechanical  Engineering,  September  1985-June  1989 
Indiana  University-Purdue  University  at  Indianapolis. 

Research  Associate,  April  1984- August  1985 

University  of  Maryland,  College  Park,  Maryland. 

Research  Student,  September  1979-December  1982 
Imperial  College,  London,  England. 

3.  Awards  and  Honors 

Rolls-Royce  Scholar  at  Imperial  College,  1979-1982. 

Awarded  a  Mark  of  Distinction  by  the  University  of  London  Senate,  1978. 

Awarded  BP  (British  Petroleum  Co.)  Scholarship  at  Thames  Polytechnic,  London, 
1972-1976. 

4.  Membership  in  Professional  Societies 

Member,  American  Institute  of  Aeronautics  and  Astronautics  (ALAA). 

Member,  American  Society  of  Mechanical  Engineers  (ASME). 

Member,  Society  for  Experimental  Mechanics  (SEM). 

Member,  International  Gas  Turbines  Institute. 
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Research  Record 


1.  Research  Interests 

Vibration  of  Structures  (including  large  scale  structures)  •  Classical  Dynamics  • 
Dynamical  Systems  (including  Bifurcation  of  Vector  Fields  and  Chaos)  •  Singularity 
Theory  and  Catastrophe  Theory  •  Computer  Aided  Testing  (including  Experimental 
Modal  Analysis)  •  Parallel  Computation  (mainframe — CrayXMP,VMP — and  desktops 
computers  —  Transputers,  Hypercube  systems,  etc)  •  Applied  Mathematics  and 
Numerical  Methods  in  Engineering  (including  Finite  Element  Method). 


2.  Publications 

a.  Articles  In  Journals,  Monographs  and  Books 

Afolabi,  D.,  and  Crandall,  S.  H.,  1991,  "The  Weyl  Groups  Ak,  Dk,  Ek  and  the  Random 
Vibration  of  Symmetric  Structures"  (forthcoming). 

Nwokah,  O.  D.  I.,  Afolabi,  D.,  and  Damra,  F.  M.,  1 990,  "On  the  Modal  Stability  of  Imperfect 
Cyclic  Systems",  in  Control  and  Dynamic  Systems:  Advances  in  Theory  and  Application, 
vols  35,  Part  2,  pp  137-164.  Edited  by  C.  T.  Leondes,  Academic  Press. 

Afolabi,  D.,  1990,  On  the  Geometric  Stability  of  Certain  Modes  of  Vibration,  NASA  Technical 
Memorandum  (forthcoming). 

Afolabi,  D.,  1989,  Effects  of  Mistuning  and  Matrix  Structure  on  the  Topology  of  Frequency 
Response  Curves,  NASA  Technical  Memorandum:  TM-102290.  National 
Aeronautics  and  Space  Administration. 

Afolabi,  D.,  1988,  "A  Note  on  the  Rogue  Failure  of  Turbine  Blades",  Journal  of  Sound  and 
Vibration,  vol  122,  pp  535-545. 

Afolabi,  D.,  1988,  "Vibration  Amplitudes  of  Mistuned  Blades",  ASME  Journal  of 
Turbomachinery,  vol  110,  pp  251-257. 

Afolabi,  D.,  1987,  "Linearization  of  the  Quadratic  Eigenvalue  Problem",  Computers  and 
Structures,  vol  25,  pp  1039-1040. 

Afolabi,  D.,  1987,  "The  Dynamic  Stiffness  Method  in  Vibration  Analysis",  International 
Journal  of  Mechanical  Engineering  Education,  vol  15,  pp  133-139. 

Afolabi,  D.,  1986,  "Natural  Frequencies  of  Cantilever  Blades  with  Resilient  Roots",  Journal 
of  Sound  and  Vibration,  vol  110,  pp  429-441. 

Afolabi,  D.,  1985,  "The  Frequency  Response  of  Mistuned  Bladed  Disk  Assemblies",  in 
Vibrations  of  Blades  and  Bladed  Disk  Assemblies,  edited  by  R.  E .  Kielband  N.  F.  Rieger, 
American  Society  of  Mechanical  Engineers,  New  York. 

Afolabi,  D.,  1985,  "The  Eigenvalue  Spectrum  of  a  Mistuned  Bladed  Disk",  in  Vibrations  of 
Blades  and  Bladed  Disk  Assemblies,  edited  by  R.  E .  Kielb  and  N.  F.  Rieger,  American 
Society  of  Mechanical  Engineers,  New  York. 

b.  Articles  In  Conference  Proceedings 

Afolabi,  D.,  1991,  "Modes  and  Quasi-Modes  in  Modal  Analysis",  Proceedings,  9th 
International  Modal  Analysis  Conference,  April  15-18, 1991,  Florence,  Italy. 

Afolabi,  D.,  and  Alabi,  B.,  1991,  "Resolution  of  Double  Modes",  Proceedings,  9th 
International  Modal  Analysis  Conference,  April  15-18, 1991,  Florence,  Italy. 
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Afolabi,  D.,  1991,  "Versal  Deformation  of  Degenerate  Frequency  Response  Functions", 
Proceedings,  4th  Workshop  i  Control  and  Dynamic  Systems,  University  of 
Southern  California,  Los  Angeles,  January,  1991. 

Afolabi,  1991,  "Small  Denominators  and  the  Problem  of  Eigenvector  Stability  in  Structural 
Dynamics"  Midwestern  Mechanics  Conference,  1991. 

Alabi,  B.,  and  Afolabi,  D.,  1990,  "On  the  Problem  of  Uniformly  Moving  Tangential  Loads 
on  a  Semi-infinite  Solid",  Proceedings  of  the  11th  U.S.  National  Congress  of 
Applied  Mechanics,  Tucson,  A Z. 

Happawana,  G.  S.,  Nwokah,  O.  D.  I.,  Afolabi,  D.,  and  Bajaj,  A.  K.,  1990,  "On  the  Dynamics 
of  Perturbed  Symmetric  Systems",  Proceedings  of  the  11th  U.S.  National  Congress 
of  Applied  Mechanics,  Tucson,  A Z. 

Afolabi,  D.,  and  Nwokah,  O.D.1. 1990,  "Effect  of  Mild  Perturbations  on  the  Dynamics  of 
Structures  with  Circulant  Matrices",  Presented  at  the  1990 
AIAA/ASME/ASCE/AHS  Structural  Dynamics  &  Materials  Conference, 
Longbeach,  CA,  April  2-4. 

Afolabi,  D.,  and  Nwokah,  O.  D.  L,  1989,  'The  Frequency  Response  of  Mistuned  Cyclic 
Systems",  Proceedings  of  the  12th  ASME  Biennial  Vibrations  Conference, 
Montreal,  Canada,  September  17-21. 

Afolabi,  D.,  and  Nwokah,  O.  D.  I.,  1989,  "On  the  Modal  Stability  of  Imperfect  Cyclic 
Systems",  Proceedings  of  the  2nd  Workshop  on  Control  Mechanics,  University  of 
Southern  California,  Los  Angeles,  January  25-27. 

Afolabi,  D.,  and  Nwokah,  O.  D.  I.,  1989,  "Dynamics  and  Control  of  Perturbed  Periodic 
Systems",  Proceedings  of  the  21st  Midwestern  Mechanics  Conference,  Michigan 
Technological  University,  August  13-16,  Houghton,  MI.  (Also  published  in 
Developments  in  Mechanics,  vol  15,  pp  539-540. 

Paydar,  N.  and  Afolabi,  D.,  1989,  "Buckling  Analysis  of  Tapered  Sandwich  Plates", 
Proceedings  of  the  21st  Midwestern  Mechanics  Conference,  Michigan 
Technological  University,  Houghton,  MI.  (Also  published  in  Developments  in 
Mechanics,  vol  15,  pp  195-196. 

Afolabi,  D.,  1988,  "Modal  Characteristics  of  Turbomachine  Blades  in  a  Multi-Stage 
Engine",  Proceedings  of  the  6th  International  Modal  Analysis  Conference, 
Orlando,  Florida,  February  1-4,  pp  521-527. 

Afolabi,  D.,  1987,  "An  Anti-Resonance  Technique  for  Detecting  Structural  Damage", 
Proceedings  of  the  5th  International  Modal  Analysis  Conference,  Imperial  College, 
London,  England,  April  u  9,  pp  491-495. 

Afolabi,  D.  and  Paydar,  N.,  1987,  "Receptance  Analysis  of  Timoshenko  Beams  with 
Various  End  Conditions",  Proceedings  of  the  20th  Midwestern  Mechanics 
Conference,  Purdue  University,  Wect  Lafayette,  Indiana,  August  31 -September  2. 
(Also  published  in  Developments  ir  Mechanics,  vol  14c,  pp  1230- 1235. 

c.  Internal  Reports 

Afolabi,  D.,  1988,  "Forced  Response  Analysis  of  Bladed  Disk  Assemblies",  Final  Report 
to:  Allison  Gas  Trubines,  General  Motors  Corp.,  Indianapolis  (20  pages). 

Afolabi,  D.,  1988,  "BDA  —  User  Guide  to  the  Computer  Program  for  Blade  Mistuning 
Analysis",  Contract  Report  to  Gas  Trubines,  General  Motors  Corp.,  Indianapolis 
(31  pages). 
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Afolabi,  D.,  1982,  "Some  Vibration  Characteristics  of  an  Aeroengine  Compressor  Fan", 
Contract  Report  No  8202,  Mechanical  Engineering  Department,  Dynamics  Lab, 
Imperial  College,  London,  England  Report  to:  Rolls-Royce  Research  Center, 
Derby,  England  (23  pages). 

Afolabi,  D.,  1982,  "Effects  of  Random  Mistuning  on  the  Vibration  of  Coupled 
Turbomachine  Blades",  Dynamics  Lab  Report  No  8201,  Mechanical  Engineering 
Department,  Imperial  College,  London,  England  (85  pages). 

Afolabi,  D.,  1982,  "A  Statistical  Interpretation  of  the  Dynamic  Response  Data  of  Bladed 
Disc  Assemblies",  Dynamics  Lab  Report  No  82016,  Mechanical  Engineering 
Department,  Imperial  College,  London,  England  (42  pages). 

Afolabi,  D.  and  Ewins,  D.  J.,  1979,  "A  Bibliography  of  Blade  and  Disk  Vibration", 
Dynamics  Lab  Report  No  7919,  Mechanical  Engineering  Department,  Imperial 
College,  London,  England  (132  pages). 

Afolabi,  D.  and  Ewins,  D.  J.,  1978,  "The  Vibration  of  Turbine  Blade  Packets",  Dynamics 
Lab  Report  No  78007,  Mechanical  Engineering  Department,  Imperial  College, 
London,  England  (156  pages). 

d.  Publications  In  Review 

Afolabi,  D.,  Modal  Interaction  in  Dynamic  Systems,  NASA  Technical  Memorandum 
(forthcoming,  1991). 

e.  Works  In  Progress 

See  attachment. 


3.  Invited  Lectures 

1991  The  Mathematical  Institute  of  Oberwolfach,  West  Germany,  "The  Wey!  Groups  Ak, 
Dk,  Ek  and  the  Random  Vibration  of  Symmetric  Structures"  ;  forthcoming,  1991. 
(Co-author:  S.  H.  Crandall.) 

1 989  NASA  Lewis  Research  Center,  Institute  for  Computational  Mechanics  in  Propulsion, 
July  1989:  "Qualitative  Characteristics  of  Mistimed  Cyclic  Systems". 

1989  Hughes  Aircraft,  Division  of  General  Motors  Corporation,  El  Segundo,  CA,  January 
28, 1989:  "Dynamics  and  Control  of  Perturbed  Periodic  Systems". 

1988  NASA  Lewis  Research  Center,  Workshop  on  Unsteady  Phenomena  in  Turbomachinery, 
Wednesday,  July  20:  "Effect  of  Aerodynamic  and  Structural  Inter-Stage  Coupling 
on  the  Vibration  of  Mistimed  Blades". 

1988  NASA  Lewis  Research  Center,  Structural  Dynamics  Branch  Seminar  Series  June  9: 
"Effect  of  Mistiming  on  Mode  Localization  and  Mode  Splitting  in  a  Multi-Stage 
Engine". 

1988  United  Technologies  Research  Center,  Pratt  &  Whitney  Aircraft,  East  Hartford, 
Connecticut,  March  24:  "A  Systems  Approach  for  Blade  Mistuning". 

1988  Purdue  University,  School  of  Mechanical  Engineering,  West  Lafayette  Indiana, 
September  29:  "Catastrophe  Theory  and  the  Dynamics  of  Imperfect  Cyclic 
Systems". 

1987  Allison  Gas  Turbines  Division,  General  Motors  Corporation,  November  "Some 
Recent  Advances  in  Turbine  Blade  Vibration". 


4.  Consulting  Activities 


1988-present:  NASA  Lewis  Research  Center,  Cleveland,  Ohio. 

1987-89:  Allison  Gas  Turbines  Division,  General  Motors  Corporation,  Indianapolis, 
Indiana. 

1988:  United  Tecnnologies  Research  Center,  East  Hartford,  Connecticut. 


4.  Research  Grants  

a.  Internal  Grants — Principal  Investigator/Grantee:  Dar6  Afolabi 

1989:  Indiana  University,  Bloomington,  Indiana,  Project  Development  Grant. 

1988:  Purdue  University,  West  Lafayette,  Indiana,  XL  Summer  Grant. 

1987:  Indiana  University,  Bloomington,  Indiana,  Overseas  Conference  Grant. 

1987:  Indiana  University,  Bloomington,  Indiana,  Project  Development  Grant. 

1986:  Indiana  University,  Bloomington,  Indiana,  Summer  Faculty  Fellowship. 

b.  External  Grants  Awarded — Principal  Investigator.  Dar6  Afolabi 

1988-90:  from  the  Air  Force  Office  of  Scientific  Research,  AFOSR,  Bolling  Air  Force  Base, 
Washington  DC.  Dynamics  and  Control  ofBladed  Disk  Assemblies  [November,  1988 
-  October,  1989]. 

1987:  from  General  Motors  Corporation,  Allison  Gas  Turbines  Division,  Indianapolis, 
Indiana.  Forced  Response  of  Mistuned  Turbine  Blades  [January,  1988  -  December, 
1988]. 

6.  Review  Activities 

Reviewer  for  the  National  Science  Foundation. 

Reviewer  for  McGraw-Hill  Book  Company. 

Reviewer  for  Journal  of  Sound  and  Vibration. 

Reviewer  for  ASME  Journal  of  Vibration  and  Acoustics. 

Reviewer  for  ASME  Journal  of  Engineering  for  Gas  Turbines  and  Power. 

Reviewer  for  ASME  Journal  of  Turbomachinery. 
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Service  Record 


1.  University  Service 

a.  Campus-wide 

Member,  Faculty  Development  Council  [1989-90]. 

Vice-President,  Black  Faculty  and  Staff  Council  [1988-89]. 

Volunteer,  'Terson-to- Person  Week",  April,  1988. 

Member,  Campus  Interrelations  Sub  Committee  [1987-88]. 

b.  School  of  Engineering  and  Technology 

Member,  Committee  for  the  Development  of  Academic  Computing  Plan,  [1987-88]. 
Member,  Faculty-Student  Advisory  Committee  [1989-90]. 

c.  Department  of  Mechanical  Engineering 

Member,  Search  and  Screen  Committee  [1987-89]. 

Coordinator,  Seminars  and  Colloquia  Series  [1986-198 7]. 

Coordinator,  Laboratory  Equipment  and  Instrumentation  [1988]. 

d.  Student  Support  Service 

Academic  Counselor,  Mechanical  Engineering  Department. 

Advisor,  IUPUI  Chapter,  National  Society  of  Black  Engineers. 

Mentor,  Minority  Engineering  Advancement  Program. 

2.  Public  Service 

Judge,  City  of  Indianapolis  Section,  National  MATHCOUNTS  Competition  for 
Pre-College  students,  February  20, 1988. 

Judge,  ACTSO-87  and  ACTSO-88  Science  and  Engineering  Competitions  for  High  School 
students,  April  1987  and  June  1988 
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